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Abstract 

Wc develop a rigorous theory of non-local Poisson structures, built on 
the notion of a non-local Poisson vertex algebra. As an application, 
we find conditions that guarantee applicability of the Lenard-Magri 
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integrability of various evolution equations, as well as hyperbolic equa- 
tions. 
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1 Introduction 



Local Poisson brackets play a fundamental role in the theory of integrable 
systems. Recall that a local Poisson bracket is defined by (see e.g. [TF86]): 

(1.1) {ui(x),Uj(y)} = H ij (u(y),u'(y), . . . ;d y )5(x - y) , 

where u = (u\, . . . , ug) is a vector valued function on a 1-dimensional man- 
ifold M, 5(x — y) is the 5-function: § M f{y)S(x — y)dy = f(x), and H(d) = 

(-Hjj(5)) i ._ 1 is an t x I matrix differential operator, whose coefficients are 

functions in u, u', . . . , . One requires, in addition, that (1.1) "satisfies 
the Lie algebra axioms". 

One of the ways to formulate the latter condition is as follows. Let 
V be an algebra of differential polynomials in ui,...,v,£, i.e. the algebra 
of polynomials in , i e I = {I, . . . ,£}, n e Z + , with up = m and the 

derivation d, defined by duf 1 ^ = u\ n+1 \ or its algebra of differential functions 
extension. The bracket (1.1) extends, by the Leibniz rule and bilinearity, to 
arbitrary f,geV: 

(1.2) {f(x), g (y)}= 2 2 y[^^;j aray{m(g),ti,-(y)}. 

i,jeIm,neZ+ V u { OU- 

Applying integration by parts, we get the following bracket on V/dV: 

(1.3) m3} ^ S £.H {S} ^, 

where $ is the canonical quotient map V — > V/dV and 5^ is the vector 
of variational derivatives J^- = XineZ (— (?) ra d ( n) . Then one requires that 

the bracket (1.3) satisfies the Lie algebra axioms. (The skewsymmetry of 
this bracket is equivalent to the skewadjointness of H(d), but the Jacobi 
identity is a complicated system of non-linear PDE on its coefficients.) In 
this case the matrix differential operator H(d) is called a Poisson struc- 
ture. (Sometimes in literature, including our previous papers, this is called 
a Hamiltonian structure, or a Hamiltonian operator, but the name Poisson 
structure seems to be more appropriate.) 

Given an element e V/dV, called a Hamiltonian functional, the Hamil- 
tonian equation associated to H{d) is the following evolution equation: 

du ,^d~h 
1.4 — = H(d)—. 

v ' dt K ' 5u 
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For example, taking H (d) = d and h = ^(u 3 + cuu"), we obtain the KdV 
equation: ^ = 3uu' + cu'" . 

Equation (1.4) is called integrable if §h is contained in an infinite di- 
mensional abelian subalgebra A of the Lie algebra V/dV with bracket (1.3). 
Picking a basis {$/i n }neZ+ of A, we obtain a hierarchy of compatible inte- 
grable Hamiltonian equations: 

— = H(d) — , neZ+. 
at n ou 

An alternative approach, proposed in [BDSK09], is to apply the Fourier 
transform F(x, y) *—> \ M dxe^ x ~ y > F(x, y) to both sides of (1.2) to obtain the 
following "Master formula" [DSK06]: 

(1-5) {/a 5 }=2 £ -%(A + 5r^(A + 5)(-A-5r-^. 

du) du- 

For an arbitrary t x i matrix differential operator H(d) this X-bracket is 
polynomial in A, i.e. it takes values in V[d], satisfies the left and right 
Leibniz rules: 

(1.6) {f x gh} = g{f x h} + h{f x g} , {fg x h} = {f x+ ggUh + {fx +8 h}^g , 

where the arrow means that A + d should be moved to the right, and the 
sesquilinearity axioms: 

(1.7) {df x g} = -\{f x g} , {f x dg} = (A + d){f x g) . 

It is proved in [BDSK09] that the requirement that (1.3) satisfies the Lie 
algebra axioms is equivalent to the following two properties of (1.5): 

(1.8) {gxf} = -{f-x-89} , 



(1.9) {fxW) = Mfxh}} + {{fxg}x+^h} . 

A differential algebra V, endowed with a polynomial A-bracket, satisfying 
axioms (1.6)-(1.9), is called a Poisson vertex algebra (PVA). 

It was demonstrated in [BDSK09] that the PVA approach greatly sim- 
plifies the theory of integrable Hamiltonian PDE, based on local Poisson 
brackets. For example, equation (1.4) becomes, in terms of the A-bracket 
associated to H: 
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and the Lie bracket (1.3) becomes 



{S/,fe} = K/A2}| A=0 - 

However the majority of important integrable equations, including the 
non-linear Schroedinger equation, is Hamiltonian with respect to a nondocal 
Poisson bracket. It has been an open problem to develop a rigorous theory 
of such brackets. The purpose of the present paper is to demonstrate that 
the adequate (and in fact indispensable) tool for understanding nondocal 
Poisson brackets is the theory of "nondocal" PVA. 

We define a non-local X-bracket on the differential algebra V to take 
its values in V((A -1 )), formal Laurent series in A~ x with coefficients in V, 
and to satisfy properties (1.6) and (1.7). The main example is the A-bracket 
given by the Master Formula (1.5), where H{d) is a matrix pseudodifferential 
operator. The only problem with this definition is the interpretation of the 
operator ^L_; this is defined by the geometric progression 

2C-i)-A~'a». 

ne£ + 

Property (1.8) of the A-bracket is interpreted in the same way, but the 
interpretation of property (1.9) is more subtle. Indeed, in general, we 
have {f x {9»h}} e V((A- 1 ))(( M " 1 )), but {g„{hh}} e V((^ 1 ))((A- 1 )), and 
{{f\g}\+nh} e V(((A + /x) -1 ))^ -1 )), so that all three terms of (1.9) lie in 
different spaces. Our key idea is to consider the space 

V A)M = V[[A- 1 ,/i- 1 ,(A + /i)- 1 ]][A,/i], 

which is canonically embedded in all three of the above spaces. We say that 
a A-bracket is admissible if 



{fx{9fi.h}} e V\,» for all f,g,heV. 

It is immediate to see that then the other two terms of (1.9) lie in Va,^ as 
well, hence (1.9) is an identity in V\^. 

We call the differential algebra V, endowed with a non-local A-bracket, 
a non-local PVA, if it satisfies (1.8), is admissible, and satisfies (1.9). 

For an arbitrary pseudodifferential operator H{d) the A-bracket (1.5) 
is not admissible, but it is admissible for any rational pseudodifferential 
operator, i.e. such that the entries of the matrix H{d) are contained in the 
subalgebra generated by differential operators and their inverses. We show 
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that, as in the local case (see [BDSK09]), this A-bracket satisfies conditions 
(1.8) and (1.9) if and only if (1.8) holds for any pair Ui,Uj, and (1.9) holds 
for any triple U{, Uj, u^. Also, (1.8) is equivalent to skewadjointness of H{d). 

The simplest example of a non-local PVA corresponds to the skewadjoint 
operator H(d) = d~ l . Then 

{u x u} = A -1 , 

and equation (1.9) trivially holds for the triple u,u,u. Note that (1.1) in 
this case reads: {u(x),u(y)} = dy S(x — y), which is quite difficult to work 
with (cf. [MN01]). 

The next example corresponds to Sokolov's operator [Sok84] H(d) = 
u'd -1 o v! . The corresponding A-bracket is 

{u\u} = u - - - v! . 

The verification of (1.9) for the triple u, u, u is straightforward. 

We say that a rational pseudo differential operator H{d) is a Poisson 
structure on V if the A-bracket (1.5) endows V with a structure of a non- 
local PVA (in other words H(d) should be skewadjoint and (1.9) should hold 
for any triple Ui,Uj,Uk)- 

Fix a "minimal fractional decomposition" H = AB^ 1 . This means that 
A, B are differential operators over V, such that Ker A n Ker B = in any 
algebra of differential functions extension of V. It is shown in [CDSK12b] 
that such a decomposition always exists and that the above property is 
equivalent to the property that any common right factor of A and B is 
invertible over the field of fractions of V. Then the basic notions of the theory 
of integrable systems are defined as follows. A Hamiltonian functional (for 

H = AB^ 1 ) is an element §h e V/dV such that = B(d)F for some 
F e V . Then the element P = A{d)F is called an associated Hamiltonian 

vector field, and we write lh «^-> P or P +±U Jft. Denote by F{H) c V/dV 
the subspace of all Hamiltonian functionals, and by H(H) a V the subspace 
of all Hamiltonian vector fields (they are independent of the choice of the 
minimal fractional decomposition for H): 

j-(tf) = (£) (imfl) cz V/dV , H(H) = A^B-^Imj-^ c V e . 

Then it is easy to show that F{H) is a Lie algebra with respect to the well- 
defined bracket (1.3), and T~L(H) is a subalgebra of the Lie algebra V e with 
bracket [P, Q] = Dq(d)P — Dp(d)Q, where Dp{d) is the Frechet derivative. 
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A Hamiltonian equation, corresponding to the Poisson structure H and 
a Hamiltonian functional e J~{H) , with an associated Hamiltonian vector 
field P e H(H), is the following evolution equation: 

du „ 

(110 > Tt= P - 

Note that (1.10) coincides with (1.4) in the local case. The Hamiltonian 
equation (1.10) is called integrable if there exist linearly independent infinite 
sequences J/i n e T(H) and P n e %{H), n e Z+, such that J/io = Po = P, 
P n is associated to §h n , and {$/i m) S^-n} = 0, [P m ,.P n ,] = for all m,ne Z + . 
In this case we have a hierarchy of compatible integrable equations 

= P n , neZ + . 

dt n 

(The P n 's are called the generalized symmetries of equation (1.10) and the 
h n 's are its conserved densities.) 

Having given rigorous definitions of the basic notions of the theory of 
Hamiltonian equations with non-local Poisson structures, we proceed to es- 
tablish some basic results of the theory. 

The first result is Theorem 5.1, which states that if H and K are com- 
patible non-local Poisson structures and K is invertible (as a pseudodiffer- 
ential operator), then the sequence of rational pseudodifferential operators 

is a compatible family of non-local 
Poisson structures. (As usual [Mag78, Mag80] a collection of non-local Pois- 
son structures is called compatible if any their finite linear combination is 
again a non-local Poisson structure.) This result was first stated in [Mag80] 
and its partial proof was given in [FF81] (of course, without having rigorous 
definitions). 

Next, we give a rigorous definition of a non-local symplectic structure 
and prove (the "well-known" fact) that, if S is invertible (as pseudodiffer- 
ential operator), then it is a non-local symplectic structure if and only if 
S^ 1 is a non-local Poisson structure (Theorem 6.2). Since we completely 
described (local) symplectic structures in [BDSK09], this result provides a 
large collection of non-local Poisson structures. We also establish a connec- 
tion between Dirac structures (see [Dor93] and [BDSK09]) with non-local 
Poisson structures (Theorems 6.12 and 6.17). 

After that we discuss the Lenard-Margi scheme of integrability for a 
pair of compatible non-local Poisson structures H, K, similar to that dis- 
cussed in [Mag78, Mag80, Dor93, BDSK09] in the local case, and give suffi- 
cient conditions when this scheme works (Theorem 7.15 and Corollary 7.16). 
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This means that there exists an infinite sequence of Hamiltonian functionals 
^h n , n e Z+, and Hamiltonian vector fields P n , n e Z + , such that we have 



(1.11) 



So 



H 



P 



K 



n 



Pi 



A 



H 



and the spans of the $/i n 's and of the P n 's are infinite dimensional. This 
produces integrable Hamiltonian equations Ut n = P n - 

Let us also mention that the Lenard-Magri scheme in the weakly non- 
local case (in the sense of [MN01]) was studied in [Wan09]. 

Now we compare briefly our approach to integrability with other alge- 
braic approaches. Probably the earliest approach is the Lax pair presen- 
tation (see the book [Dic03]). The main difficulty of this approach is to 
establish linear independence of the integrals of motion. Another popular 
approach is based on a recursion operator (see books [01v93] and [Bla98]), 
which is applied to a conserved density or a generalized symmetry to pro- 
duce a new one. Unfortunately, since R is non-local (even for the KdV) this 
approach is non rigorous and often leads to wrong conclusions (as demon- 
strated, for example, in [SW01]). A more recent approach, due to Dorfman 
[Dor93] is based on the notion of a Dirac structure. This theory, along with 
its further developments in [BDSK09] and the present paper, is a basis of 
our non-local bi-Hamiltonian approach. In fact, our approach overcomes 
the main difficulty, that of constructing a Dirac structure, in Dorfman's ap- 
proach, and that of proving linear independence of integrals of motion in 
the Lax pair approach. Also, its advantage as compared to the recursion 
operator approach is that it is rigorous. 

The applications of this theory to concrete examples are studied in Sec- 
tions 8, 9 and 10. In Section 8 we consider three compatible scalar non-local 
Poisson structures: 



and take for a compatible pair [H , K) two arbitrary linear combinations of 
these three structures: H = ^ li aiLi, K = ^ji^iLi. We study in detail for 
which values of the coefficients ctj and hi the corresponding Lenard-Magri 
scheme is integrable. 

Furthermore we study when the infinite sequence (1.11) can be extended 
to the left. The most interesting cases are those when the sequence is 
"blocked" at some step P- n , n > 0, to the left. This leads to some interest- 
ing integrable hyperbolic equations. As a result, we prove integrability of 
two such equations 



Lx = d, L 2 = d- 1 , L 3 = u'd~ l o v! ( Sokolov [Sok84]) 



(1.12) 
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where a and e are or 1, and 

(1.13) u tx = u + (u 3 ) xx . 

Of course, in the case when e = 0, equation (1.12) is the Liouville (respec- 
tively sinh-Gordon) equation if a = (resp. a = 1). For e = 1 equation 
(1.12) was studied in [Fok95]. Equation (1.13), studied in [SW02], is called 
the short pulse equation. Its integrability was proved in [SS04] 

In Section 9 we study, in a similar way, two linear combinations of the 
compatible non-local Poisson structures 

L\ = u'd^ 1 on', L2 = d^ 1 o u'd~ o u'd^ 1 ( Dorfman [Dor93]) . 

As a result we (re) prove integrability of the Schwarz KdV (also called de- 
generate Krichever-Novikov) equation 

2 U x 

and also, moving to the left, establish integrability of the following equation 

( 1 fU tx \ \ 1 / Cp + CjU + c 2 u 2 \ 

\U X \U X ) x) x U x \ U x Jx' 

where Qq,Ci,C2 are arbitrary constants. 

Finally, in Section 10 we study, in a similar way, three two-component 
non-local Poisson structures that are used in the study of the non-linear 
Schroedinger equation (NLS), see [Mag80, TF86, Dor93, BDSK09]. As a 
result, we establish integrability of the following generalization of NLS: 

iipt = i>xx + aip\tp\ 2 + if3(ip\ip\ 2 ) x , 

where a and f3 are arbitrary constants (NLS corresponds to ft = 0). This 
equation has been studied in the papers [CLL79] and [WKI79] (see also 
[KN78] and [CC87]). 

In conclusion of the introduction we would like to comment on our defi- 
nition of integrability. The existence of infinitely many linearly independent 
integrals of motion in involution J/t n , and of infinitely many linearly inde- 
pendent commuting higher symmetries P n , is only a necessary condition of 
integrability. In Section 7.3 we introduce the notion of complete integrabil- 
ity which, in our opinion, is the right necessary and sufficient condition of 
integrability. This condition requires that the orthocomplement to the span 
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H of the variational derivatives of the conserved densities £ n = , n e Z + , 
lies in the span II of the commuting generalized symmetries P n , n e Z+, 
and the orthocomplement to II lies in S. This definition is a straightforward 
generalization of Liouville integrability of finite dimensional Hamiltonian 
systems. We intend to study this notion in a forthcoming publication. 

Throughout the paper, unless otherwise specified, all vector spaces are 
considered over a field F of characteristic zero. 

We wish to thank Pavel Etingof and Andrea Maffei for (always) useful 
discussions. We are greatly indebted to Alexander Mikhailov and Vladimir 
Sokolov for very useful correspondence and discussions. We also wish to 
thank Takayuki Tsuchida for pointing out, right after the paper appeared 
in the archives, various references where some of the equations that we 
consider were previously studied. The present paper was partially written 
during the first author's several visits to MIT, the second author's several 
visits to the Center for Mathematics and Theoretical Physics (CMTP) in 
Rome, and both authors' visits to IHP and IHES, which we thank for their 
warm hospitality. 

2 Rational matrix pseudodifferential operators 

2.1 The space V x ^ 

Throughout the paper we shall use the following standard notation. Given 
a vector space V, we denote by V[A] the space of polynomials in A with 
coefficients in V, by ^[[A -1 ]] the space of formal power series in A -1 with 
coefficients in V, and by F((A -1 )) = V[[A -1 ]][A] the space of formal Laurent 
series in A -1 with coefficients in V. 

We have the obvious identifications V[A, fj] = V[\][fj] = V[n][X] and 
^[[A -1 ,/* -1 ]] = ^[[A- 1 ]]^- 1 ]] = VR/i" 1 ]]^- 1 ]]. However the space 
V((A _1 ))((// _1 ) does not coincide with V((/x ))((A )). Both spaces con- 
tain naturally the subspace T^[[A , /i ]][A,/i]. In fact, this subspace is 
their intersection in the ambient space ^[[A- 1 ,^- 1 ]] consisting of all series 
of the form 2 m ,neZ a m , n \ m fi n . 

The most important for this paper will be the space 

Vx, li :=V[[\- 1 ,n- 1 ,(\ + /M)- 1 ]][\,n] t 

namely, the quotient of the F[A, //, z/]-module F[[A _1 , fi^ 1 , z^ _1 ]][A, /i, v\ by 
the submodule [y— A— /u)y[[A _1 , fi^ 1 , i^ _1 ]][A, /i, v\. By definition, the space 
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V\ tfJi consists of elements which can be written (NOT uniquely) in the form 

MNP 

(2.1) ^=222 a m , n , p A m /i"(A + ^, 

m=— oo n=— oo p=— oo 

for some M, N,P e Z (in fact, W6 can always choose P ^ 0), and. Q"m,n,p ^ v. 

In the space y[[A _1 , /U , f _1 ]][A, fi, u] we have a natural notion of de- 
gree, by letting deg(A) = deg(/x) = deg(v) = 1. Every element A e 
^[[A -1 ,^ 1 ,^ -1 ]]^,^, i/] decomposes as a sum 

^ = £1-00^ (possibly 
infinite), where is a finite linear combination of monomials of degree 
d. Since v — A — fj, is homogenous (of degree 1), this induces a well-defined 
notion of degree on the quotient space V\ tlJr , and we denote by , for d e Z, 
the span of elements of degree ci in V\ iM . If A e Va iM has the form (2.1), then 
it decomposes as A = Xidi 4 ^o +P A^ d \ where A^ e is given by 

A^= 2 a m , n , p A"V(A + / <. 

m^M,n^N,p^P 
(m+n+p=d) 

The coefficients a m ^ n)P e V are still not uniquely defined, but now the sum 
in A" is finite (since ci — 2A sC m,n,p < A' := max(M, TV, P)). Hence, we 
have the following equality 

Fi = y[A ±1 , / , ±1 ,(A + ,.)- 1 ]< i , 

where, as before, the superscript d denotes the subspace consisting of poly- 
nomials in A- 1 ,//- 1 , (A + /x) _1 , of degree d. 

Lemma 2.1. The following is a basis of the space V d ^ over V : 

A^V, ieZ ; \ d+i {\ + (j,)-*, i e Z >0 = {1, 2, . . . } , 

in the sense that any element of the space V d ^ can be written uniquely as a 
finite linear combination of the above elements with coefficients in V . 

Proof. First, it suffices to prove the claim for d = 0. In this case, letting 
t = ///A, the elements of V® are rational functions in t with poles at 
and -1. But any such rational functions can be uniquely written, by partial 
fractions decomposition, as a linear combination of t 1 , with i e Z, and of 
(1 + t)-\ with i e Z >0 . □ 



12 



Remark 2.2. One has natural embeddings of V\ ^ in all the vector spaces 
n(A- 1 ))((^ 1 )), y((^ 1 ))((A- 1 )), y((A- 1 ))(((A + M )- 1 )), y(( M -i))(((A + 
/x)- 1 )), y(((A + ^)- 1 ))((A- 1 )), y(((A + ^)- 1 ))(( / u- 1 )), defined by expanding 
one of the variables A, /i or v = A + [i in terms of the other two. For example, 
we have the embedding 

(2-2) v,A : Vx, M -^((A- 1 ))(( / u- 1 )), 

obtained by expanding all negative powers of A + \i in the region > |A|: 

(2-3) V ,a(A + M)-"- 1 = 2 (~ V X ) Ay"""*" 1 • 

Similarly in all other cases. Note that, even though „ is naturally embed- 
ded in both spaces V((A _1 ))(QLi -1 )) and V((// _1 ))((A -1 )), it is not contained 
in their intersection V[[A -1 ,// -1 ]][A,/i/]. 

Lemma 2.3. IfV is an algebra, then is also an algebra, with the obvious 
product. Namely, if A(X, /x),J5(A, /x) e V\ ifl) then A(X, fi)B(X, [/,) e V\^. 
More generally, if S,T : V —* V are endomorphisms of V ( viewed as a 
vector space), then 

A(A + S,^ + T)B(V)ey Ai(1 , 

where we expand the negative powers of A + S and fi + T in non-negative 
powers of S and T, acting on the coefficients of B. 

Proof. We expand A and B as in (2.1): 

MNP 

A(X tf i)= 2 2 S « TO A"V(A + /^, 

m=— oo n=— QO p=— oo 
M' N' P' 

B(\,fj,)= ^ 2 2 6 m>V A™y'(A + /z)P'. 

m'=— oo n'=— oo p'=— oO' 
Using binomial expansion, we then get 

M+M' N+N' P+P' 
A(X + S,fi + T)B(X^)= 2 2 2 c w A"V(A + /iF, 

fh——ao n=— oo p=— 00 

where 

Cm,n,p =2 2 2 

m^M,m'^M',i>0 n^N,n'^N',j^Op^P,p'^P',k>0 
(m+m'—i=m) (n+n'—j=n) (p+p 1 —k—p) 

■ II • I I ; I ®"m,n,p 

(S l Ti(S + T) k b m ,y, p/ ). 
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To conclude, we just observe that each sum in the RHS is finite, since, for 



example, in the first sum we have i 
fh- M < m' < M' . 



m + m 



m, m — M' ^ m ^ M and 

□ 



Lemma 2.4. Let V be a vector space and let U be a subspace. Then 
we have: 

{ie^l^ef/ftA- 1 ))^- 1 ))} 



{^4 e V x ,, 
{A e V x ,„ 



c^AeU((^))((\-i))} 
L X+flA AeU((\- 1 ))(((X + ^- 1 ))} 
iA+w ief/(( /( - 1 ))(((A + / .)- 1 ))} 
, A , A+ ^eC/(((A + / u)- 1 ))((A- 1 ))} 



Proof. We only need to prove that {A e V x ^ \ l^ x A e U {{\~ l )){{^r 1 ))} c 
U x a- Indeed, the opposite inclusion is obvious, and the argument for the 
other equalities is the same. 

Let A e be such that its expansion i^ x A e V^((A _1 ))((/i -1 )) has 
coefficients in U . We want to prove that A lies in U x „. By Lemma 2.1 A 
can be written uniquely as 

N N 

A= ^ ^A d+ V 4 + 2^A d ^(A + / ur J , 

i=-M j=l 

with e V. Its expansion in l/((A _1 ))((/i~ 1 )) is 



.4 



2 ^+V< + 2 2 

=-Af j=ifc=o 



w 



+k —j—k 



Since, by assumption, l^ x A e U((\ ))((/U )), we have 



N 



6 C/ 



for - M < i < -1 
for s^i ^ N , 

for i > N . 



From the first condition above we have that v% lies in U for i < 0. Using the 
third condition, we want to deduce that Wj lies in U for all 1 < j < JV. It 
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then follows, from the second condition, that Vi lies in U for i > as well, 
proving the claim. 

For i > N and 1 < j s$ N we have Qj*.) = Hence, we 

will be able to deduce that lies in U for every j, once we prove that the 
following matrix 

has rank TV. First, the sign (— does not change the rank of the above 
matrix. So it sufficies to prove that the matrices 

T =(( l ~ 1 )) 

V V 7 — 1 / / h+l^i^/i+iV ' 

V V 77 l^j^N 

are non-degenerate for every h ^ 0. This is clear since the matrix Tq is 
upper triangular with ones on the diagonal, and, by the Tartaglia-Pascal 
triangle, and T^+i have the same determinant. □ 



2.2 Rational pseudodifferential operators 

For the rest of this section, let A be a differential algebra, i.e. a unital 
commutative associative algebra with a derivation d, and assume that A is 
a domain. For a e A, we denote a' = 5(a) and a" = 5 n (a), for a non 
negative integer n. We denote by fC the field of fractions of A. Then of 
course we can extend d to a derivation of JC making it a differential field. 

Recall that a pseudodifferential operator over A is an expression of the 
form 

N 

(2.4) A = A(8)= ^ a nd n , a n eA. 

If otv 0, one says that A has order N . Pseudodifferential operators form 
a unital associative algebra, denoted by A^d" 1 )), with product o defined 
by letting 

(2.5) d n o a = J] (" J a^d n ~ j , neZ,aeA. 
We will often omit o if no confusion may arise. 
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Clearly, K{{d~ 1 )) is a skewfleld, and it is the skewfield of fractions of 
^4((5 -1 )). If A e A((d~ 1 )) is a non-zero pseudodifferential operator of order 
N as in (2.4), its inverse A^ 1 e £((<? )) is computed as follows. We write 

-i 

A = a N (l + J] a^a n+N d n \ d N , 

n=— oo 

and expanding by geometric progression, we get 

00 —1 r 

(2.6) A - l = d- N oY, (- 2 a- 1 ^^) oa, 1 , 

fc=0 n=-oo 

which is well defined as a pseudodifferential operator in /C((<9 -1 )), since, by 
formula (2.5), the powers of d are bounded above by —N, and the coefficient 
of each power of d is a finite sum. 

The symbol of the pseudodifferential operator A{d) in (2.4) is the formal 
Laurent series A{\) = Xin^-oo a «^" e "4((^ _1 ))> where A is an indetermi- 
nate commuting with A. We thus get a bijective map A((d~ 1 )) — > .4((A _1 )) 
(which is not an algebra homomorphism). A closed formula for the asso- 
ciative product in A^d^ 1 )) in terms of the corresponding symbols is the 
following: 

(2.7) (AoB)(X) = A{X + d)B(X). 

Here and further on, we always expand an expression as (A + d) n , n e Z, in 
non-negative powers of d: 

GO / \ 

(2.8) (\ + d) n = 2 r\\ n ~ j y. 

3=0 \ 3 ' 

Therefore, the RHS of (2.7) means Tl,n=-^=o §)"m$ ] X m+n ^. 

The algebra A((d~ 1 )) contains the algebra of differential operators A[d] 
as a subalgebra. 

Definition 2.5. The field KL{d) of rational pseudodifferential operators is 
the smallest subskewfield of /C((5~ 1 )) containing A[d]. We denote A(d) = 
JC(d) n^4((5 -1 )), the subalgebra of rational pseudodifferential operators with 
coefficients in A. 

The following Proposition (see [CDSK12a, Prop. 3. 4]) describes explicitly 
the skewfield KL{d) of rational pseudodifferential operators. 



16 



Proposition 2.6. Let A be a differential domain, and let K. be its field of 
fractions. 

(a) Every rational pseudodifferential operator L e lC{d) can be written as 
a right (resp. left) fraction L = AS~ l (resp. L = S~ 1 A) for some 
A,Se A[d] with 5^0. 

(b) Let L = AS' 1 (resp. L = S^A), with A, S e A[d], S * 0, be a 
decomposition of L e IC(d) such that S has minimal possible order. Then 
any other decomposition L = AiS^ 1 (resp. L = S^Ai), with A\,S\ e 
A[d], we have A x = AK, Si = SK (resp. A x = KA, A x = KS), for 
some K e K\d\. 

2.3 Rational matrix pseudodifferential operators 

Definition 2.7. A matrix pseudodifferential operator A e Mat^ A^d^ 1 )) 
is called rational with coefficients in A if its entries are rational pseudodiffer- 
ential operators with coefficients in A. In other words, the algebra of rational 
matrix pseudodifferential operators with coefficients in A is Ma,t£ X iA(d). 

Let M = (AijB^j 1 ^) . j be a rational matrix pseudodifferential operator 
with coefficients in A, with Aij,Bij e A[d\. By the Ore condition (see e.g. 
[CDSK12a]), we can find a common right multiple B e A[d] of all operators 
Bij, i.e. for every i,j we can factor B = B-ijCij for some Cij e A[d]. 
Hence, AijB^ 1 = AijB -1 , where A\j = AijCij. Then, the matrix M can be 
represented as a ratio of two matrices: M = A(Bl)- 1 . Hence, 

A e Mat e x eA[d], B e A[d], \ 
AvB-teAdd-^ViJ ] ■ 

However, in general this is not a representation of the rational matrix M in 
"minimal terms" (see Definition 2.12 below). 

We recall now some linear algebra over the skewfield /C((5 -1 )) and, in 
particular, the notion of the Dieudonne determinant (see [Art57] for an 
overview over an arbitrary skewfield). 

An elementary row operation of an I x I matrix pseudodifferential oper- 
ator A is either a permutation of two rows of it, or the operation T(i, j; P), 
where 1 ^ i j ^ m and P e /C((5 -1 )), which replaces the j-th row by 
itself minus i-th. row multiplied on the left by P. Using the usual Gauss 
elimination, we can get the (well known) analogues of standard linear al- 
gebra theorems for matrix pseudodifferential operators. In particular, any 



Mat exiA(d) = \ A{B1) 
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matrix pseudodifferential operator A e Mat mX £ ^((d" 1 )) can be brought by 
elementary row operations to a row echelon form. 

The Dieudonne determinant of a A e Mat^ X £ has the form 

detA = c£ d , where ce A, £ is an indeterminate, and deZ, It is defined by 
the following properties: detA changes sign if we permute two rows of A, 
and it is unchanged under any elementary row operation T(i,j;P) defined 
above, for aribtrary i ^ j and a pseudodifferential operator P e /C((5 -1 )); 
moreover, if A is upper triangular, with diagonal entries An of order rii and 
leading coefficoent aj, then 



It was proved in [Die43] (for any skewfield) that the Dieudonne determi- 
nant is well defined and det(AB) = (det >l)(det B) for every t x t matrix 
pseudodifferential operators A,Be Mat^ X £ /C((<5 -1 )). 

The Dieudonne determinant gives a way to characterize invertible matrix 
pseudodifferential operators, thanks to the following well known fact (see e.g. 
[CDSK12a, Prop.4.3]): 

Proposition 2.8. Let T> be a subskewfield of the skewfield K^d^ 1 )), and let 
A eMat^x^D. Then A is invertible in Mat^X* if and only if det A ^ 0. 

Corollary 2.9. Let A e Mat£ X £ /C((5 -1 )) be a matrix with det A 0. Then 
A is a rational matrix if and only if A v is. 

Proof. It is a special case of Proposition 2.8 when T> is the subskewfield 
IC(d) cz /C((5 -1 )) of rational pseudodifferential operators. □ 

Remark 2.10. It is proved in [CDSK12a] that, if A e Mat^ x ^((5~ 1 )) then 
we have det^4 = c£ rf , with c e A, where A is the integral closure of A. 
Furthermore, if c is an invertible element of A, then the inverse matrix A^ 1 
lies in MaW Add' 1 )). 

Definition 2.11. Let A be a differential domain. An I x ^-matrix pseudod- 
ifferential operator A e Mat^ X £ A{[d~ v )) is called non- degenerate if it has 
non-zero Dieudonne determinant, or, equivalently, if it is invertible in the 
ring Matg X £ /C((5 -1 )) of pseudodifferential operators with coefficients in the 
differential field of fractions /C of A. 

Definition 2.12 (see [CDSK12b]). Let H e Ma,t exe lC(d) be a rational ma- 
trix pseudodifferential operator with coefficients in the differential field /C. 
A fractional decomposition H = AB^ 1 , with A,Be Mat£ X £/C[d] and B 




18 



non-degenerate, is called minimal if deg^ det B is minimal (recall that it is 
a non-negative integer). 

Proposition 2.13 ([CDSK12b]). (a) A fractional decomposition H=AB~ l 
of a rational matrix pseudodifferential operator H e Mat^ JC(d) is min- 
imal if and only if 

(2.9) Ker^nKerS = 0, 

in any differential field extension of K. 

(b) The minimal fractional decomposition of H exists and is unique up to 
multiplication of A and B on the right by a matrix differential operator 
D which is invertible in the algebra M&t£ X £ fC[d]. Any other fractional 
decomposition of H is obtained by multiplying A and B on the right by 
a non- degenerate matrix differential operator. 

Remark 2.14. Let A be a differential domain, and let /C be its field of frac- 
tions. A fractional decomposition H = AB~ l of H e Mat£ X £ A[d] over /C 
can be turned into a fractional decomposition over A by clearing the denom- 
inators of A and B. Hence, a minimal fractional decomposition H = AB^ 1 
over A, in the sense that it has minimal possible deg^ det B among all frac- 
tional decompositions of H with A,B e Mat^ X £ A[d], is also minimal over 
JC. 

3 Non-local Poisson vertex algebras 

3.1 Non-local A-brackets and non-local Lie conformal alge- 
bras 

Let R be a module over the algebra of polynomials F[5]. 

Definition 3.1. A non-local X-bracket on R is a linear map {• \ ■} : R0R — ► 

i?((A -1 )) satisfying the following sesquilinearity conditions (a, be R): 

(3.1) {da x b} = -\{a x b} , {a x db} = (A + d){a x b} . 

The non-local A-bracket {-a'} is said to be skewsymmetric (respectively 
symmetric) if (a, b e R) 

(3.2) {b x a} = -{a-x-db} ( resp. = {a_x-db}) ■ 
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The RHS of the skewsymmetry condition should be interpreted as fol- 
lows: if {a x b} = Y>n=-co c ™^"> tnen 

N N co , n. 

{a_ A _ 5 6}= 2 { - X -d) n c n = 2 ^ ( J(-l) n (^„)A"- fc 

n=— ao n=— oo fc=0 ^ ^ 

m=-oo fc=0 ^ ' 

In other words, we move —A — d to the left and we expand in non negative 
powers of d as in (2.8). 

In general we have {a\{b^c}} e i?((A _1 ))((/.f -1 )) for an arbitrary A- 
bracket {•>•}• Recall from Section 2.1 that Rx,n can be considered as a 
subspace of i?((A~ 1 ))((^- 1 )) via the embedding l^^X' 

Definition 3.2. The non-local A-bracket {■ a '} on i? is called admissible if 

{«a{^c}} e R X,n Va, 6, cei?. 
Remark 3.3. If {-a - } is a skewsymmetric admissible A-bracket on R, then 
{b[i{a\c}} e i?A,/i an d {{ fl A^}A+/t c } e -^A,^ for all a, 6, c e R. Indeed, the 
first claim is obvious since Rx,^ = R^ t x- For the second claim, by skewsym- 
metry {{axfyx+^c} = — {c-A-/_t-a{ a A^}}i an d by the admissibility assump- 
tion {c u {axb}} e Rx,v To conclude it suffices to note that when replacing 
v by —A — fj, — d in an element of R\ v = R[[X ,u ,(X + f ) _1 ]][A, z^], 
we have that v~ l is expanded in negative powers of A + /i and (A + 
is expanded in negative powers of fx. As a result, we get an element of 
R[[X-\ f i-\(X + ^- 1 ]][X, / j] = Rx,^. 

Definition 3.4. A non-local Lie conformal algebra is an F[d]-module R 
endowed with an admissible skewsymmetric A-bracket { • a • } : R® R —* 
-R((A -1 )) satisfying the Jacobi identity (in Rx,^)'- 

(3.3) {a A {6 M c}} - {b^axc}} = {{axb} x +^c} for every a,b,ce R. 

Example 3.5. Let R = (F[<9] <g> V) FC, where V is a vector space with 
a symmetric bilinear form (• | •). Define the (non-local) A-bracket on R by 
letting C be a central element, defining 

{a x b} = (a^CA" 1 {ora,beV, 

and extending it to a A-bracket on R by sesquilinearity. Skewsymmetry for 
this A-bracket holds since, by assumption, (• | •) is symmetric. Moreover, 
since any triple A-bracket is zero, the A-bracket is obviously admissible and 
it satisfies the Jacobi identity. Hence, we have a non-local Lie conformal 
algebra. 
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3.2 Non-local Poisson vertex algebras 

Let V be a differential algebra, i.e. a unital commutative associative algebra 
with a derivation d : V — » V. As before, we assume that V is a domain and 
denote by /C its field of fractions. 

Definition 3.6. (a) A non-local X-bracket on the differential algebra V is 
a linear map {-a - } : V ® V — > V((A -1 )) satisfying the sesquilinearity 
conditions (3.1) and the following left and right Leibniz rules: 

, 4 s {aA&c} = 6{aAc} + c{a A ^} , 

{ab x c} = {a x+ dc}^b+ [b x+d c}^a . 

Here and further an expression {a\ + gb}^c is interpreted as follows: if 
{a\b} = S n= _ Q0 c n X n , then {a A +5^}^c = S n= _ Q0 c n (A + <9) n c, where we 
expand (A + d) n c in non-negative powers of d as in (2.8). 

(b) The conditions of (skew) symmetry, admissibility and Jacobi identity for 
a non-local A-bracket {• a •} on V are the same as in Definitions 3.1, 3.2 
and 3.4 respectively. 

(c) A non-local Poisson vertex algebra is a differential algebra V endowed 
with a non-local Poisson X-bracket, i.e. a skewsymetric admissible non- 
local A-bracket, satisfying the Jacobi identity. 

Example 3.7 (cf. Example 3.5). Let V = F[u| \i = \, ...,£, n e Z + ] be 
the algebra of diffenertial polynoamials in £ differential variables u±, . . . ,U£. 
Let C = [c%j) i =1 be an £ x £ symmetric matrix with coefficients in F. The 
following formula defines a structure of a non-local Poisson vertex algebra 
on V: 

m,neZ + i,jeZ+ C^j 

For example, {ui\Uj} = CyA -1 but, for P,Q e F[ui, . . . , u^] a V, we get an 
infinite formal Laurent series in A -1 : 

lj=l J 

!j = ln=0 J 
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We will prove that this is indeed a non-local Poisson A-bracket in the next 
section, where we will discuss a general construction of non-local Poisson 
vertex algebras, which will include this example as a special case (see The- 
orem 4.8). 

Proposition 3.8. Let {■ \ -} be a non-local Poisson vertex algebra structure 
on the differential domain V. Then there is a unique way to extend it to a 
non-local Poisson vertex algebra structure on the differential field of fractions 
/C, and it can be computed using the following formulas (a,b e K\{0}): 

(3.5) {axb- 1 } = -b- 2 {a x b} , {a7» = -{ax+gb}^ 2 . 

Proof. It is straightforward to check that formulas (3.5) define a non-local 
A-bracket on the field of fraction /C, satysfying all the axioms of non-local 
Poisson vertex algebra. In particular, admissibility of the A-bracket can be 
derived from Lemma 2.3. The details of the proof are left to the reader. □ 

Thanks to Proposition 3.8 we can extend, uniquely, a non-local Poisson 
vertex algebra A-bracket on V to its field of fractions /C. The following results 
are useful to prove admissibility of a non-local A-bracket. 

Lemma 3.9. Let V be a differential algebra, endowed with a non-local A- 
bracket {■ \ }. Assume that V is a domain, and let K. be its field of fractions. 
Let S = {Sij) i: j eI £ Mat^ X £ (/C((5 -1 ))) be an invertible £ x t matrix pseu- 

dodifferential operator with coefficients in K. Letting Sij = 2^=-oo s ij;nd n , 
the following identities hold for every a e /C and i,j e /; 



(3.6) 



and 



(3.7) 



t N 
r,t=l n=— oo 

{a x s rt , n }(fi + dT(S~X(fi) e /C((A- 1 ))(( / u- 1 )) , 



N 



{(S 1 )ii(A) A+M a} = - 2 £ i S rt;n X+fl+d a}^ 
r,t=l n=—x> 

(A + 5)"(5- 1 ) tJ (A)), A , A+M (5*- 1 )„(^) e /C(((A + ^)" 1 ))((A- 1 )) 



where t M)A and /-a,a+^ are as in (2.3). In equation (3.7) S* denotes the 
adjoint of the matrix differential operator S (its inverse being (S^ 1 )*). 
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Proof. The identity S o S = 1 becomes, in terms of symbols, 



I 

^ S r j(fl + d)(S' 1 ) t j(fJ.) = S rj . 

t=x 

Taking A-bracket with a, we have, by sesquilinearity and the (left) Leibniz 
rule, 

e 

= 2{a A 5 rt (/i + d)(5-%-(/x)} 
t=l 

t N 

t=l n=— oo 
I N 

t=l n=-oo 
i 

+ 2 V,A<S ri (A + fi + d){a x {S- 1 ) tj {fi)} . 
t=i 

Note that l^.\S(\ + \i + d) is invertible in Mat^ (£[<5]((A _1 ))((/i -1 ))), its 
inverse being (A + pi, + d). We then apply L flj \{S^ 1 )i r (X + /i + d) on 

the left to both sides of the above equation and we sum over r = 1, ... ,£, 
to get 

i 

t=l 

l l N 
r=l f=l n=— oo 



proving equation (3.6). 
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Similarly, for the second equation we have, by the right Leibniz rule, 



e 

t=i 

l N 

t=l n=— go 
I N 

(=1 n=— go 
£ AT 

+ 1j 1j {( 5 ~%'( A )a+m+5«}^a,a+ m (A - A - n - d) n s rt - n 

t=l 71= — GO 

£ iV 

t=l 71= — GO 
I 

+ Zj {( 5 ' _1 )y(A)A+ M +5a}^'-A,A+ M 5'* r (//) . 
t=i 

We next replace in the above equation /x (placed at the right) by /i + 5, and 
we apply the resulting differential operator to i\ t \+ il {S*~ 1 ) r i(fi) . As a result 
we get, after summing over r = 1, . . . , £, 

e 

Yi {( 5_1 )^( A )A+M+a«}^4i 

l l N 

= "S £ E {^;n A+M+(5 a}^((A + 5)"(5- 1 )^(A)), A , A+M (5*- 1 )„( / ,), 

r=l t=l n=0 



proving equation (3.7). □ 

Corollary 3.10. Let V be a differential algebra, endowed with a non-local A- 
bracket {■ A •}. Assume that V is a domain, and let fC be its field of fractions. 
Let S = (SijJijzi E M.a.te x e (/C[<5]) be non- degenerate (cf. Definition 2.11). 
Then the following identities hold for every a e /C and i,j e /: 
(3.8) 

{axiS- 1 )^)} 

l N 

= " 2 £(Oir(A + ii + d){a x s rt]n }(fi + e/C AiM , 

r,t=l n.=0 
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and 
(3.9) 

{(,S- 1 ) iJ (A) A+At a} 

i N 

r,t=l n=0 
where Sij = 2n=0 s ij;nd n . 

Proof. It is immediate from equations (3.6) and (3.7). □ 

Corollary 3.11. Let V be a differential algebra, endowed with a non-local X- 
bracket {■ x ■}. Assume that V is a domain, and let /C 6e its /zeZd of fractions. 
Let A e V(d) = /C(<9)nV((<9 -1 )) be a rational pseudodifferential operator with 
coefficients in V . Then {a x A([i)} and {A(X) x+IM a} lie in V Aj(U for every a e V. 
/n particular, the X-bracket is admissible. 

Proof. First, note that if the pseudodifferential operators A, B e K,((d~ 1 )) 
satisfy the conditions 

{a x A{fi)} , {A{X) x+fx a}{a x B{fi)} , {B(X)\+ fl a} e K x ,» , 

for every a e /C, so does Ai?. Indeed, by the Leibniz rule, 

{a x (AB)(n)} = {a x A(fi + d)B(n)} 

= {a x A(fi + d)}^B(ix) + A(X + f i + d){a x B([i)} , 

and both terms in the RHS lie in /C Ai/i by the assumption on A and B, 
thanks to Lemma 2.3. Similarly, by the right Leibniz rule, 

{{AB)(X) x+ ^a} = {A(X + d)B(X) x+tl a} 

= {B(X) x+tl+8 a}^L XMti A*(n) + {A{X + d) x+ ^ +d a}^B{X) , 

and both terms in the RHS lie in K\n, (rather in the image of fC x ^ in 
/C(((A + /i) _1 ))((A -1 )) via t A ,A+ M ) by Lemma 2.3. By Corollary 3.10 we 
have that, if 5 e V[d], then {a x S~ 1 ([i)} and {S~ 1 (X) x+fJ ,a} lie in K X;fl for 
all a e /C. Hence, by Definition 2.5 and the above observations, we get 
that, if A e V{d) = K(d) n V^ 1 )), then {a A A(/z)} and {^(A) A+M a} lie in 
/C AjAt for all a e /C. On the other hand, if a e V, we clearly have {a A .A(//)} e 
V((A" 1 ))((/i- 1 )) and {A(X) x+fl a} e V(((A + ^)" 1 ))((A- 1 )). The claim follows 
from Lemma 2.4 applied to V = fC and U = V. □ 

Remark 3.12. In the case when S e V(d) is a rational pseudodifferential 
operator with coefficients in V, thanks to Corollary 3.11, we can drop t„ A 
and i AjA +n respectively from equations (3.6) and (3.7), which hold in the 
space V Xjil . 
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4 Non-local Poisson structures 



4.1 Algebras of differential functions 

Let Re = F[u^ | i e I, n e Z+] be the algebra of differential polynomials 
in the I variables U{, i e / = {1,...,£}, with the derivation d defined by 
d(Uj"^) = ttj n+ . The partial derivatives — are commuting derivations of 
Ri, and they satisfy the following commutation relations with d: 



(4.1) 



d 8 



du<; n) 



3 

(the RHS is if n = 0) 



Recall from [BDSK09] that an algebra of differential functions is a dif- 
ferential algebra extension V of Re, endowed with commuting derivations 

extending the usual partial derivatives on R%, such that only a finite number 
of d/ are non-zero for each / e V, and such that the commutation rules 

(4.1) hold on V. 

Given an algebra of differential functions V, its localization by a multi- 
plicative subset is again an algebra of differential functions. Also, we can 
add to V solutions of algebraic equations over V, or functions of the form 
F(ipi, . . . , (p n ), where F is an infinitely differentiable function in n variables 
and (pi, . . . , cp n lie in V, to obtain again an algebra of differential functions. 
(Note, though, that in general we cannot add to V solutions of linear differ- 
ential equations. For example, a solution of the equation /' = fv! is / = e u , 
which can be added, while a non-zero solution of the equation /' = fu can 
never be added, due to simple differential order considerations.) We will use 
these facts in the examples further on. 

Remark 4.1. An algebra of differential functions V can be equivalently de- 
fined as follows. It is a commutative associative algebra extension of Re, 
endowed with commuting derivations ^ and ^-f^y, i e I,n e Z + , such that 

d 



du 



fa acts trivially on Re, - ^ extends the usual action of partial derivatives 
on Re, and for every / e V we have S f n) = for all but finitely many choices 

of indices (i, n). In this case, the total derivative d : V — > V defined by the 
formula 



df Jn+l) _j_ df_ 



gjn) Ui " dx ' 



2(3 



satisfies equation (4.1) 

Note that if the algebra of differential functions V is a domain, then its 
field of fractions K, is again an algebra of differential functions in the same 
variables u\ , . . . , U£, with the maps ° (n) : /C — ► K. defined in the obvious 

8u\ 

way. When V = K, we call it a field of differential functions. 

We denote by C = {ceV|5c = 0} V the subalgebra of constants, and 

by 

J = |/eV — — 7-r = for all i e J, n e Z+j c V 

the subalgebra of quasiconstants. It is easy to see that C a T . 

Given / e V which is not a quasiconstant, we say that is has differential 
order N if 8 f N \ =£ for some i e I, and 8 f n) = for every j e I and n > N. 

We also set the differential order of a quasiconstant element equal to —00. 
We let Vat be the subalgebra of elements of differential order at most N . 
This gives an increasing sequence of subalgebras 

(4.2) CcJ= V-00 cVocVic-cV, 

such that BVn c Vn+i- Clearly, if V is a field of differential functions, then 
this is a tower of field extensions. 

It is easy to show, using (4.1), that 

(4.3) dV n V N = dV N -i for N > 1 , and dV n V = dT . 
4.2 Normal extensions 

We refine the filtration (4.2) to a filtration V m = Vm.fi c Vm,i c ■ ■ ■ c 
V m j = V m+ i, where 



(4.4) y m>i = / 6 V, 



for all j > i [ c V„ 



Clearly, each subspace V m « is preserved by all partial derivatives s , n) for 



8uf 



(n,j) < (m,i) (in lexicographic order), and it is annihilated by — f-y for 
(n,j) > (m,i). 
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Definition 4.2. The algebra of differential functions V is said to be normal 
if the map — : V m i — ► V m j is surjective for every ie I,rn e Z + . 
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Lemma 4.3. Any algebra of differential function V can be extended to a 
normal one, which can be taken to be a domain provided that V is. 



Proof. Given an algebra of differential functions V and an element a e V m . 

i 

du 



which is not in the image of ^ , one can construct an algebra of differential 



functions V extension of V with an element A e V m i such that (m \ = a. For 
example, we can take the algebra of polynomials in infinitely many variables 

^fc+feo,lH \-k m ,i-lJ^ 



V = V 



dx k du^ ' ... du[ m \ 



m.i—l 



k,k 01 ,. . . ,k m ,i-i £ ^+ 



and define on it a structure of algebra of differential functions by letting J| 
and — fry for (n, j) ^ (m, i — 1) act on — - , V , m '\t A - ; in the obvious 

way (suggested by the notation used to denote the new variables), letting 
—Try for (n, j) > (m,z) act as zero on — ■ —rr— — r — —, and letting 

ouj dx K du\ ' ...au- 1 

^fc+fc ,iH bfc mii _i^ . 0k+k O jA \-k m ,i-\ a 



^ K dx^f 0A . . . du^f— 1 J dx^f ' 1 . . . dutf^- 1 ' 

The lemma follows by standard arguments using Zorn's Lemma. □ 

Example 4.4. The algebra .F[ii^, i e I,n e Z+] of differential polynomials 
over a differential algebra J 7 is a normal algebra of differential functions in 
the variables e I (we can always integrate polynomials). 

Example 4.5. The algebra of differential functions V = JF\u^ , vS n > , n > 1], 
is not normal, since u^ 1 is not in the image of 4^. A normal extension of 
it is V = J-lu- 1 , u( n \ n > 1, log it]. Indeed a preimage via ^ of u m (logit) n , 
m e Z\{— 1}, n e Z + , is obtained by integration by parts 

\duu m (log u) n = —^—u m+l (\ogu) n —\duu m (\ogu) n ~ l , 

m + 1 m + 1 

and a preimage via ^ of u -1 (log u) n is ^j-(log u) n+1 . Similarly, J-[tS n \ n e 
Z +) it^ 1 ,logu^ s ^] is normal for every s. 

Example 4.6. Other examples of normal algebras of differential functions are 
J-[u( n \ n e Z+, e- u ], since we can always integrate by parts P(u)e nu , n e Z, 
for a polynomial P(u). 
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4.3 Variational complex 

For / e V, as usual we denote by §f the image of / in the quotient space 
V/dV. Recall that, by (4.1), we have a well-defined variational derivative 
£ : V/dV -> V®^ given by 

Su- 1 °> a (n)' % b 1 ■ 

Given a set J and an element X e V J , we define the Frechet derivative 
of X as the differential operator Dx{d) : V IL — » V J given by 

p. x ■ 

(4-5) (Dx(d)P) j = 2 Sfg)^, jeJ. 

neZ + iel v u i 

Its adjoint operator is the map D\{d) : V® J —*■ V®^ given by: 



(4.6) {D* X W\= J XH)1§^) ■ 



Here and further, for a possibly infinite set J, V® J denotes the space of 
column vectors in V J with only finitely many non-zero entries. (Though in 
this paper we do not consider any example with infinite £, we still distinguish 
V e and V®^ book-keeping device.) 

The following identity can be checked directly and it will be useful later: 

(4.7) $X ■ D Y (d)P + $Y ■ D x (d)P = $P ■ ^(X ■ Y) , 

for all XeV J Je V® J , P e V e . 

The above notions are linked naturally in the variational complex: 

Tl^-T -> V/dV Av^^Sf-... 

ox 

where is the space of skewadjoint £ x £ matrix differential operators over 
V, and 5{F) = D F (d) - D* F {d), for F e V® e . The construction of the whole 
complex can be found in [DSK09], but we will not need it here. In [BDSK09] 
it is proved that the variational complex is exact, provided that the algebra 
of differential functions V is normal. In particular, if V is normal we have 
that Ker ( -£-) = P+dV, and that F e V®^ is closed, i.e. its Frechet derivative 
Dp{8) is selfadjoint, if and only if it is exact, i.e. F e j^V® 1 . 
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4.4 Differential orders 

Given an arbitrary k x ^-matrix A with entries in V, we define its differential 
order, denoted by dord(^4), as the maximal differential order of all its entries. 

Given a matrix differential operator D = Xir=o ^id 1 e Matfc x ^V[5], we 
define its differential order as 

(4.8) dord(D) = max{dord(^i), . . . , dord(yl n ,)} , 
which should not be confused with its order, defined as 

(4.9) \D\=n if An ^ 0. 

(Note that the notion of order carries over to matrix pseudofferential oper- 
ators, while the differential order is not defined in general.) 

Lemma 4.7. Let D e Matfc x ^ V[o] be a matrix differential operator over V 
and let F eV e . Then: 

(a) dord(DF) sc max{dord(£>), dord(F) + \D\}. 

(b) If D has non- degenerate leading coefficient (meaning that its determi- 
nant is not a zero divisor inV) and it satisfies dord(F) + \D\ > dord(D), 
then doid(DF) = dord(F) + \D\. 

(c) If D has non- degenerate leading coefficient and it satisfies dord(DF) > 
dord(D), then dord{DF) = dord(F) + \D\. 

Proof. Let D = ^™= A s d s . Clearly, for / e V and s e Z+, we have 
dord(/W) = dord(/) + i. Hence, If h > max{dord(£>), dord(F) + \D\}, 
we have 

3=1 s=0 



for every i = l,...,k, proving part (a). Furthermore, if \D\ = n and 
dord(F) + n > dord \D\, we can use (4.1) to get 



3 3 

^ n (dord(F)+n) i DF ) i = Yj 2 ^(dord(F)+n) (A s )ijFj 



j=l s=0 

e dF (n) 1 dF- 

Since, by assumption, the leading coefficient A n e Mat^V of D is non- 
degenerate, the RHS above is non-zero for some i. Hence, dord(DF) = 
dord(i ? ) + n, proving part (b). Part (c) follows from parts (a) and (b). □ 
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4.5 Construction of non-local Poisson structures 

Let V be an algebra of differential functions in u±, . . . , U£. Assume that V is 
a domain, and let K, be the corresponding field of fractions, which is a field 
of differential functions. Let H = (Hij). ^ e Mat^ "^((d" 1 )) be an t x i 
matrix pseudodifferential operator over V, namely 

N 

Hu = ^ H ir,nd n e V((r 1 )) , i,jel. 

n=— oo 

We associate to this matrix H a non-local A-bracket on V given by the 

following Master Formula (cf. [DSK06]) 

(4.10) 

{fxg}H= 2 -%(A + 5) n ^(A + 5)(-A-5r-^r eVftr 1 )). 

i,jel OUj 0U\ 

In particular 

(4.11) {u ix Uj}H = Hji(X) , i,j e I. 

The following result gives a way to check if a matrix pseudodifferential 
operator H e Mat^ x e ^(i^ 1 )) defines a structure of non-local Poisson vertex 
algebra on V. The analogous statement in the local case was proved in 
[BDSK09]. 

Theorem 4.8. Let V be an algebra of differential functions, which is a 
domain, ane let K, be its field of fractions. Let H e Mat^ X £ V((5 -1 )). Then: 

(a) Formula (4.10) gives a well-defined non-local \-bracket on V. 

(b) This non-local X-bracket is skewsymmetric if and only if H is a skew- 
adjoint matrix pseudodifferential operator. 

(c) If H = (Hij).. eI e Mat^ x ^V(5) is a rational matrix pseudodifferen- 
tial operator with coefficients in V, then the corresponding non-local 
X-bracket {-\-}h ■ V x V — > V((A -1 )) (given by equation (4.10)) is 
admissible. 

(d) Let H = {Hij) i j e j e Mat^ x ^V(d) be a skewadjoint rational matrix 
pseudodifferential operator with coefficients in V. Then the non-local 
X-bracket {■ \ defined by (4.10) is a Poisson non-local X-bracket, i.e. 



31 



it satisfies the Jacobi identity (3.3), if and only if the Jacobi identity 
holds on generators (i,j, k e I): 

(4.12) {u ix {u jfi u k } H } H - {u jfi {u ix u k } H }H ~ {{ui\Uj} Hx+ll u k } H = , 
where the equality holds in the space V\ tfl . 

Proof. For the proofs of (a), (b) and (d) one does the same computations 
as in the proof of [BDSK09, Thm.1.15] for the local case. So, we only prove 
part (c). Let a, f,g e V. By the Master Formula (4.10) and the left Leibniz 
rule, we have 

mmh}h = 2 {«A-%y(^ + Brad? + - sr^L )h 

i,jei ou) du\ 
i,jel ou i 

+ S -4(A+/i+s) B w^+^j-M-r4r 

i,jel OU) du\ 



All sums in the above equations are finite. Therefore, all three terms in the 
RHS lie in thanks to Corollary 3.11 and Lemma 2.3. □ 

Definition 4.9. Let V be an algebra of differential functions. A non-local 
Poisson structure on V is a skewadjoint rational matrix pseudodifferential 
operator with coefficients in V, H = \Hij) i j e j £ Mat^ x ^V(5), satisfying 
equation (4.12) for every i,j,k e I. 

Remark 4.10. It is easy to show that, if L e )C(d) is a rational pseudodiffer- 
ential operator, then it can be expanded as 

co N 

(4.13) L =2I! 2 Pl d- 1 o P2 d- 1 o...d- 1 op s d n , 

s=l n=0p!,...,p s eV M 
(finite sum) 

for some fixed M, N e Z+. To see this, write L = AS , where A, S e V[d] 
and S = Yjn=o s nd n nas non-zero leading coefficient sm, and expand S^ 1 



32 



using geometric progression: 



(4.14) < S- 1 = r Ar £( 



-l 



O S 



-1 

N 



On the other hand, it is not hard to see that if L admits an expansion 
as in (4.13), then {a\L(fi)}H £ K\,tt f° r every a e K, and every matrix 
pseudodifferential operator H. As a consequence, if all the entries of a 
matrix pseudodifferential operator H admit an expansion as in (4.13), then 
the corresponding A-bracket {• \ -}h on K, is admissible. 

Remark 4.11. It is claimed in the literature (without a proof) [DN89] that, 
in order to show that a skewadjoint operator H defines a (local) Pois- 
son structure, it suffices to check the Jacobi identity for the Lie bracket 
{• ; ~}h = {• A 'IhI^q in V/dV on triples of elements of the form \fui, where 
/ e T is a quasiconstant. This is indeed true, provided that the algebra of 
quasiconstants T is "big enough", by the following argument. By a straight- 
forward computation, using the Master Formula, we get 



Clearly, this is zero for all f,g,h e J- and all i,j,k e / if and only if H 
is a Poisson structure, provided that the algebra T satisfies the following 
non-degeneracy conditions: 



(ii) if P{d)f = for some differential operator P e V[d] and for all /e J 7 , 
then P = 0. 

Obviously, T fulfills these conditions if it contains the algebra of polynomials 
¥[x]. Often in the literature this criterion is used also for non-local Poisson 
structures, which does not seem to have much sense, since in the non-local 
case V/dV does not have a Lie algebra structure. 

4.6 Examples 

Example 4.12. Let V be any algebra of differential functions in £ differen- 
tial variables, with subalgebra of quasiconstants J- a V. Any skewadjoint 
rational matrix pseudodifferential operator with quasiconstant coefficients, 




(i) if ^ha = for some a e V and all h e J 7 , then a = 0, 
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H = yHij(d)) i - eI e Mat^ x ^ T{d), is a Poisson structure. Indeed, by askewad- 
jointness of H the A-bracket {• \-}h is skewsymmetric, and by the Master 
Formula, all triple A-brackets are zero. Note that, if H e Mat^ -^((i? -1 )) is 
skewadjoint, even if it is not a rational matrix, the corresponding A-bracket 
{■ x ~}h is still admissible, hence it defines a non-local Poisson vertex algebra 
on V. 

In the special case when -ffy(A) = CjjA -1 , and C = (cjj)i,j=l ^ s a sym- 
metric matrix with constant coefficients, we recover the non-local Poisson 
vertex algebras from Example 3.7. When C if a symmetrized Cartan ma- 
trix or extended Cartan matrix of a simple Lie algebra, we get the Poisson 
structure for a Toda lattice (see [Fr98]). 

Example 4.13. The following three operators form a compatible family of 
non-local Poisson structures (i.e. any their linear combination is a non- 
local Poisson structure) on the algebra R± = ¥[u,u' ,u" , . . .] of differential 
polynomials in one variable: 

(i) K\ = d (GFZ Hamiltonina structure), 

(ii) K-\ = d (Toda non-local Poisson structure), 

(iii) H = u'd^ 1 o u' (Sokolov non-local Hamitonian structure), 

First, any linear combination over C of K\ and is a non-local Pois- 
son structure, as discussed in Example 4.12. Next, it is easy to show (cf. 
[BDSK09, Example 3.14]) that if _1 is a symplectic structure over the field of 
fractions /Ci = FracRi, known as the Sokolov symplectic structure, [Sok84]. 
Hence, by Theorem 6.2 below, we deduce that If is a non-local Poisson 
structure. To conclude that K\, K-x, H form a compatible family, it suffices 
to check that 

{u x H(fj,)} K±1 - { Ufl H(X)} K±1 = {H(\) x+ ^u} K±1 , 

where H(\) = u'(d + A) -1 ?/ e V((A -1 )). This is straightforward. 

Example 4.14. Dorfman non-local Poisson structure on the algebra of dif- 
ferential polynomials R\ = ¥[u, u' , u" , . . . ] is: 

h = d- 1 ou'd- 1 ou'd- 1 . 

One easily shows (cf. [BDSK09, Example 3.14]) that H^ 1 is a symplectic 
structure over the field of fractions /Ci = FraciZi, known as Dorfman sym- 
plectic structre, [Dor93], hence H is indeed a non-local Poisson structure. 
Furthermore, one can show, by a lengthy calculation, that Sokolov's and 
Dorfman's non-local Poisson structures are compatible. 
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Example 4.15 (cf. [Dor93]). Another triple of compatible non-local Poisson 
structures on R\ = ¥[u, u',u", . . .] is: 

(i) K\ = d (GFZ Poisson structure), 

(ii) K—i = d^ 1 (Toda non-local Poisson structure), 

(iii) H = d^ 1 ou' + u'd~ 1 (potential Virasoro-Magri non-local Poisson struc- 
ture) . 

Example 4.16 (cf. [Mag80]). There is yet another triple of compatible non- 
local Poisson structures on R± = ¥[u, u', u" , . . . ]: 

(i) K\ = d (GFZ Poisson structure), 

(ii) K 3 = d\ 

(iii) H = d o ud~ l o ud (modified Virasoro-Magri non-local Poisson struc- 
ture) . 

Example 4.17 (cf. [Mag78, Mag80]). The following is a triple of compatible 
non-local Poisson structures on R2 = F[u, v, u' , v', . . . ]: 

(i) K\ = dl (GFZ Poisson structure), 

(iii) H = I V ' C ' _i ^ ^1 U ] (NLS non-local Poisson structure). 

y — uo ovue ouj 

5 Constructing families of compatible non-local 
Poisson structures 

As in the previous sections, let V be an algebra of differential functions in 
the variables ui, . . . , U£, we assume that V is a domain, and we let /C be its 
field of fractions. As in the local case, two non-local Poisson vertex algebra 
A-brackets on V (respectively two non-local Poisson structures) are said to 
be compatible if any their linear combination is again a non-local Poisson 
vertex algebra structure (resp. a non-local Poisson structure). Such a pair 
is called a bi-Poisson structure. More generally, a collection of non-local 
Poisson structures {H a } ae _^ on V, is called compatible if any their (finite) 
linear combination is a non-local Poisson structure over V. 
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Recalling the Jacobi identity (4.12), we introduce the following notation. 
Given rational i x ^-matrix pseudodifferential operators K, H e Mat^ x ^ V(d), 
we let J{H,K) = ^(H^K) - J 2 {H,K) - J 3 {H,K), where J a {H,K) = 
(j°- k (H, K)(\, /x)) . keI , for a = 1,2,3, are the arrays with the following 
entries in Va,^: 

^{H.K^^n) = {u ix {uj^u k } H }K , 

(5.1) J 2 (H,K) ijk (X,fj,) = {u jfi {uixu k } h) k , 
J 3 (H,K) ijk (X,fi) = {{ui X Uj}Hx +lM u k}K ■ 

Consider a collection {H a } aeJ ^ of skewadjoint rational non-local matrix pseu- 
dodifferential operators. By definition, H a is a Poisson structure if and only 
if J{H a , H a ) = 0. It is easy to see that the H ai s form a compatible family 
of Poisson structures if and only if each pair is compatible, i.e. 

(5.2) J(H a ,H l3 ) + J(H^,H a ) = , Va^eA. 

Theorem 5.1. Let H, K e Mat^ X £V(<9) be compatible non-local Poisson 
structures over the the algebra of differential functions V, which is a domain. 
Assume that K is an invertible element of the algebra M&ti x iV(d) . Then 
the following sequence of rational matrix pseudodifferential operators with 
coefficients in V: 

= K , tfM := {HK^Y^H e M&t M V{d) , n>l, 

form a compatible family of non-local Poisson structures overV. 

Remark 5.2. It is stated in [FF81] that H^ n \ n ^ 0, are non-local Poisson 
structures, but the prove there is given only under the additional assupmtion 
that H is invertible as well. In this case the proof becomes much easier since 
is invertible, therefore one needs to prove that (£TW)-i is a symplectic 
structure. 

Following the idea in [TT11], we will reduce the proof of Theorem 5.1 to 
the following special case of it: 

Lemma 5.3. Let H, K e Mat^ X £V(5) be compatible non-local Poisson 
structures over V , and assume that K is an invertible element of the algebra 
Mat^V(d). Then H{d)K- l {d)H{d) e Mat ex£ V{d) is a non-local Poisson 
structure over V . 

Proof. To simplify notation, in this proof we denote H by H, and we let 
R = HK- 1 so that R* = K~ X H. Let = EK~ x h{ = RH = HR*), 
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and let {-A'}2 = {' A '} be the non-local A-bracket on V associated to 
ffP] E MaWV(d) via (4.10). We need to prove the Jacobi identity, i.e. 
using the notation in (5.1), that J(H^,H^) = 0. 

We need to compute all three terms J a = J a (H^ , H^)ijk(X, /x), for 
a = 1,2,3, of the Jacobi identity. First, if / e V and i e I, we have, in 

(5.3) H A /} 2 = £{twa/W&(A), 

(5.4) {«, v /} 2 = 

(5-5) {/A+^fch = ^] i?fe r (A + ^ + 5) {/ A+At «r-}H . 

re/ 

Both the above equations follow immediately from the Master formula (4.10) 
and the definition of H^. The following identities are proved in a similar 
way, using that K o K^ 1 = 1, 

(5.6) {u lX f} H = X>*A + 3/k^ s *(A), 

sel 

(5-7) {«j M /}ff = XK+s/}^'^)' 

tel 

(5.8) {f\+fiUk}H = ^] -Rfcr(A + /i + 5) {f\ +fJ ,U r } K ■ 

rel 

Next, it is not hard to chek, using the left and right Leibniz rules and Lemma 
3.9, that, given an admissible non-local A-bracket {• \ •} on V, the following 
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identities hold in V\n, for every i,j, k e I: 



(5.9) {u ix H^J{n)} = Y u { u i\{ u ty u k]H 

tel 

- 2j -Rfer(A + ^ + 5) {uix{u t yU r }K} 
r,tel 

+ 2 -Rfcr(A + [I + 5) {Uj A {'Uj At 'Ur}H 
reT" 

(5.10) {^-^g ] (A)} = 2{« iM {« M « fc }H 

- X! -Rfcr(A + /U + 5) 
r,sel 

+ ^] -Rfcr(A + 11 + 8) {Uj^{UixU r } H 



MM) 

y=H+d 
y=H+8 



, ^i(A) 

a;=A+<5 



re/ 



(5.11) {Hf{\) x+ ,u k } = Y,{{u sx u J } H 



sel 



A+jU+5 



c=A+<? ' 



x=X+d 



y=H+d 



~ Yj {{ u sxUt}K X +n+d u k} 
s,tel 

tel 

Here and further we use the following notation: given an element 

N 



m,n,p=— oo 



and f,geV, we let 



P(x,y) 

(5.12) 

m,n,p—— oo 



x=A+3 



/ 



i V 

2 P W (A + /i + d) p ((A + 5) m /) (Ox + d) n s) e V A)/1 • 

;P= — 00 

In equation (5.11) we used the assumption that H and X are skewadjoint. 
Combining equations (5.3) and (5.9), equations (5.4) and (5.10), and equa- 
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tions (5.5) and (5.11), we get, respectively, 
(5.13)/ 1 = {iHx{uj pUkjzh 



:=A+5 



J] i U sx{utyU k }H}H 
s,tel 

^ R kr {\ + ^ + d){u sx {utyU r } K }H 
,s,tel 

^ R kr (X + /z + 5) {u sx {u.,yz r }ir}# 



y=n+8 



R$M 



x-\+d 



Km 



y=fi+d 



r,sel 



;=A+<9 



R*si{X) 



(5.14)/ 2 = { Uj Ju ix u k } 2 }2 



Y,{Uj {u sx Uk}H}H[ R*si( X ))( R*M) 

t ~ 1 tV V x=\+d / V y=(j,+d 

V R kr {\ + ^ + d){u t {u sx u r } K } H ( R* sl {\))( R* t j{lA 



r.s.tel 



^] i?fc r (A + /i + 5) {^{^aM^ }# ( y _^ +d R tj{^)) i 



r,te7 



(5.15)/ 3 = {{n iA -u j } 2A+At n fc } 2 



J] i?fe r (A + it + <5) {{n sa ,n i } / / A+/j+5 u r } i /_ 
^] fi fcr (A + /x + d) {{u sx u t } Kx+fi+d u r } H 



, RtiW 

x=X+d 



r,s,tel 



x=X+8 



R 



w (A))i^(/i) 



J] # fcr (A + + d) {{UixU^Hx+^+B^H^Rtjilj) 



■tel 



We need to prove that J 1 - J 2 - J 3 = 0. The first term of the RHS of (5.13) 
combined with the first term of the RHS of (5.14) gives, by the Jacobi 
identity for H and by equation (5.8), 
(5.16) 

Yj [{Usx{utyU k } H }H ~ {Uj y {u sx U k } H }H\ -Rjj(A))( = .g^tjW) 



s,tel 



Y,{{u sx U t }H x+ yU k } H }( R*si( X ))( Rtj(V-) 

y V x=\+d / V v=u+d 



s,tel 



£ R kr (X + ti + d){{u sx ut} Hx+ y u r } K ( i?i(A)Y A*- (m) 



r,s,tel 



=A+3 
y=H+d 
x=\+d 



Similarly, the third term of the RHS of (5.13) combined with the first term 
of the RHS of (5.15) gives, by the Jacobi identity for H and by equation 
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(5.7), 
(5.17) 



^] R kr (X + fJ, + 8) ({u sx {u jfl U r } H }H ~ {{UsxUj}H X+f , + d U r}H^j 



r.sel 



(I x RtiW) = Y, RkriX + n + d)^ Au sx u t } h }h( x R* si (X) 



r,sel 



V # fcr (A + //+<9){«t }ff}jr( ,-R*i(A))( A*,(/«) 



r,s,tel 



In the same way, the third term of the RHS of (5.14) combined with the 
third term of the RHS of (5.15) gives, by the Jacobi identity for H and by 
equation (5.6), 

- Yj R kr( X + H + d) ({u ty {u iX U r } H } H + {{ui\Ut} H X+y U r} H 



r,tel 



(5.18) 



R*M 



x=X+8 



r,tel 

2] Rkr(^ + H + d){u sX {ut y u r } H } K 
r,s,tel 



y=fi+d 



Finally, combining the second term in the RHS of (5.13), (5.14) and (5.15), 
together with the RHS of equations (5.16), (5.17) and (5.18), we get 

J 1 - J 2 - J 3 = ^ R kr (X + fi + d) {u sx {u ty u r } K } H 
r,s,tel 

+ {u ty {u sx U r } K }H + {{u sx U t }K x+v U r }H + {{u sx U t }H x+v Ur}K 



ty l "Si "T J J\ J -« i l l J -fv x-\-y 

{ut y {u sx U r } h) K ~ {u s \{u t yU r } H }K 



x—X+d 



y=H+d 



RtM 



which is zero since, by assumption, H and K are compatible. 



□ 



Remark 5.4. The proof of Lemma 5.3 does not use the assumption that K 
is a Poisson structure. 

Proof of Theorem 5.1. We prove, by induction on n ^ 1, that the rational 
matrix pseudodifferential operators 

ijM = K, fl-M =H,..., fl-M e Mat^V(d) , 

form a compatible family of non-local Poisson structures over V. For n = 1, 
this holds by assumption. Assuming by induction that the statement holds 
for 77, ^ 1, we will prove that it holds for n + 1. Namely, thanks to the 
observations at the beginning of the section, we need to prove that 
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(i) J(#[ n+1 ], = 0, 

(ii) J(iTM , fffa+i] ) + J(#[™+!] , fl-M ) = for every m = 0, . . . , n. 

By the inductive assumption, H = 5j^=o ^i-H"^ is a Poisson structure for 
every xo, . . . ,x n e F. Hence, by Lemma 5.3, we get the following Poisson 
structure for every point (xo, . . . , x n ) e ¥ n+ : 

n In 

HK- X H = ^ x iXj H^K- l H^ = ^ Q P (xo, ■ ■ ■ ,*n)# W , 

i,j=0 p=0 

where, for p = 0, . . . , 2ra, 



(5.19) Q p (x ,...,Xn) = ^ XiXj . 

(i+j=p) 

We thus get 

2n 

= J(HK- l H,HK- l H) = Y i Q 2 P (x ,...,x n )J(H^,H^) 

2n 

+ 2 Q p (x ,...,x n )Q q (x ,...,x n )(j(HW,HW) + J(HM,Hto)), 

p,q=0 

(p<q) 

for every (xo, . . . ,x n ) e F n+1 . Note that, by the inductive assumption, 
J(H&\ fl"M) = for every < p < n and J(ff[ p U M ) + J(# M , # [p] ) = 
for every < p < q < n. Hence the above equation gives 
(5.20) 

2n 

2 Q^ ,...,x n )J(FW,FW) 

p—n+l 

n 2n 

+ S 2 Q p (x ,...,x n )Q q (x ,...,x n ){j(H^,H^) + J(H^,H^)) 

p=0 q=n+l 
2n 

+ 2 Qp(xo,...,x n )Q,(x 0j ... J ^)(J(HW,HM) + J(fr[«l,frW))=0 

p,Q = Tl + l 

(p«?) 

for every (xq, . . . ,x n ) e F n+1 . Next, we introduce a grading in the algebra 
of polynomials in xo,...,x n , letting deg(xi) = i. Then Q p (xq, . . . , x n ) is 
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homogeneous of degree p. By looking at the terms of degree d = 2n + 2 in 

equation (5.20), we get 

(5.21) 

n 

Ql +1 (x ,..., x n )J(H^ , H ) + 2 Q p (x , . . . , x n ) 
Q2n + 2- P {xo, • • • , x n )(J(Hb\ Hl 2n+2 ~ri) + J(H^ 2n+2 -P\H^)) = , 

while, by looking at the terms of degree d = m + n + 1 with me {0, . . . , n} 
in equation (5.20), we get 

773 

, . y, Qp{ x 0: • • • i x n)Qm+n+l—p( x 0i ■ ■ ■ j x n) 

( 5 - 22 ) p% 

(j(tfM,if[ m+n+1 - p ]) + j(H^ m+n+1 - p \H^)) =0, 

for every (xo, . . . , £ n ) e F n+1 . To conclude the proof, we only need to show 
that equations (5.21) and (5.22) imply respectively relations (i) and (ii) 
above. This is a consequence of the following lemma. 

Lemma 5.5. (a) For every n ^ 1, 
(5.23) 

Q 2 n+l {xo, ...,x n )$ Span F j Q p (x Q , x n )Q 2 n+2-p(xo, ...,x n )\ 

< J 2^p^n 



(b) For every n ^ 1 and me {0, . . . ,n}, 

(5.24) QmQn+i i Span F ] QpQm+n+l-p \ 

Proof. Note that, 



Q p (0, ...,0,x k ,.. 


• ; 2-n) 






i ° 2 


if p < 2k 


= ^ X{Xj = j 




\ip = 2k 


i,j=k 


[ 2xfcXp_fc + . 


if p > 2k 


(i+j=p) 







We prove part (a) separately in the cases when n is even and odd. If 
n = 2k — 1 is odd, letting xq = ■ ■ ■ = x k _\ = we have Q n+ i = x\ 0, 
and Q p = for all p = 2, . . . , n = 2k — 1. This implies (5.23) for odd n. If 
n = 2, we have Q2 = 2x0^2 + xf , Q3 = 2xix 2 , Q4 = x 2 , hence (J§ ^ FQ2<?4- 
If n = 2k with k ^ 2, letting xo = • • • = £fc-i = we have Q p = for all 
p = 2, . . . , n - 1, Q n = x|, Q n+ i = 2x fe a; fe+ i, Q n+2 = 2x fc x fe+2 + Since 
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Qn+i = ^ x k x k+i is not a mu ltipl e °f QnQn+2 = 2x\xk+2 + x 2 k x 2 k+1 , (5.23) 
holds for even n. 

Similarly, we prove part (b) separately in the cases when m is even and 
odd. If m = 2k is even, letting xq = ■ ■ ■ = Xk-\ = we have Q m Q n +i = 
x kQn+i 0, and Q p = for all p = 2, . . . , m — 1. Hence (5.24) holds for even 
m. For m = 1 < n, we have Qo = Xq, Q\ = 1xqX\, Q n +i = 2x\x n + . . . . 
Therefore, QoQn+2 is divisible by Xq, while Q\Q n +\ = 2x^x\{2xix n + . . . ) 
is not. Finally, if m = 2k + 1 is odd, with k ^ 1, letting xq = ■ ■ ■ = Xk~i = 
we have Q p = for all p = 2, ... ,m — 2, Q m -i = x l, Qm = 2x k x k+ i, 

Qn+i = 2x k x n+ i- k + 2xk+ix n _ k + Hence, Q m -iQn+2 = x 2 k Q n+2 is 

divisible by x\, while Q m Q n +i = ^x 2 k x k 

not, proving (5.24) for odd m. □ 

□ 

Example 5.6. Let K = d 3 , H = d 2 o o \d 2 . These are compatible Hamil- 
tonian structures (see [DSKW10]). Hence, by Theorem 5.1, 

ffH = (HK-^^H = d 2 o(-o 3) 2n o B , n e Z+ , 

are compatible Poisson structures. This was proved in [DSK W10] by direct 
verification, and deduced from Theorem 5.1 in [TT11]. 



6 Symplectic structures and Dirac structures in 
terms of non-local Poisson structures 

6.1 Simplectic structure as inverse of a non-local Poisson 
structure 

As in the previous sections, let V be an algebra of differential functions in 
the variables u\, . . . , ug, which is a domain, and let K, be its field of fractions. 

Recall that (see e.g. [BDSK09]) a (local) symplectic structure on V is an 
t x I matrix differential operator S = \Sij(d))ij e j E Mat^ x ^V[d] which is 
skewadjoint and satisfies the following symplectic identity: 

(6-1) £ (^#A« - + (-A - „ - 5r^#) = . 

We can write the symplectic identity (6.1) in terms of the Beltrami A- 
bracket <■*■>: V x V ^ V[A], introduced in [BDSK09]. It is defined as the 
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symmetric A-bracket such that (ui X Uj) = Sij, and extended by the Master 
Formula (4.10): 



Then, the symplectic identity (6.1) becomes 

(6.2) (uj x {u ilx u k }s) - (u ifl {u jx u k }s) + ({ u j\ u t}s x+fj: u k ) = , 

where, recalling (4.10), we let {uj x Ui}s = Sij(X). 

Note that, if S e Mat^ x ^V(5) is a rational matrix pseudo differential 
operator with coefficients in V, then, by Corollary 3.11, all three terms in 
the LHS of equation (6.2) lie in Va,^. Hence, equation (6.2) still makes sense 
(as an equation in Va,^). 

Definition 6.1. A non-local symplectic structure on V is a skewadjoint 
rational matrix pseudodifferential operator S = (SV,(d)) . . gJ e M&t£ X £V(d) 
with coefficients in V, satisftying equation (6.2) in V\ n for all k e I. 

Theorem 6.2. Let S e Mat^ X £ V(d) be a skewadjoint rational matrix pseu- 
dodifferential operator with coefficients in the algebra of differential functions 
V. Assume that S is an invertible element of the algebra Mat^ X £ V(d). Then, 
S is a non-local symplectic structure over V if and only if S^ 1 is a non-local 
Poisson structure over V . 

Proof. Clearly, S is skewadjoint if and only if is skewadjoint. Hence, 
recalling the Definition 4.9 of non-local Poisson structure, we only need to 
show that equation (6.1) in Va,^ is equivalent to the Jacobi identity (4.12), 
again in V x>fl , for H = S^ 1 . By equation (3.6), Remark 3.12, and the Master 
Formula (4.10), we have, letting Sij(d) = X! P =-oo s ij;pd P i 




(6.3) 



{uix{uj /J u k } H }H = {ui\{S ^kji^s-l = 



I N 



-2 2 {s-^ux + ^ + d^srtjs-^ + dYis- 1 )^) 



r,t=l p=— co 



= - 2 {S- l ) kr {\ + M + d) ((A + d) n {S- l ) si {\)) 



r,s,tel, neZ-j. 
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Exchanging % with j and A with fi, we get 

{uj^{u ix u k } H } H = - Yj (S~ 1 ) kr (\ + fj, + d) 

, . r,s,teI,neZ + 

fa + gr(5-%^))( g^^ (5-^(A)). 

Similarly, by equation (3.7) and Remark 3.12, we have, using the assumption 
that S 1 is skewadjoint, 



{{ui\Uj}H x+h u k }H = {(S 1 )ji(A) A+At n fc } 5 -i 

l N 

= 22 {^; PA+/1+c5 n fc }5-^(5- 1 )y(^)(A + ^(5- 1 ) sl (A) 



s,t=lp=— 00 

(6 - 5) = 2 (5-V(A + M + 5)(-A- M -ay 



r,s,tel, msZ + 



((g- 1 )^(M) ^; ( t; g) (5- i ).(A)) 



Combining equations (6.3), (6.4) and (6.5), we get that the LHS of the 
Jacobi identity (4.12) is 

, RR * = 2j ( 5 ^r(X + ^ + d)\-—^ r x+—^ r y 

I 6 - 6 ) r,s,teJ,neZ + \ du\ 

5r ^W) (L =A+5 ( S -)..(A)) (l^/s-'leW) , 

where we used the notation introduced in (5.12). Clearly, the RHS of (6.6) 
is zero, provided that the symplectic identity (6.1) holds. For the opposite 
implication, we have, by (6.6), 

Yj S lk {x + y + 8) I {u ix {u jy u k } H }H ~ {uj y {u ix u k } H }H 

i,j,kel \ 



-{{u ix Uj}H x+y U k }H | ( x=x+d S ia( X )) ( 



y=H+d 



" Zj " - (n) A + Z (n) ^ - C— A - A* - C; („) • 

nez+ \ ^4 <9u^ Buy J 

Hence, equation (4.12) implies equation (6.1). □ 
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6.2 Dirac structure in terms of non-local Poisson structure 

Let V be an algebra of differential functions, which is a domain, and let K, 
be its field of fractions. 

We have the usual pairing V® e x V 1 — » V/dV given by (F\P) = §F ■ P. 
This pairing is non-degenerate (see e.g. [BDSK09, Prop. 1.3(a)]). We extend 
it to a non-degenerate symmetric bilinear form 

(6.7) <•!•>: (V® £ © V e ) x (V® e © V e ) - V/dV , 

given by (F © P\G © Q) = \{F ■ Q + G ■ P). 

The Courant-Dorfman product o is the following product on the space 

(6.8) (F©P)o(GffiQ) = (D G (d)P+D* P (d)G-D F (d)Q+D* F (d)Q)®[P,Q], 
where, for P, Q e V , we let 

(6.9) [P,Q]=D Q (d)P-D P (d)Q. 

By definition, if F e V® £ is closed (cf. Section 4.3), we have Dp(d)Q — 
D F (d)Q = 0. Moreover, it is straightforward to check that, for arbitrary 
G e V® e and P e V e , we have 

D* G {d)P + D P {d)G=^\P-G. 

Hence, formula (6.8) takes a simpler form when F and G are closed elements 
of V® e : 

(6.10) (F®P)o(G(BQ) = l-^p.G^e[P,Q]. 

Remark 6.3. All the above notions have a natural interpretation from the 
point of view of variational calculus. Indeed, the space V is naturally iden- 
tified with the Lie algebra of evolutionary vector fields g s , and the space 
V® is naturally identified with the space of variational 1-forms J) 1 . Then 
the contraction of variational 1-forms by evolutionary vector fields gives the 
inner product (6.7); the Courant-Dorfman product corresponds to the de- 
rived bracket [• , -]d, where [• , •] is the Lie superalgebra bracket on the space 
of endomorphisms of the space of all de Rham forms over V, and d = ad(#), 
where 5 is the de Rham differential, [BDSK09, Prop. 4. 2]. 
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Definition 6.4 ([Dor93, BDSK09]). A Dirac structure is a subspace C c 
V®^©V^, which is maximal isotropic with respect to the inner product (6.7), 
and which is closed under the Courant-Dorfman product (6.8). 

Remark 6.5. If £ <= V® e © V e is a Dirac structure, then the subspace 

g = [F © P e C\ F is closed } c £ 

is a Lie algebra with respect to the Courant-Dorfman product, and its de- 
rived Lie algebra lies in the subalgebra 



Indeed, the Courant-Dorfman product (6.8) on V® £ © V satisfies the left 
Jacobi identity, cf. [BDSK09, Sec. 4. 2]. Moreover, by the isotropicity as- 
sumption on L we have \P G = — \Q ■ F for F © P, G©Qe£, so that the 
Courant-Dorfman product restricted to g, given by formula (6.10) is also 
skewsymmetric. Therefore q is a Lie algebra with respect to o, and, again 
by formula (6.10), f) contains the derived subalgebra gog. 

Given two £ x t matrix differential operators A,B e Mat^ X £ V[d] consider 
the following subspace of V® e © V e : 



If, moreover, B is non- degenerate, then (6.12) holds if and only if AB^ 1 e 
Mat^ X £ /C((5 -1 )) is skewadjiont, while if A is non- degenerate then (6.12) 
holds if and only if BA^ 1 e Mat£ X £ /C((5 -1 )) is skewadjiont. 

Proof. For X,Y e V® e we have 



{B{d)X®A{d)X\B{d)Y®A{d)Y) = \Y ■ (A*(d)B(d) + B* (d)A(d)) X . 



Hence, due to non-degeeracy of the pairing (F\P) = §F-P, the space Ca,b is 
isotropic if and only if (6.12) holds. The remaining statements are straight- 





(6.12) 



A* 



B + B*A = 0. 



forward. 



□ 
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Example 6.7. Letting A e Mate x e V and B = ILj 5, condition (6.12) holds if 
and only if A is a symmetric matrix with entries in C a V (the subring of 
constant functions). In this case AB" 1 is a skewadjont matrix pseudodiffer- 
ential operator and C a ,b = {AX © dX \ X e V® e } is an isotropic subspace 
of V®^ © V . It is not hard to show directly that Ca,b is maximal isotropic 
if and only if the matrix A is non-degenerate. When V = fC is a differential 
field, this is a corollary of the following general result: 

Proposition 6.8 ([CDSK12b]). Let K. be a differential field, and let H = 
AB -1 be a minimal fractional decomposition of the skewadjoint rational 
matrix pseudodifferential operator H e Mati x glC(d). Then the subspace 
C A b c © K e is maximal isotropic with respect to the inner product 
(6.7). 

Proposition 6.9. Suppose that A, B e Mat^ X £ V[<9] satisfy equation (6.12). 
Then the following conditions are equivalent: 

(i) (X o Y, Z) = for all X,Y,Ze C a ,b- 

(ii) for every F,G e one has: 

A*(d)D mG (d)A(d)F + A*(d)D* A(d)F (d)B(d)G 
(6.13) -A*{8)D B{8)F {8)A{d)G + A*(d)D* B{8)F (d)A(d)G 

+B*(8)D A{d)G (8)A(8)F - B*(8)D A(d)F (8)A(8)G = 0. 

(Hi) for every i,j,k e I, one has in the space V[X,ff\: 
(6.14) 

2 U s (\ + v + d) f^£(A + 8) n Ati(X) - + 9TAf (M)) 

+sjE.(A + v + d) ( d -^-(x + 8TMX) - + drA tj („)) 

+At s {X + ^ + 8){-\-^-dT^^BM + d -^A tj {^ = 0. 

\ OU s OU s J J 

Proof. Letting X = B{8)F © A(8)F, Y = B{3)G © A(8)G, Z = B{d)E © 
A(d)E, condition (i) reads 

${A(d)E) ■ (D B{d)G {d)A{d)F + D* A{d)F (8)B{8)G - D B(8)F (8)A(8)G 
+^* (a)F (5)A(5)G) + (S(5)£;)-(£> A(5)G (5)A(5)F- J D A(5)F (5)A(5)G)=0. 
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Since the above equation holds for every E e V®^ it reduces, integrating by 
parts, to equation (6.13). For this we use the non-degeneracy of the pairing 
(6.7). This proves that conditions (i) and (ii) are equivalent. 

We next prove that conditions (ii) and (iii) are equivalent, provided 
that (6.12) holds. For a = 1, ... ,6, let (6.13) a be the fc-entry of the a-th 
term of the LHS of (6.13): for example (6.13)i = (A*(d)D B{8)G (d)A(d)F) k . 
We have, by the definition of the Frechet derivative and some algebraic 
manipulations (similar to those used in the proof of [BDSK09, Prop. 1.16]), 

(6.13)! = 22 (A^f^^-G^Mm 

i,j,s,tel neZ+ ^ du t 

+AUd)B sj {d)^-d n A ti {3)F i 



(6-13) 2 = 22 [At s m-a) n {?^Fi)BtMG. 

i,j,s,tel neZ+ ^ du s 
dF- 

+ A* s (d)(S) n ^A* t (d)B tj (d)G j 
(6-13)3 = 2 2 {AUd^^^F^A^G 



3 



i,j,s,tel neZ + y 

du t 
dF- 

+ AUd)Bsi(d)— ( ^d n A tj (d)G j 



(6-i3) 4 =22 (^sm-dr{ d -^F t )A tJ (d) Gj 



i,j,s,tel neZ+ v ' du s 

+AtM(-d) n ^B*(d)A tj (d)G j 
dul 

(6-13) 5 =22 (BK^^^-G^Mm 

' GUI ' 



i,j,s,teIneZ+ x ua t 

dG 



+B* ks (d)A SJ (d)^d n A ti (d)F i 
du: 



(6-13) 6 = 2 2 (Bt(8)(^§-F^A tj (d) Gj 

i,j,s,tel neZ+ ^ du t 

dF- N 
+Bt s {d)A m {d)-^ ) d n A t] {d)G J 



du\ 



Combining the second terms in the RHS of (6.13)i and (6.13)s we get zero, 
thanks to equation (6.12). Similarly, we get zero if we combine the second 
terms in the RHS of (6.13)2 and (6.13)4, and if we combine the second 
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terms in the RHS of (6.13)3 and (6.13)6- Equation (6.14) thus follows from 
equation (6.13) once we replace d acting on Fi by A, and d acting on Gj by 
fj,. □ 

Remark 6.10. It follows by Definition 6.4 and Proposition 6.9 that a Dirac 
structure is a maximal isotropic subspace C of V®^ © V e satisfying one of 
the equivalent conditions (i)-(iii) above. 

Proposition 6.11. Suppose that A,B e M&ti X £ V[d] satisfy equation (6.12) 
and assume that B is non- degenerate. Suppose, moreover, that the (skewad- 
joint) rational matrix pseudodifferential operator H = AB^ 1 has coefficients 
in V, i.e. He Mat^ V{d). Consider the corresponding non-local \-bracket 
{• A ~}h given by the Master Formula (4.10). Then the Jacobi identity (4.12) 
on {'A'}ff is equivalent to equation (6.14) on the entries of matrices A and 
B. 

Proof. Letting A st (d) = Y^ =0 a st , m d m and B st {d) = Sf =0 6 lj;n 5 n . By for- 
mula (4.11) and the left Leibniz rule (3.4) we have, 

{Ui\{u jfl U k } H }H = Yj{ u i\ A kr{V + d)B~ 1 (fl)} H 
rel 

M 

= S 2 { u i\ a kr;m}H(n + S) m B~\fJ,) 
rel m=0 

+ Y i A kr {\ + [i + d){u iX B~^(fi)} H . 

rel 

By equation (3.8) we have 
{^xB^(fi)} H 

M 

= ~ S 2 ( i? " 1 )-( A + V + S){Ui X bst;m}H^ + d) m (B- r ) tj (,l) . 
s,tel m=0 

Combining the above two equations we then get, using the Master Formula 
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(4.10), 

{uix{u jlx U k } H }H 



M 



da 



S 2 E ^(^ + ^V))((A + 5)M st (A + ^-HA)) 

Af 

2 S 2 ^r(A + ^ + 5)(S- 1 ) r . s (A + /i + ( ?) 
p,q,r,s,tel m—O neZ+ 
SbgUm 



cu 



2 S 

r,s,tel nsZ+ 



- 2 2 ^r(A + /i + 5)( J B- 1 ) rs (A + / u + 5) 

p,q,r,s,tel neZ + 

OUp 

Next, we apply B^, k (X + \i + d) to both sides of the above equation (on the 
left), replace A by A + d acting on Bn>(X), replace fj, by fi + d acting on 
Bjj/((j,), and sum over i,j,k e I. As a result we get, using the assumption 
(6.12) (see notation (5.12)), 
(6.15) 

J] B^ k (X + n + d){u ix {uj u k } H } H ( J3«>(A))( Bjj'ip) 



i,j,kel 



+< 

2 2 (s fc %(A + M+5)^^(A + 5r^(A) 
+^, fc (A + /, + 5)^^(A 



(n) 



Exchanging i' with j' and A with fi in (6.15), we get 
(6.16) 

T, B k'k( x + H + 8){uj {u lx u k } H }H{ Bui(X))( B jf (n)) 

. T~r _ a \ x—X+o / V y=fj,+d / 



i,j,kel 



■ Y, B k *>k( X + V + d ){ u ix{ujU k } H }H[ Bij'(ji) 
i,j,kel 



k,jel neZ^ 



2 2 Kva+m+3) 



X = ^L + 

dA ki ,{\) 



cu 



y=X+8 

( f i + d) n A jjl ( f i) 



\-AU(\ + m + 5)^rr(^ + 5) n %'(/^) 
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We are left to study the third term in the Jacobi identity. By the right 
Leibniz rule (3.4) we get, 



{{uiXU^Hx^U^H = 2{A?>( A + d)B ri 1 (X)x+ (i Uk}H 
rel 

M 

= 2 2 {a 3 r, mx+ , +d u k } H J\ + d) m B-^\) 

rel m=0 
M 

+ 2 2 {B~ l W\+fi+dUk\H^{-H - d) m a jr ^ m . 



By equation (3.9) we have 

I M 

= ~ 2 2^m A+/1+a ^}^o((A + 5r(fi-%(A))(B*- 1 ) sr ( M + 5). 



,i=l m =0 



Combining the above two equations and using the Master Formula (4.10) 
we then get 

{{UixUj} Hx+fi U k } H 

2 A kt (X + » + d)B u \X + M + d){—\ — M — g)" g^+^ B-i(A) 



r,s,tel 



- 2 A fc ,(A + A t + 5)B- 1 (A + M + 5)(-A-M-5) n 

r,p,q,s,tel 
neZ + 

OUp 

Hence, if we apply, as before, B?, k (X + fi + d) on the left, replace A by 
A + d acting on B U :(X), replace fi by fj, + d acting on Bjji(pL), and sum over 
i,j,k e /, we get, using (6.12), 



(6.17) 



2 B k'k( X +V+ d ){{ u ixUj}H x+ yU k } H y ^ Bii>{X) i 
i,j,kel 

B n>(»)) =-22 A ^ X + m + 5)(-a - m - ey 

j,kel neZ + 



y=H+d 
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Combining (6.15), (6.16) and (6.17), we get that the expression 



B k'ki x + M + d) [{ u ix{uj y u k } H }H ~ {u jy {u ix u k } H } H 

j,kel 

-{{UixUjjH^yU^H) \\ x=x+a B iv( X )) (| = +d B 3i'^)) 



is the same as the LHS of (6.14). The claim follows from the assumption 
that the matrix B e Mat£ X £ V[d] has non-zero Dieudonne determinant. □ 

Theorem 6.12. Let V be an algebra of differential functions in the variables 
u±, . . . ,U£, which is a domain, and let K, be its field of fractions. Let H = 
AB- 1 , with A, Be M.&t M V[d], B non-degenerate, be a minimal fractional 
decomposition (cf. Definition 2.12 and Remark 2.14) of the rational matrix 
pseudodifferential operator H e Mat^ V(d). Then the subspace 

(6.18) C A>B {K) = {B(d)X © A(S)A \ X e /C® £ } c K,® 1 © K l , 

is a Dirac structure if and only if H is a non-local Poisson structure over 
V. 

Proof. It immediately follows from Remark 2.14 and Propositions 6.6, 6.8, 
6.9 and 6.11. □ 

Remark 6.13. We may define a "generalized" Dirac structure as a subspace 
C of V® 1 © V e , such that £ c (i.e. C is isotropic), and £ o £ <z C 1 
(i.e. condition (i) in Proposition 6.9 holds), where C~*~ is the orthogonal 
complement to C with respect to the inner product (6.7). Note that a 
Dirac structure is a special case of this when C is maximal isotropic. If 
A, B e Mat^V[<9], with B non-degenerate, then Ca,b is a generalized 
Dirac structure if and only if H = AB" 1 is a non-local Poisson structure 
over V (not necessarily in its minimal fractional decomposition). Note also 
that any subspace of a generalized Dirac structure is a generalized Dirac 
structure. 



6.3 Compatible pairs of Dirac structures 

The notion of compatibility of Dirac structures was introduced by Gelfand 
and Dorfman [GD80], [Dor93] (see also [BDSK09]). In this paper we intro- 
duce a weaker, but more natural, notion of compatibility, which still can 
be used to implement successfully the Lenard-Magri scheme of integrabil- 
ity, and which is more closely related to the notion of compatibility of the 
corresponding non-local Poisson structures. 
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Given two Dirac structures C and £ a V®^© we define the relations 
.19) 



Mc,c> = {P®P' 
Mc,c> = {F®F' 



F®PeC,F@P' eC for some F e V® e } cV'eV*, 
FQPe C, F'@Pe C for some PeV'Jc v®^©V® £ . 



Definition 6.14. Two Dirac structures C, £' a V®^ © V are said to be 
compatible if for all P, P' ,Q,Q' e V e , F, F' , F" e V® £ such that 

PffiP', QffiG/ e M c ,a and F © F', F' © F" e Mc,c , 

we have 

(6.20) (F\[P,Q]) - (F'|[P,Q']) - (F'|[P',Q]) + (F"\[P',Q']) = 0, 
where, as before, (F|P) = $F-P, and, for P,Q e V e , [P, Q] is given by (6.9). 

Remark 6.15. The original notion of compatibility, introduced by Dorfman 
in [Dor93], is similar, except that Mc,a is replaced by the "dual" relation 

Af* £, = {F © F' e V® £ © V® £ | \F ■ P = JF' • P' for all P © P' e A/" A /y} . 

Since £ and £' are isotropic, we have, for F © F' e Mc,c'i an d for Q®Q'e 
M"c,c; $ F ■ Q = -JG • P = JF' • Q', where PeV'andGe V®^ are such that 
F © P, G © Q e C, F' © P, G © Q' e C. Hence, A/i> c ^c,a ■ 

Even with the weaker notion of compatibility, the following important 
theorem still holds (cf. [BDSK09, Thm.4.13]). 

Theorem 6.16. Let (C,C) be a pair of compatible Dirac structures. Let 
F ,F 1 ,F 2 e V® e be such that: 

(i) Fq and F\ are closed, i.e. Dp (d) = Dp n (d) for n = 0, 1; 

(ii) F © Fi , Fi © F 2 e Afc,c ■ 

Then, for all P®P',Q®Q'e Mc,c, we have 

(6.21) lQ'.(D F2 (d)-D F2 (d))P' = 0. 
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Proof. By the assumption (6.20), we have 
= (F \[P,Q\) - (Fi\[P,Q']) - (Ft\[P',Q]) + (F 2 \[P',Q']) 
= J (F ■ D Q {d)P - F ■ D P (d)Q - F l ■ D Q ,{d)P + Ft ■ D P {d)Q' 

-Ft ■ D Q (d)P' + Ft ■ D P ,{d)Q + F 2 ■ D Q ,{d)P' - F 2 ■ D P ,(d)Q'] 

=\p- - - \ Q ■ t^ p) - 

- J p' • ^WQ) - (W)) + JV ■ - (^1^0) 

- J Q ■ D Fo (d)P + J P ■ D Fo (d)Q + J Q' ■ D Fl (3)P - J P ■ D Fl {d)Q' 

+ ^Q-D Fl (8)P' ~jp'- D Fl (3)Q -jtf. D F2 (8)P > + J P' ■ Dp 2 {d)Q' . 

In the second identity we used the definition (6.9) of the Lie bracket on 
V , and in the last identity we used equation (4.7). Since, by assumption, 
F ®Ft e Mc,c and Q®Q' e N c ,c>, we have (by Remark 6.15) that (F \Q) = 
(Ft\Q'). Hence the first term in the RHS above is zero, and, by the same 
argument, the first four terms are zero. The following six terms are also 
zero since, by assumption, Dp (d) and Dp 1 (d) are selfadjoint. In conclusion, 
equation (6.21) holds. □ 

6.4 Compatible non-local Poisson structures and the corre- 
sponding compatible pairs of Dirac structures 

In Theorem 6.12 we proved that to a non-local Hamiltonian structure H e 
Mat^ x ^ V(d) in its minimal fractional decomposition H = AB -1 , with A, B e 
Mat^Vfd], there corresponds a Dirac structure Ca,b(K) over the field of 
frations /C. In this section we prove that to a compatible pair of non-local 
Poisson structures H = AB^ 1 , K = CD^ 1 , in their minimal fractional 
decompositions, there corresponds a compatible pair of Dirac structures 
£a,b(^C), £c,D(fc) over /C. This is stated in the following: 

Theorem 6.17. Let V be an algebra of differential functions in ut, ■ ■ ■ ,ug, 
which is a domain, and let fC be its field of fractions. Let H, K e Mat^ V(<9) 
be compatible non-local Poisson structures over V. Let H = AB^ 1 , K = 
CD^ 1 be their minimal fractional decompositions (cf. Definition 2.12 and 
Remark 2.14)- Then /^^(/C) and Cc.d(K>) are compatible Dirac structures 
over /C. 

By Theorem 4.8, the Poisson structures H and K over V are compatible 
if and only if we have the following "mixed" Jacobi identity on generators 
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(i,j,ke I): 
(6.22 



{ui\{uj u u k } H } k - {u jix {u iX u k } H } K - {{u iX Uj} Hx+fj u k } K 



^-{Uixiuj^Uhj^H - {u jtl {Ui X U k } K }H - {{Ui\Uj}K x+f ,U k } H = 0, 

In order to relate the above condition to the compatibility of the cor- 
responding Dirac structures Ca,b and £>c,d, we need to compute explicitly 
each term of the above equation. This is done in the following: 

Lemma 6.18. Suppose that the pairs (A,B) and (C,D), with A,B,C,D e 
Mat£ X £ V[d], satisfy equation (6.12): 

(6.23) A*B + B*A = , C*D + D*C = 0. 

Assume that B and D are non- degenerate, and that the (skewadjoint) ra- 
tional matrix pseudodijjerential operators H = AB^ 1 and K = CD" 1 have 
coefficients in V, i.e. H,K e Mat^ V(d). Consider the corresponding non- 
local X-brackets {-\-}h and {-\-}k given by the Master Formula (4.10). 
Then, in terms of notation (5.12), we have the following identities for every 
i',f,k'el: 
(6.24) 

Yi B k'k( X + M + d ){ u ix{Uj y U k } H }K{\ x=x+d D i¥ (\)) (\ y= ^ +g B jf (fj)) 
i,j,kel 



2 S Bt^x+^+df-^iix+drc^x) 



2 S At, k (X + » + d)^>^(X + drC ii ,(X) 



a ( n ) 

i,kel neZ + VU { 



(6.25) 

Yi D k'k( x + V+ d ){ u ix{uj y u k } K }H{\ x ^ x+d B iV {X)) {\ y=fl+d D jf {fx)) 

i,j,kel 

2 2 DUX + » + d) 8 -^^(X + drA u ,(X) 



i,kel neZ+ du. 



■S 2 c^x+v+df-^^ix+drAi'W, 

i,kel neZ+ 



5G 



(6.26) 

2 B k'k( X + V + d ){ u iy{ u ixUk}H}K{\ x=x+d Bn:(X)) {\ y=ix+g D jr {n)) 
i,j,kel 

= 2 2 B^ix+^+df-^^-^+drc^) 

j,keIneZ + OU- 

+ 22 Aux+v+df-^^^+drc^), 

j,kel neZ + V u j 

(6.27) 

2 D k'k( X +V+ d ){ U jy{ U ix U k}K}H(\ x=X+d Dii,(X)) (\ y=fl+g Bjf(n)) 

i,j,kel 

= 2 2 D^ + V + Q^^fr + WAM 



j,keIneZ+ uu j 

+ 2 2 cz fk (\ + v+d) d -^^(v+dTA^), 

j,keIneZ + ^U- 



(6.28) 

i,j,kel 

= -22 C ^( A + m + 5)(-a - m - ar^^BtfM 

- 2 2 ^ A + M + 5)("A - M - W^&AjfM , 

j,keIneZ + Su k 

(6.29) 

2 B k'k( x + M + ^{{"^^lic^+^fclifd^A+eAi'CA)) (I^^^w'Cm)) 

= -22 ^( A + a* + 5)(-a - M - ffr^^-D^) 

j,kel neZ+ 0U k 

- 2 2 ^ A + /* + 5)(-a - - sr d -^^-c ir {p) . 

j,keIneZ+ 0U k 

Proof. For equation (6.24), we can use the Leibniz rule and equation (3.8) 
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to get 

{ui x {uj y u k }H}K = J^{u ix A kr (y + d)(B~ 1 ) rj (y)} K 

rel 

= Yl { U ixakr;m}K{y + d) m {B- 1 ) rj (y) 

rel meZ + 

^ A kr (x + y + d)(B- 1 ) rp (x + y + d){u ix b m - m } K (y + d) m {B' 1 ) qj (y) . 

r,p,qel 

We can then use the Master Formula (4.10) to get 

{u tx {u 3y u k } H } K =Y, 2 ( dAkr[V ^ ^ (B-%(y)) (x + d) n K sl {x) 

r,seIneZ + Su s 

^ A kr (x + y + d){B- x ) rp {x + y + 3) ( ^^^ (iT 1 )^)) 

neZ+ 

x(x + d) n K si {x). 

If we now replace x with A + d acting on Dai (A) and y by fi + d acting on 
Bjji(ii), and we apply By k {\ + /i + d), acting from the left, to both sides of 
the above equation, we get, after using the assuption (6.23), that equation 
(6.24) holds. Equation (6.25) is obtained from (6.24) by exchanging the roles 
of H and K. Equation (6.26) is obtained from (6.24) by exchanging A with 
/i and i and i' with j and f respectively, and equation (6.27) is obtained 
from (6.26) by exchaing the roles of H and K. Finally, equations (6.28) and 
(6.29) can be derived with a similar computation, which involves the right 
Leibniz rule (instead of the left) and equation (3.9) (instead of (3.8)). □ 

Let us next describe the relations (6.19) associated to Dirac structures 
C-A,B{fc) and Cc,d{1£) defined in (6.18). We have 
(6.30) 

Mc AB {K),c c nW = {A(d)X © C(d)X> \X,X'e K® 1 , B{d)X = D{d)X'} , 
Mc A , B (ncc.nVQ = {B{8)Z®D{S)Z'\Z,Z'eK®^ A{d)Z = C{d)Z'} . 

Hence, by Definition 6.14, the Dirac structures Ca,b and Cc,d are compat- 
ible if and only if, for every X, X', Y, Y', Z, Z', W, W' e V® 1 such that 

B(d)X = D(d)X' , B(3)Y = D(d)Y' , B(d)W = D(8)Z' , 

(6.31) 

A{d)Z = C{d)Z' , A(S)W = C(d)W , 
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we have the following identity: 

[B(d)Z\[A(d)X,A(d)Y]) - (D(d)Z'\[A(d)X,C(d)Y']) 

6.32) 

- {B(d)W\ [C(d)X', A(d)Y]) + (D(d)W'\ [C(d)X', C(d)Y']) = . 

Lemma 6.19. Suppose that H = AB~ l and K = CD^ 1 are non-local 
Poisson structures, and that conditions (6.31) hold. Then equation (6.32) 
is equivalent to the following equation: 
(6.33) 

- l(A(d)Y) ■ D mx (d)C(d)Z> + l(A(d)Y) ■ D% {d)x (d)C(d)Z' 
+ l(B(d)Y) ■ D c(8)z ,(d)A(d)X - \(B(d)Y) ■ D mx {d)C{d)Z' 
+ l(C{d)Y') ■ D D{d)z ,(d)A(d)X + l(C(d)Y') ■ D* A{8)x (d)D(d)Z' 
+ S(A(d)Y) • D B{8)w {d)C{d)X' + \(A(d)Y) • D* (8)x ,(d)B(d)W 

- l(C(d)Y') ■ D D{d)x ,{d)A{d)W + l(C(d)Y f ) ■ D* mx ,{d)A{d)W 

+ S(D(d)Y') ■ D A{ff)W {d)C{d)X' - S{D{d)Y>) ■ D c{d)x ,(d)A(d)W = 0. 

Proof. By (6.9) and (4.7), we have 

{B(d)Z\[A(d)X,A(d)Y]) = \{B{d)Z) ■ D A{8)Y {d)A{d)X 
(6.34) - \{B{d)Z) ■ D A(8)x (d)A(8)Y = \(A(d)X) • j-{B{d)Z\A{d)Y) 
- \{A{d)Y) ■ D B{8)z (d)A(d)X - $(A(d)Y) ■ D A(d)x (d)B(d)Z . 

Similarly, we have 

{D(d)Z'\[A(d)X,C(d)Y']) = \(A(8)X) ■ ^{D{d)Z'\C{d)Y>) 

- \{C{d)Y>) ■ D mz> (d)A(d)X - l{C{d)Y') ■ D A{d)x {d)D{d)Z> , 

(B(d)W\[C(d)X',A(d)Y]) = \(C(8)X') • ^{B{d)W\A{d)Y) 



\{A{d)Y) ■ D B{d)w {d)C{d)X' - l{A{d)Y) ■ D* c(d)x ,(d)B(d)W , 



and 
(6.37) 



(D(d)W'\[C(d)X',C(d)Y']) = \(C(d)X') ■ ^-(D(8)W'\C(8)Y') 

ou 

- ${C{d)Y') ■ D D{s)w ,{d)C{d)X' - \{C{d)Y') ■ D* cm ,(d)D(d)W . 
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By the skewadnointness of H and K, which translates to (6.23), and by 
conditions (6.31), we have 

(B(d)Z\A(d)Y) = -{A{d)Z\B{d)Y) 

= -{C{d)Z'\D{d)Y') = {D{d)Z'\C{d)Y') , 

hence the first terms in the RHS of (6.34) and (6.35) cancel. Similarly for the 
first terms in the RHS of (6.36) and (6.37). Therefore, combining equations 
(6.34)-(6.37), we get that equation (6.32) is equivalent to 
(6.38) 

- \{A{d)Y) ■ D B{d)z (d)A(d)X - l(A(d)Y) ■ D* A{8)x {d)B{d)Z 

+ S(C(S)y') • D D{8)z ,{d)A{d)X + HC{d)Y') ■ D* A{8)x (d)D(d)Z' 
+ S(A(d)Y) • D B[d)w {d)C{d)X> + \(A(d)Y) • D* {8)x ,(d)B(d)W 

- m^Y') ■ D D{d)w ,{d)C{d)X' - \{C{d)Y ! ) ■ D* c{8)x ,(d)D(d)W' = . 

Next, since by assumption H = AB^ 1 is a non-local Poisson structure, it 
follows by Propositions 6.11 and 6.9 that equation (6.13) holds. In particu- 
lar, 

- $(A(d)Y) ■ D B(d)z (d)A(d)X - S(A(d)Y) ■ D A[g)x (d)B(d)Z 
(6.39) = -\(A(d)Y) ■ D B{d)x {d)A{d)Z + \(A(d)Y) • D* B{8)x (d)A(d)Z 
+ \(B(d)Y) • D A(d)z {d)A{d)X - l{B{d)Y) ■ D A{8)x {d)A{d)Z . 

Similarly, using the assumption that K = CD^ 1 is a non-local Poisson 

structure, we get 

(6.40) 

- J(C(5)F') • D mw ,{d)C{d)X> - \{C{d)Y>) ■ D* c{8)x ,{d)D{d)W 
= -\{C{d)Y') ■ D mx ,{d)C{d)W + l(C(d)Y f ) ■ D* mx ,{d)C{d)W 
+ S(D(d)Y') ■ D c(d)w ,{d)C{d)X' - \{D{d)Y') ■ D c{e)x ,{d)C{d)W' . 

Combining equations (6.38), (6.39) and (6.40), we get (6.33). □ 

Proof of Theorem 6.17. By Lemma 6.19, we only need to prove that, if con- 
dition (6.22) holds, then equation (6.33) holds for every X, X' ,Y,Y' ,W, Z' 
satisfying the first three identities in (6.31). It follows by some straightfor- 
ward computation that we can rewrite each term in the LHS of (6.33) as 
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follows 
(6.41) 

- $(A(d)Y) ■ D B(d)x (d)C(d)Z' = -$(A(d)Y) ■ B(d)D x (d)C(d)Z' 

f £ F^ fe (A+/,+5)^^(M + 5)^, y ( M )(| A=5 X i ,)(| M=5 4), 

J i' U^J OU\ 



j,k,el, neZ + 



(6.42) 

\(A(d)Y) ■ D* mxl (d)C(d)Z> = $(A(d)Y) ■ D x ,{d)D*{d)C{d)Z' 

f £ ^M* fe (A+/,+5)(-A- / .-5)"^^C jy ( / x)(| A ^)(| M= ^,) 

i',j',k'el du k 



j,k,eI,neZ + 



(6.43) 

K*W) • D cm ,{d)A[d)X = \(D(d)Y') ■ C(d)D z ,(d)A(d)X 

+ \ 2 ^'^( A +/^+ 5 )^^( A + 5 ) n ^'(A)(| A=a ^)(L =a 4) 



i',j',k'el 
i,k,el ,neZ+ 



(6.44) 

- KiJ(d)y) • ^ (5) x(5)c(5)z' = -K^(a)y) • a(5)l> x (5)c(5)z' 

-f 2 n^ fe %(A+ M +5)^^( M + 5)"C7 i/ ( M )(| A ^)(|^4) 



i',j',k'el va j 
j,k,eI,neZ + 



(6.45) 

KC(a)Y') • Z> D(a)z ,(d)A(d)X = KC(5)n • D{d)D z ,{d)A{d)X 
+ \ 2 ^^a(A+/^+5)^^(A + 5)M i HA)(| A=9 ^)(L =a 4 



i',j',k'el 
i,k,el ,neZ+ 



(6.46) 

KC(5)y') • D1 (5)X (5)B(5)W = K<W) • D* x (d)A*(d)B(d)W 



i',j',k'el ua k 
j,k,el ,neZ + 



61 



(6.47) 

${A(d)Y) ■ D B{d)w (d)C(d)X> = l(A(d)Y) ■ B(d)D w (d)C(d)X> 
+ f ^ n'^ fe (A+/i+5)^^(A + 5)"QHA)(| A= ^,)(| /i=5 ^), 

J i',j',k'el du i 



i,k,el ,neZ + 



(6.48) 

\(A(d)Y) ■ D* m ,(d)D(d)Z> = \(A(d)Y) ■ D x ,(d)C*(d)D(d)Z> 
+ f ^ n'^A+M+5)(-A-/x-5)"^^^ y ( / i)(| A= ^0(L=^ 



i',j',k'el ua k 
j,k,eI,neZ + 



(6.49) 

- KC(d)y') • D D(5) x' (d)A(d) W = -\{C{d)Y') ■ D(d)D x ,(d)A(d)W 



i',j',k'el vu j 
j,k,eI,neZ + 



(6.50) 

KC(5)y') • D* B{3)X {d)A(d)W = \{C(d)Y>) ■ D x {d)B\d)A(d)W 
+ f 2 ^C fc %(A + ^ + 5)(-A-/,-5)"^^%K^)(| A ^0(l^)> 

■V n> h'^T GUI 



i',j',k'el ua k 
j,k,el ,neZ + 



(6.51) 

KB(5)y) • D A(8)w (d)C(d)X' = $(B(3)Y) ■ A(d)D w (d)C(d)X' 



Wei 
i,k,el ,neZ + 



(6.52) 

- yp(d)r') • I> c(a) jp(5)A(d)W = -$(D(d)Y') ■ C(d)D x ,(d)A(d)W 

J ,•/ h'^J ou) 



i',j',k'el ua j 
j,k,el ,neZ+ 



It follows from the skewadjointness conditions (6.23) that the first term in 
the RHS of (6.41) cancels with the first term in the RHS of (6.44), the first 
term in the RHS of (6.42) cancels with the first term in the RHS of (6.48), 
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the first term in the RHS of (6.43) cancels with the first term in the RHS 
of (6.45), the first term in the RHS of (6.46) cancels with the first term in 
the RHS of (6.50), the first term in the RHS of (6.47) cancels with the first 
term in the RHS of (6.51), and the first term in the RHS of (6.49) cancels 
with the first term in the RHS of (6.52). Furthermore, combining the second 
terms of the RHS's of (6.47) and (6.51), we get, thanks to (6.24), 

f 2 {B kk id)Y k ,){u lX {u^u k } H } K {\ x ^D u/ (d)X^(\^ g B jf (d)W f ) . 

i',j',k'el 
i,j,kel 

Combining the second terms of the RHS's of (6.43) and (6.45), we get, 
thanks to (6.25), 

[ 2 (l) fcfe K5)n , 0{^A{^ M n fe } x }^(| A=a ^(5)X l ,)(|^^ y (5)Z J ',) . 

i',j',k'el 
i,j,k,el 

Combining the second terms of the RHS's of (6.41) and (6.44), we get, 
thanks to (6.26), 

i',j',k'el 
i,j,k,el 

Combining the second terms of the RHS's of (6.49) and (6.52), we get, 
thanks to (6.27), 

i',j',k'el 
i,j,k,el 

Combining the second terms of the RHS's of (6.46) and (6.50), we get, 
thanks to (6.28), 

- f 2 {D kkl {d)Y^){{u lX u 3 } HX+ ^u k } K {\ x=d B v ^ 

J i' ! j',k'el 
i,j,k,el 

Finally, combining the second terms of the RHS's of (6.42) and (6.48), we 
get, thanks to (6.29), 

- J] ( B kkid)Y k/ ){{u l xu j } Kx+ ^u k } H {\ x=d D ii/ {d)X' i ,)(\ fi=g D jjl (d)Z' jl ). 

t,j,k,el 
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Putting together all the above results, we conclude that the LHS of (6.33) 
is equal to 

J] {B{d)Y) k ({ui X {u jfi u k } H }K + {uix{u jfl u k } K }H 
i,j,kel 

~ {u jfl {UixU k } H }K ~ {Uj^{Ui\U k }K}H ~ {{Ui\Uj} Hx+IJ/ U k } K 

-{{u tX u 3 } Kx+ u k } H )(\ x ^B(d)X) i (\^ d B(d)W) j , 

which is zero by (6.22). □ 

In view of Theorem 6.17, we can translate Theorem 6.16 in terms of 
compatible non-local Poisson structures. 

Theorem 6.20. Let V be an algebra of differential functions in m, . . . ,ug, 
which is a domain. Let H,K e Mat^ x ^V(<9) be compatible non-local Poisson 
structures over V. Let H = AB -1 , and K = CD~ l be their minimal frac- 
tional decompositions (cf. Definition 2.12), with A,B,C,D e Mat^ x ^V[5], 
B, D non- degenerate (cf. Definition 2.11 and Remark 2.14)- Let Fq = 
B{d)Z, Fx = D(d)Z' = B(d)W, F 2 = D{d)W , with Z,Z',W,W e V® e , be 
such that 

(6.53) D{3)Z' = B(d)W , A{3)Z = C{d)Z' , A(d)W = C(d)W' , 
and Fq and F\ are closed, i.e. 

D%(d) = D Fo (8) , D* Fl (8) = D Fl (d). 

Then: 

(a) For all X,Ye V® e , such that D(d)X, D(d)Y e lm(B), we have 
(6.54) ■ C*{d)(D F2 (d) - D%{d))C{d)X = 0. 

(b) Lf we also assume that K is non- degenerate, then F 2 is closed. 

(c) Moreover F 2 is exact in any normal differential algebra extension V of 
V: F 2 = 5 -&, where J/ 2 e V/dV. 

Proof. By Theorem 6.17, Ca,b{K^) and Cc,d{K^) are compatible Dirac struc- 
tures over /C, the field of fractions of V. Recalling the expressions (6.30) of 
Mc,c and Nc,a f° r these Dirac structures, we get by Theorem 6.16 that 
equation (6.54) holds over JC, hence over V, proving (a). Let us prove part 
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(b). It is proved in [CDSK12b] that any two non-degenerate matrix differ- 
ential operators B{d) and D{d) have a right common multiple B{d)D\{d) = 
D(d)Bi(d), where Bi(d),Di(d) e Mat£ X £ IC[d] are non-degenerate. By clear- 
ing denominators, we can assume that B\{jS) and D\(d) have coefficients in 
V. Hence, if X,Y e Im^), we have D(d)X,D(d)Y e Im(J3). Therefore, by 
part (a) we have 

^G-Bt{d)C%d){D F2 {d) - D F2 {d))C{d)B 1 {d)F = 0, 

for all F,G e V®^. Since, by assumption, C{8) and B\(d) are non-degenerate, 
it follows that Dp 2 (d) = Dp 2 (d), as we wanted. Finally, part (c) follows by 

the fact that, under the assumption that V is normal, the variational complex 
is exact (see [BDSK09, Thm.3.2]). □ 

7 Hamiltonian equations associated to a non-local 
Poisson structure 

7.1 A simple linear algebra lemma 

Let U, V, W be vector spaces over F. 

Definition 7.1. Given linear maps A : U — » W and B : U — > V, we say 

that v e V and w e W are (^4, B)-associated (over £/), and we denote this by 

v <^' B > w or w <^' B > v, if there exists u e U such that v = Bu and io = Au. 

Let (• | •) : U x U — * G be a symmetric bi-additive form with values in 
an abelian group G, and, by abuse of notation, let also (• | •) : V x W — > G 
be a bi-additive form with values in G. Given a subspace V\ a V, we define 
its orthogonal complement Vi~ a W as 

Vi = {w e W | (v\w) = for all v e Vi) 

and similarly, given a subspace W\ a W, we define its orthogonal comple- 
ment Wi a V as 

Wt = {veV\ (y\w) = for all w e Wi} . 

We say that linear maps A : U — > W and A* : V —> U are adjoint if we 
have 

(t>|Au) = , 
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for every u £ U and v e V, and similarly we say that linear maps B : U — > V 
and B* : W — > f7 are adjoint if we have 

for every u £ U and io e W. (The adjoints A* and B* are unique if the 
inner product (• | •) : U x U — * G is non-degenerate.) 

Lemma 7.2. Let A : f7 — VF, S : 17 — V, C : U ^ W , D : U — V ', 6e 

/inear maps, and Zei A* : V U , B* : W -+ U , C* : V — U , D* : W -+ U 
be their adjoint. Assume that 

(7.1) A*B + B*A = , C*D + D*C = 0. 

Let {i>n}nL_i c V" and {w; n }^L cz W be finite sequences such that 

r» c\\ ( A > B ) (C,D) 

(7.2) w n _i W Wn W Un , 
/io/ds /or every n = 0, . . . , N. 

(a) Then we have (v m \w n ) = for every m = —1, . . . , N, n = 0, . . . , N. 

(b) Suppose, moreover, that the the following orthogonality conditions hold: 

(7.3) ^Span F {?; n }^ = _ 1 ^ clmC, ( Span F {?z/ n }^L ) cImJ3. 

Then, we can extend the given finite sequences to infinite sequences 
{v n }n=—i c V > { w n}n=o ^ such that the association relations (7.2) 
hold for every n e Z + . 

Proof. By assumption, for every n = 0, . . . , N there exist u n ,u' n e U such 
that 

v n -i = Bu n , w n = Au n , v n = Du' n , w n = Cu' n . 

Hence, by definition of adjoint operators and assumption (7.1), we have, for 
every m,n, 

(%K) = ( Du m\ Cu n) = ( u 'm\ D *CK) = -(u' m \C*Du' n ) 
= -(C*Du' n \u' m ) = -(Du'JCu'J = -{v n \w m ) , 

and similarly 

(%|%) = (Bu m+1 \Au n ) = {u m+ i\B* Au n ) = -(u m+ i\A* Bu n ) 
= -{A*Bu n \u m+ i) = -(Bu n \Au m+1 ) = -(v n -x\w m+1 ) . 
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Hence, 



(7.4) (v m \w n ) = -{v n \w m ) = -(v n -i\w m+ x) . 

Letting m = n in equation (7.4) we get {v n \w n ) = for every n = 0, . . . , N, 
while letting m = n — 1 in (7.4) we get (u n _i \w n ) = for every n = 0, . . . , N. 
Equations (7.4) imply {v m \w n ) = (v m+ i\w n ~i), and therefore, by induction 
on n - ra, we get that (v m \w n ) = for n ^ m. On the other hand, by the 
first identity in equation (7.4) it follows that (v m \w n ) = also for n < m. 
This proves part (a). 

By part (a), we have that (v]y\w n ) = for every n = 0, ...,N, and 
therefore by the second orthogonality condition (7.3) we get that vn = 
Bun+i for some un + \ e U. We then let wn+i = Aujy+i and we get, by 

construction, that vjy ^' B > w^+i- By the same argument as in the proof of 
part (a), we have that (v n \wN+i) = for every n = —1, . . . , N, and therefore 
by the first orthogonality condition (7.3) we get that wn + i = Cu' N+1 for 
some u' N+1 e U. We then let vjy + i = Du' N+1 and we get, by construction, 

that vn + i wn+i- Hence, we prolonged the original finite sequences 

{v n }n=-i and { 

w n}n=o by one step- The claim follows by induction. □ 



7.2 Hamiltonian functionals and vector fields, and Poisson 
bracket 

This section serves as a motivation to the Lenard-Magri scheme of inte- 
grability, discussed in the following sections. We assume that the algebra of 
differential functions V is a domain, and we denote by /C its field of fractions. 

Let H e Mat£ X iV(d) be a non-local Poisson structure over V. Recall 
that, since H has rational entries, it admits fractional decomposition H = 
AB~ l , where A,B e Mat ixi V[d] and B is non-degenerate (cf. Definitions 
2.11 and 2.12). 

Definition 7.3. Elements §h e K/dK and P e /C^ are H-associated, and we 
denote this by J/t * — * P, if 

(7.5) S ^ = B(d)F, P = A(d)F, 

for some fractional decomposition H = AB^ 1 , with A,Be Mati x iV[d] and 
B non-degenerate, and some element F e /C . In this case, we say that $/ 
is Hamiltonian functional for H, and P is a Hamiltonian vector field for H. 
We denote by F{H) a K/dlC the subset of all Hamiltonian functionals for 
H, and by %{H) a K, the subset of all Hamiltonian vector fields. 
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Remark 7.4. Note that in the definition (7.3) we can fix a minimal fractional 
decomposition H = A\B^ 1 for H, with Ai,B\ e M.&t£ X gV[d] and B\ non- 
degenerate of minimal possible order. Indeed, by Proposition 2.13(b) if 
H = AB^ 1 is any other fractional decomposition, then there exists D e 
Mati X ilC[d] such that A = AiD, B = B\D. Therefore, if equations (7.5) 
hold for some F e K® 1 , then we have ^ = -Bi(d)Fi, P = A±(d)F 1 , where 
F\ = D(d)F e K® . It follows that F{H) and H{H) are vector spaces over 
F. In fact, they are given by the following formulas 

HH) = (j^)~\ B A^)) c K{H) = A^^j-K/dK)) cz . 

Remark 7.5. Consider Definition 7.1 with U = V = K® e and W = /C £ . 
Comparing this with Definition 7.3, we have that \h < — > P if and only if 

< A ' >^ P for some fractional decomposition = AB^ 1 . 

In the local case, when e Mat^ x ^V[5] is a (local) Poisson structure 

over V, then J/ e JF(H) = K/dK and P e %{H) = H(d)[lmj^j cz K e are 

associated if and only if P = H(d)^-. 

Lemma 7.6. (a) If If A P and\g A Q, then\{af+bg) A {aP+bQ) 
for every a,b e C (the subfield of constants in K). In particular, F{H) 
and H(H) are vector spaces over C. 



(b) If J/ A P, then I \g e F(H) \g A P\ = J/ + F (H), where 



(7.6) Fo{H) = [\geF{H) \g A J 

(c) If If A P, then \Q e U{H) J/ A q\ = P + H (H), where 
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(7.7) H (H) = [QeU{H) 

Proof. Obvious, using Remark 7.4. □ 

Lemma 7.7. (a) The space K is a lie algebra with bracket (6.9), and 
H{H) cz K is its subalgebra. 

(b) We have a representation <fi of the lie algebra K on the space K/dK 
given by 
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and the subspace T(H) a K/dK is preserved by the action of the Lie 
subalgebra T-L(H) a K . 

(c) If \h < — ► P and \h < — > Q for some §h e P(H), then the action of 
P,Q e 7i(H)on J-{H) is the same: 

Proof. It follows immediately from [BDSK09, Lem.4.7-8], using the fact that 
C-a,b{K>) is a Dirac structure if H = AB^ 1 is a minimal fractional decom- 
position for H. □ 

Thanks to Lemma 7.7, we have a well-defined map {• , -}h ■ J~(H) x 
P(H) T{H) given by 
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where PeU{H) is such that \f p 



Proposition 7.8. (a) The bracket (7.8) is a Lie algebra bracket on the 
space of Hamiltonian functionals P{H) . 

(b) The Lie algebra action of rl{H) on P{H) is by derivations of the Lie 
bracket (7.8). 

(c) The subspace 

A(H) = {G/,P) e JF{H) x H(H) | J/ ♦ p] 

is a subalgebra of the direct product of Lie algebras P{H) x T-L{H). 

Proof. It follows immediately from [BDSK09, Prop. 4. 9, Rem. 4. 6], using the 
fact that Ca,b{K-) is a Dirac structure. □ 

7.3 Hamiltonian equations and integrability 

Let V be an algebra of differential functions. We have a non-degenerate 
pairing (• | •) : V x V® £ — > V/dV given by 

(7.9) (P|£) = $P.£. 

(See e.g. [BDSK09] for a proof of non-degeneracy of this form.) 
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Let H e Mat£ X ^V(5) be a non-local Poisson structure. If H = AB^ 1 
is a fractional decomposition of H, with A, B e Mat£ X £V[5] and B non- 
degenerate, according to Definition 7.1 with U = V = V®^ and W = V , we 
have that £ e V® f and Fe V* are (j4, P)-associated, 

(7.10) £ ( ^3p, 

if there exists F e V® £ such that £ = BF, P = AF. 

Let $/i e V/5V and P e V e be such that fjj P for some fractional 
decomposition H = AB~ l . The corresponding Hamiltonian equation is, by 
definition, the following evolution equation on the variables u = {ui) ieI - 

By the chain rule, any element / e V evolves according to the equation 

al ie/neZ+ 5«i 

and, integrating by parts, a local functional J/ £ V/5V evolves according to 

$L.[p.!l (=(p\ s -l)). 

dt J 6u \ y l 5u J J 



Definition 7.9. The Hamiltonian equation (7.11) is said to be integrable if 
there exist sequences {Cn}neZ + c V®^ and {P n }neZ + ^ V such that: 

(i) elements £ n 's and P n 's span infinite dimensional (over the subalgebra 
C a V of constants) subspaces of V®^ and respectively; 

(ii) for every n £ Z + we have the association relation £ n <^' B > P n , for some 
fractional decomposition P = AP -1 ; 

(iii) the elements £ n 's are closed, i.e. they have self-adjoint Frechet deriva- 
tives: D &i (d) = Dl(d); 

(iv) the elements P n 's commute with respect the the Lie bracket (6.9): 
[P m , P n ] = for all m,n e Z + ; 

(v) (P m | £ n ) = for all m,n e Z + . 
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In this case, we have an integrable hierarchy of Hamiltonian equations 

da 



dt 



P n , ne 



n 



Remark 7.10. Recall from [BDSK09, Prop. 1.5] that if £ e V® e is closed, then 
it is exact in any normal algebra of differential function extension V of V: 
£ = Tu f° r some \h E V/dV. If V is a domain, then by Lemma 4.3 it can 
be extended to a normal algebra of differential functions V which is still a 
domain, and we can consider its field of fractions /C. Then we have £ n = 
where {§h n } ne z cz F{H) cz IC/dJC form an infinite sequence of Hamiltonian 
functionals in involution: {§h m ,^h n } h = for every m, n e Z + (cf. Section 
7.2). 

In analogy to Liouville integrability of finite dimensional Hamiltonian 
systems, we should require in addition some completeness property of the 
span H cz V®^ of the variational derivatives of the conserved densities £ n = 
^ IL , n e Z+, and of the span n of the generalized symmetries P n , n e Z+. A 
natural condition, analogous to Liouville integrability, in the general setup 
of non-local Poisson structures, is the following. 

Definition 7.11. A completely integrable system associated to the non-local 
Poisson structure H = AB^ 1 , in its minimal fractional decomposition, is a 
pair of subspaces 5 = B(U) cz V® £ and n = A(U) cz V , for some subspace 
U cz V® e , such that 

(i) S consists of closed elements in V®^: D^(d) = D^(d) for all (eS; 

(ii) n is an abelian subalgebra of V with respect to the Lie bracket (6.9); 

(iii) n 1 = 3 and H 1 - = n with respect to the pairing (7.9). 

In this case, for every Pellwe get a completely integrable Hamiltonian 
equation ^jr = P, and all local functionals $/ e V/dV such that |^ e H are 
its integrals of motion in involution. 

Remark 7.12. In the local case we let E be the span of £ n = n e Z + , 
where h n are the conserved densities, and we let n = H(E). Then the 
above condition H 1 - = S is equivalent to the condition that H is a maximal 
isotropic subspace of = V®^ with respect to the skewsymmetric bilinear 
form Q\ x — > V/5V given by = (-£f£|??). Indeed, £ e fii satisfies 

(£|£ n > = -(f|JTf n ) = for all n if and only if £ ± H(E) = II. In this 
case, the J/i n 's are automatically in involution and n consists of commuting 
higher symmetries. 
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Remark 7.13. We can generalize Definition 7.11 of complete integrability 
to the case of an arbitrary Dirac structure C a V® 1 © V e as a subspace 
A a C such that S = vri(A) a V® 1 and II = ^(A) a V , the projections 
of A c V® £ © in the first and second components respectively, satisfy 
conditions (i)-(iii) above. These conditions are equivalent to require that 
7Ti (A) consists of closed elements, that Ao = vri(A) ©7r2(A) is a maximal 
isotropic subspace of V® e © V with respect to the symmetric bilinear form 
(• | •) defined in (6.7), and the Courant-Dorfman product o defined in (6.8) 
restricted to Ao is zero. In other words, Ao is a Dirac structure with zero 
Courant-Dorfman product. 

Example 7.14. It is not hard to check, using arguments similar to those in 
[BDSK09], that the KdV equation is completely integrable in the sense of 
Definition 7.11. 



7.4 The Lenard-Magri scheme of integrability 

Theorem 7.15. Let V be an algebra of differential functions, and let H = 
AB^ 1 and K = CD^ 1 be rational skewadjoint pseudodifferential operators, 
with A,B,C,D e Mat^ x ^V[5] and B,D non- degenerate. Let {£n}^L_i a 
V® , {P n )n=o ^ ^ e fi n tt e sequences such that 

(7.12) <- L -^ P n <- L -^ £ n , 

holds for every n = 0, . . . , N. Then: 
(a) We have 

(7.13) (P n \Cm) = 0, m>-l,n>0. 



(b) Assume that V is a domain and H and K are compatible non-local Pois- 
son structures with K non- degenerate, and assume that £_i and £o are 
closed, i.e. their Frechet derivatives are selfadjoint. Then the elements 
£ n , n ^ 1, are closed as well, and we have 

(7.14) [? m ,P n ]eKerB*nKerZ)* , m, n ^ . 

(c) Assume that the matrices A,B,C,D have non- degenerate leading coef- 
ficients, and that 

dord(P n ) > max { dord(A) - \H\ + \K\,dord(B) + \K\, 

(7.15) 1 , 

dord(C , dord(D) + \K\L 
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for some n > 0. Then 

dord(£ n ) = dord(P n ) - \K\ , dord(P n+1 ) = dord(P n ) + \H\ - \K\ . 
In particular, if \H | ^ \K\, then 

dord(P i ) = dord(P n ) + (j - n)(\H\ - \K\) = dord(^) + \K\ , 
for every j ^ n. 
(d) Assume that the following orthogonality conditions hold: 

(7.16) (Span c {£ m }£ = _ 1 ) ± cImC' , f Spane^}^) 1 c ImB , 



where the orthogonal complements are with respect to the pairing de- 
fined in (7.9). Then we can extend the given finite sequences to infinite 
sequences {£m}m=-i cz V®^ ; {P n }n=o c V sitc/i i/iai i/ie association 
relations (7.12) /ioW /or every n e Z+. 

Corollary 7.16. Let V be an algebra of differential functions, which is a 
domain, let H = AB~ 1 and K = CD^ 1 be compatible non-local Poisson 
structures, where A,B,C,D are £ x t matrix differential operators with non- 
degenerate leading coefficients and such that \H\ > \K\. Let {£ n }n=— l c 
V® e , {P n }n=o c * be fi n ^ e sequences such that £_i and £o ar ^ closed, 
conditions (7.12) hold for every < n < N , condition (7.15) holds for 
some < n < N , and the orthogonality conditions (7.16) hold. Then the 
given finite sequences can be extended to infinite sequences {£n}n=-i cz V®^, 
{P n }" =0 cz V , such that the differential orders of the £ n 's and the P n 's tend 
to infinity as n — > oo, all £ n 7 s are closed, and equations (7.12), (7.13) and 
(7.14) hold. Consequently, 

is an integrable bi-Hamiltonian equation for every n e Z+. //, moreover, 
Ker B* n KerL>* = 0, all equations (7.17) /orm a (compatible) integrable 
hierarchy of bi-Hamiltonian equations. 

Proof of Theorem 7.15. Parts (a) and (d) are special cases of Lemma 7.2(a) 
and (b) respectively, since the assumption that H and K are skewadjoint is 
equivalent to equations (7.1). Part (b) follows Lemmas 7.17 and 7.18 below. 
Finally, part (c) follows form Lemma 7.19 below, with £ = £ n , P = P n and 

Q = Pn+l- □ 
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Lemma 7.17. Let H = AB^ 1 and K = CD^ 1 be compatible non-local 
Poisson structures, with K non- degenerate, over the algebra of differential 
functions V, which is a domain. Let £o>£i be closed elements of V®^ , £2 £ 
V® e , and Pi, P 2 e V e be such that 

(7.18) eo^V^V^V^. 

Then £2 is closed. 

Proof. If H = AB -1 and K = CD~ l are minimal fractional decompositions, 
then the statement follows from Theorem 6.20(b). Indeed, conditions (7.18) 
imply the existence of Z,Z',W,W e V® e such that £ = B(d)Z, £1 = 
D{d)Z' = B{d)W, £ 2 = D(d)W, and solving equations (6.53). 

In general, the fractional decompositions H = AB~ l and K = CD -1 
are not necessarily minimal. Let fC be the field of fractions of V. By Propo- 
sition 2.13 we have A = A\P, B = B X P, C = CiQ, D = D X Q, with 
Ai,Pi,Ci,Z)i,P,Q e M&t M K.[d], where H = AiB^ 1 and K = CiPf 1 
are minimal fractional decompositions, and P, Q are non-degenerate. Ob- 
viously, by the definition (7.10) of (A, P)-association, if £ < A '~ B > P holds over 

V, in the sense that £ = BF, P = AF for some F e V® £ , then £ P 
holds over /C, indeed £ = P^Pi, P = AiF u where P x = PP e /C®*. Hence, 
conditions (7.18) hold (over /C) with A,B,C,D replaced by A\, B%, C±, D\. 
Then, by Theorem 6.20(b) we get that £2 is closed over /C, hence over V. □ 

Lemma 7.18. Let H = AB^ 1 and K = CD^ 1 be compatible non-local 
Poisson structures, with K non- degenerate, over the algebra of differential 
functions V, which is a domain. Let {£n}^L_i be closed elements of V® £ , 
and {P n }n =0 be elements ofV 1 satisfying conditions (7.12). Then 

[P m ,P n ] e KerP* n KerP* , m,n^0. 

Proof. Let JC be the field of fractions of V, and let H = AiPf 1 , K = CiPf 1 
be their minimal fractional decomposition over /C. As observed in the proof 
of Lemma 7.17, all association relations (7.12) hold, over /C, after replacing 
A,B,C,D with A\,Bi,C\,D\ respectively. By Theorem 6.12, £ai,Bi(^C) 
and Canity are Dirac structures in /C®^ © K. , in particular they are 
closed with respect to the Courant-Dorfman product (6.8). By the def- 
inition (6.11) of the Dirac structures Ca 1 ,b 1 (K-) and £>C\,Di{)€)> we have 
€n-i © Pn e Ca 1 ,b 1 (1C) and £ n © P„ e Cc^dA^) for every n > 0. By the 
assumption that all the £ n 's are closed, we can use formula (6.10) to de- 
duce that ^(P m I £ re _i) © [P m , P n ] e £ AljBl (/C) and ^(P m | £„) © [P m , P n ] e 
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£>At,B\(K>)- Hence, by Theorem 7.15(a) we conclude that 0© [P m ,Pn] £ 
£ai,£i(£) n ^c 1) x) 1 (/C). Namely, there exist F\ e KerSi cz K and G\ e 
KerDi cz /C £ such that [P m ,P n ] = A^Fi = C x {d)G e Ai(KerSi) n 
Ci(KerL'i). By skewadjointness of -ff and ET, we have B*A\ = — A*B\ 
and D*C\ = —C*D\, which immediately implies Ai(Ker B\) cz Ker B* and 
Ci(KerL»i) cz KerDJ. Therefore, [P m ,P n ] e KerBfnKevDf cz K l . On the 
other hand, since B and -D are right multiples of B\ and D\ respectively, 
we have Kerl?* n V f c Ker B* and KerD* n V* c KerD*. Thererore, 
[P m , P n ] e Ker B* n Ker D* cz V f , as we wanted. □ 

Lemma 7.19. Let V be an algebra of differential functions, let A, B,C,D e 
Mat£ X ^V[5] be matrices with non- degenerate leading coefficients. Denote by 
H = AB^ 1 and K = CD^ 1 the corresponding rational matrix pseudodiffer- 
ential operators. Let P,QeV e and £ e V®^ satisfy the following association 
relations (cf. Definition 7.1) 

(7.i9) p ( ^e ( ^Q, 

and assume that 

dord(P) > max{dord(A) -\H\ + |i^|,dord(B) + \K\, 



(7.20) 

dord(C),dord(L>) + \K\}. 

(Here we use the notation introduced in (4.8) and (4.9) J Then 

dord£ = dord(P) - \K\ and dord(Q) = dord(P) + \H\ - \K\ . 

Proof. By definition, the relations (7.19) amount to the existence of elements 
F, G e V e such that 

CG = P , DG = £ , BF = £ , AF = Q . 

Since C has non-degenerate leading coefficient and dord(P) = dord(CCr) > 
dord(C), we get by Lemma 4.7(c) that 

dord(G) = dord(CG) - |C| = dord(P) - |C| . 

Next, since by assumption D has non-degenerate leading coefficient and 
dord(G) + \D\ = dord(P) - |C| + \D\ = dord(P) - \K\ > dord(D), we get 
by Lemma 4.7(b) that 

dord(£) = dord(DG) = dord(G) + \D\ = dord(P) - |C| + \D\ 
= dord(P) -\K\. 
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Similarly, since, by assumption, B has non-degenerate leading coefficient 
and doxA{BF) = dord(£) = dord(P) — \K\ > dord(-B), we get by Lemma 
4.7(c) that 

dord(F) = dord(PF) - \B\ = dord(£) - \B\ = dord(P) - \K\ - \B\ . 

Finally, since, by assumption, A has non-degenerate leading coefficient and 
dord(F) + |A| = dord(P)-|K| = \B\ + \A\ = dord(P) -\K\ + \H\ > dord(^), 
we get by Lemma 4.7(b) that 

dord(Q) = dord(AF) = dord(F) + \A\ = dord(P) - \K\ - \B\ + \A\ 
= dord(P) -\K\ + \H\. 

□ 

Proof of Corollary 7.16. The statement of the corollary is basically a sum- 
mary of parts (a)-(d) of Theorem 7.15, except that we need to explain why, 
by the condition (7.14), it follows that each P n lies in an infinite dimen- 
sional abelian subalgebra contained in Span{P n }^ =0 . This follows from the 
observation that Ker(i?*) n Ker(D*) is finite dimensional over C, and the 
following result. 

Lemma 7.20. Let U be an infinite dimensional subspace of a Lie algebra 
such that [U, U] is finite dimensional. Then any element of U is contained 
in an infinite dimensional abelian subalgebra of U . 

Proof. Let a\ be a non-zero element of U. The centralizer C\ of a\ in U is 
the kernel of the map ad a : U — » [U, U], hence, it has finite codimension in 
U. Next, let 0,2 be an element of C\ linearly independent of a±, and let C2 
be its centralizer in C\. By the same argument, C2 has finite codimension in 
C\. In this fashion we construct an infinite sequence of linearly independent 
commuting elements of U. □ 

□ 

Remark 7.21. Suppose that the sequences {^ n }" = _ 1 cz V® e and {P n }^ =0 
satisfy relations (7.12) for each n e Z + with respect to the compatible non- 
local Poisson structures H = AB^ 1 and K = CD^ 1 , and assume that 
their spans S = Span{£ n } nS ._i and IT = Span{P n } define a completely in- 
tegrable system in the sense of Definition 7.11. Then the orthogonality 
conditions (7.16) automatically hold for some N (possibly infinite). Indeed, 
by the relations (7.12) we have, in particular, that 5 cz Im(i?) n im(D) and 
II cz Im(.A) n Im(C). Therefore conditions (7.16) follow by axiom (hi) in 
Definition 7.11. 
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Remark 7.22. Unfortunately we are unable to prove a stronger form of equa- 
tion (7.14), namely that the generalized symmetries P n obtained by the 
Lenard-Magri scheme commute. The usual "proof" of this fact using the 
recursion operator (see e.g. [Bla98, p. 64] or [01v93, p. 317]) is not rigorous. 

The recurrence relations (7.12) are usually represented by the following 
diagram, called the Lenard-Magri scheme: 
(7.21) 

P Pi P2 

(A,B) /• \ (C,D) (A,B)/^ \ (C,D) (A,B)/^ \ (C,D) 

£-1 

Explicitly, diagram (7.21) holds if there exists a sequence {F n } n> _! in V® 1 

such that the following equations hold (n e Z + ): 

(7.22) 

B{d)F- X = £_ x , C(8)F 2n = A(d)F 2n _ 1 = P n , B(d)F 2n+1 = D(d)F 2n = £ n . 

7.5 Notation, terminology and assumptions 

In the following sections we apply the machinery developed so far in explicit 
examples. 

As stated in Theorem 7.15, if the algebra V is a domain, the coordinates 
of all £ n 's and P n 's solving the recurrence equation (7.12) lie in V. However, 
for notational purposes it is convenient to go to the field of fractions of V, 
so we will assume that V is a field. Under this assumption, the notation 

$/i *—* P introduced in Definition 7.3 is consistent with the notation £ <^' B > 
P in (7.10). Namely, the following conditions are equivalent: 

(i) $h <*+ P, 

(ii) xjj <-^— > P for some fractional decomposition H = AB^ 1 , 

(iii) |^ ^ ,B > P for the minimal fractional decomposition H = AB^ 1 . 
Hence, in the rest of the paper we will use the more suggestive notation 

Furthermore, in all examples we will begin the Lenard-Magri scheme 
with £_i = 0. In this case, by Lemma 7.17, all the elements £ n are closed, 
provided that £o is closed. 
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Recall also that if all £ n 's are closed, they are exact in some algebra of 
differential function extension of V: £ n = ^f. In this case, and using the 
more suggestive notation above, the diagram (7.21) takes the form 

(7 .23) JO JU P $h ^J\^J/n^... 

and it is equivalent to the existence of {-F n }n>-i in V®^ such that the fol- 
lowing equations hold (n e Z+): 
(7.24) 

B(8)F- 1 = 0, C(d)F 2n = A(d)F 2n - 1 = P n , B(d)F 2n+1 = D(d)F 2n = -± . 

ou 

Remark 7.23. In general, if V is an arbitrary algebra of differential functions, 
{Un=-i <= V© £ , {P„}£L c V«, and A,B,C,D e MaWV[d], then, the 
whole infinite sequences {£ n } n j>_i and {-P n }n^0i constructed using Theorem 
7.15, have coordinates in V. 



S-type vs C-type Lenard-Magri schemes 

Consider a Lenard-Magri scheme as in (7.23). We say that it is finite if it 
can be extended indefinitely, but in any such infinite extension the linear 
span of {\h n } ne i + or of {P n }neZ + is finite dimensional. We say that the 
Lenard-Magri scheme (7.23) is blocked if it cannot be extended indefinitely, 

namely, for some n, there is no §h n such that P n * — * \h n , or there is no 

P n+ i such that \h n < — > P n +i- 

For an integrable Lenard-Magri scheme (7.23), we say that it is of S-type 
if the differential orders of the elements P n grow to infinity, and it is of C-type 
if the differential orders of the P n 7 s are bounded. It is easy to see that for an 
integrable Lenard-Magri scheme of S-type the order of the pseudodifferential 
H should be greater than the order of the pseudodifferential K. Indeed, since 
we have P n < — ► \h n < — > P n +ii if Pn has differential order large enough, 
then, by Lemma 7.19, dord(P n+ i) = dord(P n ) + ord(-ff) — ovd(K). 

Remark 7.24. This terminology in inspired by the terminology of Calogero, 
who calls an integrable hierarchy of "S-type" if the differential orders of 
the canonical conserved densities are unbounded, and of "C-type" otherwise 
(see [MSS90, MS12]). Note that, though these two terminologies are close, 
they do not coincide. For example, the linear equation 4| = u'" is C- 
integrable in Calogero's terminology, but the corresponding Lenard-Magri 
scheme, considered for example in [BDSK09], is integrable of S-type. 
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8 Liouville type integrable systems 



In this section V is a field of differential functions in u, and we assume that 
V contains all the functions that we encounter in our computations. As 
before, we denote by C c V the subfield of constants, and by T cz V the 
subfield of quasiconstants. We shall denote by x the element of T such that 
dx = 1. 

Recall from Example 4.13 that we have the following triple of compatible 
non-local Poisson structures: 

L t = d, L 2 = d- 1 , L 3 = u'd- 1 o v! . 

Given two non-local Poisson structures H and K of the form 

(8.1) H = a\L\ + a 2 L 2 + a 3 L 3 , K = b\L\ + b 2 L 2 + b 3 L 3 , 

with a,i,bi e C, i = 1, 2, 3, we want to discuss the integr ability of the corre- 
sponding Lenard-Magri scheme. 

8.1 Preliminary computations 

First, we find a minimal fractional decomposition for the operators H and 
K. 

Lemma 8.1. For X = x\L\ + x 2 L 2 + x 3 L 3 , with x\,x 2 ,x 3 e C, we have 



.2) X 



>\2 



_n 1 X 2 + X 3 (u') . 

x\S o —o H o — x 3 u 



u" 



-1 



The above fractional decomposition is minimal only for x 2 x 3 ^ 0. For x 2 =£ 
0, x 3 = 0, the minimal fractional decomposition is 

(8.3) X = ( Xl d 2 + x 2 )d~ 1 , 
for x 2 = 0, X3 ^ it is 

(8.4) J=( Sl aoIa + S3«')(^) _1 , 
and for x 2 = x 3 = it is X = x±d. 

Proof. Straightforward. □ 

Later we will need the following simple facts concerning the numerators 
of the fractional decompositions for X. 
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Lemma 8.2. (a) For xi,X2, X3 e C ; 11 / 0, consider the equation 

(8.5) xi5o— 5+ 0-s 3 u)F = /, 

in F e V and f e T . If X3 ¥= 0, then all the solutions of equation (8.5) 
are 

F = au , f = x^ol for some a e C , 
while if X3 = 0, then all the solutions of equation (8.5) are 

F = au 1 + f3(xu' — u) + 7 , f = X2(a + fix) for some a, j3,j e C . 

(b) For x\,X2 £ C, x± =£ 0, an element F e V satisfies 

(8.6) (xid 2 + x 2 )F e T 
if and only if F e T . 

(c) For xi, X3 e C, x\ 7^ 0, an element F e V satisfies 

(8.7) (xido —d + x 3 u'^Fe Vi 

if and only if F e Vo and F 1 = . 

Proof. If n ^ 2 and F e V solves equation (8.5) and has differential order 
less than or equal to n, then, using (4.1), we have 

which implies that ^J n ^ = 0. Hence F must have differential order at most 
1. Then we have 

= J—LHS{S.h) = n - — 1-^' 
^0,(4) v 7 1 v -11" ;w r,/'A2 



dit( 4 ) W" <9u' [W t 

so that = j^jru" . But then equation (8.5) becomes 

/ f)F\" f)F 
(8.8) Xl [— ) +(x 2 +x 3 (n') 2 )^ 7 -^ / F = /. 

If has differential order n ^ 0, then applying g ^ ( f +2) to both sides of 
equation (8.8) we get Sl f( n fg u / = 0- Hence, |^ = ip e T . In other words, 
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F = ipu' + fo, where /o e Vo has differential order less than or equal to 0. 
But then the condition F' = §^u" becomes (p'u' + /q = 0. This implies, 
using (4.3), that (/o + (p'u)' = ip"u e dV n Vo = dF. So, necessarily, y/' = 0. 
Hence, if = a + (3x and fo = — /3u+7, for some constants a, (3, 7 e C. Putting 
these results together, we have 

F = {a + fix)u' - /3u + 7 , 

and plugging back into equation (8.8) we get 

x 2 (a + fix) + x 3 j3uu' — xsju' = f . 

Since, by assumption, / e F, we obtain f3 = 7 = if £3 ^ 0, completing the 
proof of part (a). 

For part (b) we just observe that, if F e V n for some n ^ satisfies 
condition (8.6), then 

3 OF 
0= w^t(^" + ^)=x 1 - 



Hence, i 7 must be a quasiconstant. 

Similarly, if F e V n for some n ^ 1 satisfies condition (8.7), then 

5 / F' ,\ xi dF 



Hence, F must lie in Vo- Furthermore, 

8 ( „F' ,„\ ( F' 1 dF\ 

= d^[ Xld V + X3uF ) =Xl {-W 2 + ^r 

Hence, F must be such that F' = W-u' . □ 
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Next, we compute the spaces Fq(X) and Ho(X) defined in (7.6) and 
(7.7). Here and further we use the following notation: given two constants 
Xi, xj e C such that Xi =£ 0, we let 



(8.9) Xij = J-^ eC. 



(We assume that the field C contains all such elements.) 
Lemma 8.3. For X = x\L\ + X2L2 + X3L3, we have: 
(a) F (X) = Ker (£) if Xl x 2 x 3 * 0; 
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MX) 


= C<\e Xl2X u + C\eT x ^ x u + Ker (^) if x x x 2 + 0, x 3 = 


MX) 


= C$e Xl < iU + C$e~ Xl3U + Ker (^) if Xl x 3 * 0, x 2 = 0; 


MX) 


= C\u + Ker (^-) if x\ and x 2 = x 3 = 0; 


MX) 


= C\«Jx 2 + x 3 (u') 2 + Ker (^) , if x x = §,x 2 x 3 / 0; 


MX) 


= Ker (Jjj) if x\ = and x 2 = or x 3 = 0. 


(b) n (x) 


= C © Cu' i/ X2X3 =£ 0; 


H (X) 


= C if x 2 ^ 0,x 3 = 0; 


n (x) 


= Cu' if x 2 = 0, x 3 0; 


H (X) 


= if x 2 = x 3 = 0. 



Proof. First, let us find all elements P e T-Lq(X). By Remark 7.4, if X = 
YZ -1 is a minimal fractional decomposition, we need to solve the following 
equations in F, P e V: 

(8.10) ZP = 0, P = YP. 

By Lemma 8.1, if x 2 = x 3 = 0, then Y = x\d and Z = 1, so the only solution 
of (8.10) is given by F = 0, P = 0. If x 2 ^ 0, x 3 = 0, then Y = x x d 2 + x 2 
and Z = d, so we get P e C and P e C. Similarly, if X2 = 0, x 3 =/= 0, 
then y = x\d o ^-<9 + 2:31/ and Z = ^5, so we get F e C and P e Cu'. 
Finally, if x 2 =£ 0, x 3 ^ 0, then Y = xi<9 2 o ±d + X2+ ^(»') 2 g _ ^ and 
Z = d o ^pd. Hence, the solutions of (8.10) are F = a + flu' e C ©Cm', and 
P = YP = x 2 (3 — x 3 au' e Cu' . This proves part (b). 

Next, we find all elements J/ £ -Po(^O) namely all solutions of the fol- 
lowing equations in P e V and $/ £ V/5V: 

(8.11) yp = o, 6 J- = zf. 

ou 

If xi = 0, £2 ^ 0, X3 = 0, we have Y = x 2 is invertible, and similarly, if 
x\ = 0, X2 = 0, X3 ^ 0, we have Y = x 3 v! is invertible too. In both these 
cases we thus have P = 0, and hence \f e Ker (^) . If x\ = 0, x 2 + 0, x 3 i= 0, 
then Y = X2+x ^ 2 d - x 3 u' and Z = d o ^<5. The equation YP = has a 
one-dimensional (over C) space of solution, spanned by P = \jx 2 + x 3 (u') 2 . 
Hence, all elements §f £ ^(X) are obtained solving the equation 

5 J- = ado \d^/x 2 + x 3 {u') 2 = a( XiU =) , 
ou u" y -\/x 2 + x 3 (u') 2) 
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for a £ C. Its solutions are of the form J/ = —a\/x2 + x 3 {u') 2 + k, where 
k e Ker (4-) . Next, if x\ ^ 0, X2 = x 3 = 0, then Y = x\d and Z = 1, so the 
equations (8.11) give FeC and J/ £ C\u + Ker (J^) . If x\ 0, X2 ^ 0, X3 = 
0, then Y = x\d 2 + X2 and Z = d. In this case, the first equation in (8.11) 
reads 

x x F" + x 2 F = 0. 

By Lemma 8.2(b), it must be F e J 7 , and it is easy to see that the space of 
solutions is two-dimensional over C, consisting of elements of the form 



with a+ e C. Then, the second equation in (8.11) gives 

^ = a + x 12 e Xl2X - a„x 12 e~ xl2X , 
ou 

so that J/ = a + x\2^e xl2X u — a^x\2\e~ xl2X u + k, where k e Ker(Jj). 
Similarly, we consider the case x\ ^ 0, X2 = 0, x 3 =£ 0. In this case 
Y = x\d o ^jd + x 3 u' and Z = -yd. The first equation in (8.11) reads 



F'\' 

+ x 3 u'F = . 



Xi 

\ w y 

By Lemma 8.2(c), we must have F e Vo such that F' = It is easy to 
see that the space of solutions is two-dimensional over C, and it consists of 
elements of the form 

F = a+e Xl3U + a^e~ Xl3U , 
with q+ e C. Then, the second equation in (8.11) gives 

^ = a + x l3 e x ^ u - a^x 13 e' x ^ u , 
Ou 

and its solutions for J/ are of the form J/ = ct+$e Xl3U + a-\e~ Xl3U + k, for 
k € Ker (yjj). Finally, we are left to consider the case when x\ 0,X2 
0, x 3 =£ 0. In this case Y = x x d 2 a ±d + X2+ ^(»0 2 g _ x ^ u > and z = g J_g_ 
The first equation in (8.11) reads 

/F'\" F' 
x n~^) + {x 2 + x 3 {u') 2 )— -x 3 u'F = 0. 
\ u / u 

By Lemma 8.2(a), the only solution of this equation is F = 0. But then the 
second equation in (8.11) gives \f e Ker (jr)- □ 
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In the statement of Lemma 8.3 and further on in this section, we assume 
that V contains all the elements which appear in the statement, namely 

Next, for each element J/ e J-q(X), we want to find an element P e H(X) 
which is X-associated to it, and for each element P e T~Lq(X), we want to 
find an element J/ £ H(X) which is X-associated to it. Recall, by Lemma 

7.6, that if J/ < — ► P, then all elements in 1~L{X) which are X-associated to 
$/ are obtained adding to P an arbitrary element of Wq(X), and all elements 
in J~(X) which are X-associated to P are obtained adding to J/ an arbitrary 
element of Fq{X). 

Lemma 8.4. Let X = x\L\ + X2L2 + x 3 L 3 , and let 02,03,7 £ C\{0}. We 
have: 

(i) $P e H(X) such that \e^ x u P, if x 3 =£ 0; 
\e lx u <^-» ^j{x\"f 2 + X2)e~ iX , if x 3 = 0. 

(ii) $P £ H(X) such that ^ u ^ P, if x 2 ^ 0; 
^ u (x l7 2 + x 3 )e~« u u', ifx 2 = 0. 

(Hi) $u < — > (X2X + x 3 uu'). 

(iv) KA*2 + a 3 K) 2 ^ - (x.d 2 o ±d + **+**F)' d - x 3 u') V^ + asM 2 . 

(v) J0^1,i/x 2 ^0; 

^/ e ^(X) suc/i that < — > 1, i/x2 = 0, X1X3 7^ 0; 
$2^7 X ' ^/^i = 0,2; 2 = 0,x 3 j= 0; 
jf ^Ul, i/xi ^0,x 2 =0,x 3 = 0. 

(vi) J0~u', tfs 3 #0; 

$^f e T(X) such that J/ «— — » u' , if x 3 = 0, X\Xi =£ 0; 

^ u', ifxi 0,x 2 = 0,x 3 = 0. 
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Proof. The condition [e^u < — ► P means that, for some fractional decom- 
position X = YZ- 1 , there exists F e V such that P = YF and ZF = . 
Let us consider first the case x 3 7^ 0. In this case, by Lemma 8.1 a minimal 
fractional decomposition for X is (8.2) if X2 ¥= 0, and (8.4) if X2 = 0. In the 
former case Z = d o \d, hence the equation ZF = e yx reads d^r = e yx , 
namely 

F' = -e rx u" + au" , 

7 

for some a e C. This equation has no solutions since, integrating by parts, 
we get \{^e ix u" + a«"j = r yS i e' yx u, and this is not zero by (4.3). Similarly, 
in the case X2 = we have Z = —fd, hence the equation ZF = e 7X reads 

F' = e 7 V , 

which has no solutions since, integrating by parts, §e yx u' = —'y§e yx u 0. 
To conclude the proof of part (i), we consider the case X3 = 0. By Lemma 
8.1 a fractional decomposition for X is X = YZ" 1 given by (8.3). Hence, a 
solution F £ V to the equation ZF = e~ fX is F = ^e^ x , and in this case we 
have P = YF = {xid 2 )^ x = (xa + ^)e^. 

Next, let us prove part (ii). The condition ^e 7 " < — ► P is equivalent 
to the existence of F e V such that P = YF and ZF = je yu , where 
X = YZ~ l . Let us consider first the case X2 0. In this case, by Lemma 
8.1 a minimal fractional decomposition X = Y Z^ 1 for X has Z = d o 
if X3 0, and Z = d if X3 = 0. In both cases the equation ZF = je" 1 "" 1 
would imply ^e lu e dV, which is not the case by (4.3). In the case X2 = 0, 
a fractional decomposition for X is X = YZ^ 1 given by (8.4). Hence, a 
solution F e V to the equation ZF = r ye' yu is F = e 71t , and in this case we 
have P = YF = ( Xl d o ±d + x 3 u')e^ u = (xij 2 + x 3 )e 7 V. 

For part (hi) it suffices to check, using the fractional decomposition (8.2), 
that F = xv! — u e V is a solution of the equations 

YF = (x^ 2 o + ^±^Dl d _ X3U ^ F = X2X + X3UU > . 
Similarly, for part (iv), letting F = — \/ a 2 + cl 3 {u') 2 £ V, we have 
ZF = d^ = lWa2 + a 3 (ur, 

YF = — ( Xl d 2 o + x > +x u f )2 d - x 3 u>) V«2 + a 3 K) 2 • 
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Next, let us prove part (v). For x 2 ^ 0, consider the fractional decom- 
position X = YZ^ 1 given by (8.2). It is easy to check that, letting F = — , 
we have 

YF = (x x 5 2 o + **+**W 2 d - x , u ') £ = 



Hence, $0 * — ► 1, as we wanted. If x\ 0, x 2 = 0, X3 # 0, a minimal 
fractional decomposition for X is X = YZ" 1 given by (8.4). Therefore the 

relation J/ < — ► 1 is equivalent to the existence of F e V such that 

(8.12) ZF = ^ = ^, YF = (xido ±-d + x 3 u')F = 1. 

w ou V u' / 

By Lemma 8.2(c), the second equation in (8.12) implies that F e Vo is such 
that F' = §^u'. In this case, the second equation in (8.12) reads 

dF 

x\d— + x 3 Fu' = 1 , 
ou 

which, by the commutation relation (4.1), is equivalent to 

/ d^F \ 

But obviously, the above equation is never satisfied. If x\ = 0, X2 = 0, x% j= 

0, it is easy to check that F = —1 solves 

' j X311' 



X3U 1 v! x 3 u' 5u 2x 3 u' ' 

proving that \ 2x311' ^ Finally, for x\ =£ 0,X2 = 0,2:3 = 0, a minimal 

fractional decomposition X = Y Z~ l is given by Y = x\d and Z = 1. In 
this case, it is immediate to check that F = — solves 

X\ 

„ x x 5 n xu 

YF = Xl d— = 1 , ZF = — = —\— , 

X\ X\ OU X\ 

proving that Jff * — * 1- 

We are left to prove part (vi). For x 3 ^= 0, consider the fractional decom- 
position X = YZ^ 1 given by (8.2). It is easy to check that, letting F = 
we have 

zf = do\d^- = 0, 



YF 



(x^o^ d+ ^Dl d - X3U ')^ = u'. 
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Hence, $0 < — > u', as we wanted. If X\ 0, X2 0, x$ = 0, the minimal 
fractional decomposition for X is X = Y Z~ l given by (8.3). Therefore the 

relation J/ < — ► u' is equivalent to the existence of F e V such that 
(8.13) ZF = F' = -/- , YF = ( Xl d 2 + x 2 )F = 



Su 

If F e V n , for n ^ 0, we get, applying ^ ( f +2) to both sides of the second 



" 7 + Xl)F = u 

to both sidi 

equation in (8.13), that g ^ n) = 0. Hence, it must be F e J 7 . But in this 
case, the second equation in (8.13) has clearly no solutions. If x\ = 0,X2 
0, X3 = 0, it is easy to check that F = ^ solves 

if = rE 2 — = u , zF = d — = —)- 



' X2 X2 5u 2X2 



proving that § ^ «^-» u'. Finally, if X\ =£ 0,X2 = 0,2:3 = 0, a minimal 
fractional decomposition X = Y Z~ l is given by Y = x\d and Z = 1. In 
this case, it is immediate to check that F = ^- solves 

^ ^ - u , „„ It 5 r v? 
YF = x 1 d—=u', ZF= — = —\—, 

X\ X\ OU iX\ 

proving that j^- < — ► u'. □ 

In order the check orthogonality conditions (7.16) for Liouville type in- 
tegrable systems we will use the following results. 

Lemma 8.5. (a) (CI) 1 = Im(d). 

(b) (Cu^ = lm(±d). 

(c) (Span c {l,^}) 1 = Im(5oJ 77 5). 

(d) For 62^3 e C\{0}, we have 

Proof. Parts (a) and (b) are immediate. Let us prove part (c). It is im- 
mediate to check, integrating by parts, that $(a + f3u')d^n = for every 
a, (3 e C. Hence, lm(d a \d) cz ( Span c {l, n'})" 1 . On the other hand, if 
/ E ( Span c {l, n'}) 1 , it must be 

f = dg = — dh for some g,h e V . 
w 
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But then dh = u'dg = d(u'g) — u"g, which implies g = ^Trd(u'g — h). Hence, 

f = d^d(u'g-h)elm(dol J d). 
u" u" 

We are left to prove part (d). We have 



2 



hSu J V & 2 + b 3 {u') 

The inclusion Im ( b2 +W 2 g _ j^A c ( Qd^^L \*~ follows by inte- 

gration by parts, and the following straightforward identity 

150 » + hu )d , ; == = 0. 



We are left to prove the opposite inclusion. If f e (Cd , = 

V V 6 2+63(«' 

have 



_L 

, we 



. = dg (b 2 + b 3 (u') 2 )2 

1 8 , u ' b 2 u" ° 9 ■ 

for some be V. Letting g = , . we then get: 

-y/62 +6 3 («') 

(&2 + 6 3 K) 2 )i /» 

■J L _ . // 



& 2 u" V^ + bsK) 2 

(6 2 + 6 3 (n') 2 )i / 



-/? 



M" V (b 2 + b 3 (u') 2 )l V&2+W) 



=h>) 

A2 ) 



-d-hu)-. 



□ 



8.2 Integrability of the Lenard-Magri scheme: 61 = 

In this and the next two Sections we consider the case when b\ = 0, for 
which we get integrable Lenard-Magri schemes of S-type, in the terminology 
introduced in Section 7.5, in the case a± (described in Section 8.3 below), 
and of C-type in the case a\ = (described in Section 8.4 below). In Sections 
8.5 we will consider the remaining case, when b\ 0, for which we again get 
some integrable Lenard-Magri schemes of C-type. 
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According to Theorem 7.15 (and Remark 7.10), in order to apply suc- 
cessfully the Lenard-Magri scheme of integrability, we need to find finite 
sequences {P n }^ =0 , {^n}^Lo> satisfying the relations (7.12) for £ n = 
or equivalently the relations (7.23), and the orthogonality conditions (7.16). 
For b\ = we display below such sequences separately in all possibilities for 
the coefficients a 2 ,a 3 ,b 2 ,b 3 being zero or non-zero, and a± arbitrary. (Note 
that, since we are assuming b\ = 0, we don't need to consider the case 
6 2 = 63 = 0.) 

, , , ^ , ^ H -, K r ^ H , K 



[111 
(iv 
(v 
(vi 
(vii 
(viii 
(be 

(x 
(xi 
(xii 



6263 ± 0,a 2 a 3 ± 0: $0 A 1 A $0 A v! A W b 2 + b 3 (u') 2 . 
b 2 b 3 ± 0,a 2 # 0,a 3 = 0: $0 ^ 1 A \^b 2 + b 3 {u'f. 
b 2 b 3 # 0, a 2 = 0, a 3 7^ 0: $0 ^ «' ^> W b 2 + b 3 {u') 2 . 
b 2 b 3 ^0,a 2 = a 3 = 0: ^AoA Wb 2 + b 3 (u') 2 . 
b 2 * 0,63 = 0,a 2 a 3 ^ 0: JO ^ 1 ^ JO ^ n' ^ 
62 ^ 0, 63 = 0, a 2 7^ 0, a 3 = 0: $0 ^ 1 ^ JO. 
6 2 ^0,6 3 = 0,a 2 =0,a 3 7^0: $0 ^ u' ~ 
b 2 ^ 0, 63 = 0, a 2 = a 3 = 0: $0 ^ A Jo. 
62 = 0,63 * 0,a 2 a 3 7^ 0: $0 ^ u' $0 A 1 ♦ J^. 
62 = 0,6 3 ^0,a2^0,a 3 = 0: loii^Jj^. 
62 = 0, 63 7^ 0, a 2 = 0, a 3 7^ 0: $0 u' ^> JO. 
62 = 0, 63 7^ 0, a 2 = a 3 = 0: $0 ^ A Jo. 



All the above iT- and i^-association relations hold due to Lemmas 8.3 and 
8.4. Moreover, using Lemmas 8.1 and 8.5 we check that both orthogonality 
conditions (7.16) hold. Hence, by Theorem 7.15 and Remark 7.10 all the 
above sequences can be continued indefinitely, possibly going to a normal 
extension V of V, to an infinite sequence 

50 A p A fa A Pl A $ hl . . . . 
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Note that, by Lemma 7.6, at each step the subsequent term is unique up to 
a linear combinations of the previous steps. 

Next, we want to discuss integrability of the corresponding hierarchies 
of Hamiltonian equations 4r- = P n , n e Z + . Namely, according to Defini- 
tion 7.9, we need to see when the vector spaces Span c {$/i n } c V/dV and 
Span c {P n } cz V are infinite dimensional. 

First, we consider the cases (vi), (viii), (xi) and (xii), where we show 
that integrability does not occur (regardless of a\ being zero or non-zero) 
since the Lenard-Magri scheme repeats itself. In case (vi), by Lemmas 8.3 

and 8.4 we have 7io(H) = C and, for every a e C, {$/ e F(K) \ $/ < — * a} = 
Fo{K) = Ker(^). Hence, any infinite sequence extending the given finite 
one will have §h n e Ker (^) and P n e C, for every n £ Z+. Similarly, in case 

(xi) we have Hq(H) = Cu' and, for every au' e Cv! , {J/ e T{K) \ J/ 
au'} = J-q(K) = Ker (-£-) . Hence, any infinite sequence extending the given 
finite one will have \h n e Ker (^) and P n e Cu', for every n e Z + . In 
cases (viii) and (xii) we have T-Lq(H) = and Fq{K) = Ker(^). Hence, 
^h n e Ker and P n = 0, for every n e Z + . In conclusion, in all these 
cases Span c {P n } is finite dimensional, and integrability does not occur. 

For the remaining 8 cases, we will prove in Section 8.3 that when a\ =£ 
we get some integrable Lenard-Magri scheme of S-type, and we will prove in 
Section 8.4 that when cq = we get some integrable Lenard-Magri scheme 
of C-type (in the terminology of Section 7.5). 

8.3 Integrable Lenard-Magri schemes of S-type in the case 
b\ = and a\ ^ 

Cases (i), (u), (Hi), (iv) 

In all the sequences (i)-(iv), after one or two steps, we arrive at (after shifting 
indices) J/i-i = 5V^2 + bs(u') 2 . The next term in the sequence, which we 
denote Po, is obtained by solving the following equations for F and Po i n V: 

BF = dol^dF = JL$yb2 + hW, 
u" 5u JJ 

Po = AF = (a lS * o \ id + Q2 + Q ;^ 8 - a 3 u>) F . 

V u u" / 
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It is easy to check that a solution is given by F = —yb? + b 3 (u') 2 , and 

aib?,u' \" (a 3 b 2 - a 2 b 3 ) 



Pn 



3.14) 



a\b 2 b 3 v!" < a^bju' \u") 2 < (a 3 b 2 - a 2 b 3 



(62 + b 3 {u') 2 )2 (62 + b 3 {u') 2 f2 V & 2 + 6 3 («') 2 

The above computation works regardless whether a% is zero or not. But 
to prove that the Lenard-Magri scheme is integrable of S-type we need to 
assume a\ =fc 0, in which case dord(Pn) = 3 is greater than max{dord(>l) — 
\H\ + \K\, dord(P) + \K\, dord(C), dord(D) + \K\}, which is less than or 
equal to 2 (for all the cases (i)-(iv)). Therefore, by Corollary 7.16, each 
Hamiltonian PDE ^ = P n , n e Z+, is integrable, associated to an integrable 
Lenard-Magri scheme of S-type. (Note that, since Ker(P*) n Ker(L>*) = 
C ®Cv! j= 0, we cannot conclude that [P m ,P n ] = for every m,n e Z + , 
and therefore that the sequence of equations 4^- = P n , n e Z + , form a 
compatible hierarchy.) 

After rescaling x and t appropriately in the equation % = Pn, we con- 
clude that the following bi-Hamiltonian equation is integrable, associated to 
an integrable Lenard-Magri scheme of S-type: 

, . du u'" u'(u") 2 a 
8.15 — = r -3 V 5 + r, aeC. 

dt (l + (u') 2 )t (l + (-u') 2 )f {l + {u') 2 )2 

This is an equation of the form [MSS90, eq.(41.5)] with 03 = (1 + (u') 2 )5. 
This particular 03 does not appear in their list, but as V. Sokolov informed 
us, a simple point transformation reduces (8.15) to an equation from their 
list. 

Cases (v), (vii) 

In the sequences (v) and (vii), after one or two steps, we arrive at (after 
shifting indices) J/t-i = \ ^ ■ The next term in the sequence, which we 
denote Pn, is obtained by solving the following equations for F and Pq hi V: 

u" 5u JJ 2b 2 

Po = AF= UP o L k B + * 2 + -f) 2 d _ aW ) p . 

It is easy to check that a solution is given by F = jjjg, and 
(8.16) P = -u + _„+_(„) . 
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As before, if a% 0, then dord(Po) = 3 is greater than max{dord(^4) — 
\H\ + \K\, dord(P) + \K\, dord(C), dord(L>) + \K\}, which is at most 2. 
Therefore, by Corollary 7.16 each Hamiltonian PDE ^ = P n , n e Z+, is 
integrable, associated to an integrable Lenard-Magri scheme of S-type. 

After rescaling x and t appropriately in the equation ^ = Pq, we con- 
clude that the following bi-Hamiltonian equation is integrable, associated to 
an integrable Lenard-Magri scheme of S-type: 

(8.17) ^ = u'" + ev! + a{v!f , 

where e is 1 (in case (v)) or (in case (vii)) and a e C. By a Galilean 
transformation we can make e = 0. The resulting equation is called the 
potential modified KdV equation (equation (4.11) in the list of [MSS90]). 

Cases (ix), (x) 

In all the sequences (ix) and (x), after one or two steps, we arrive at (after 
shifting indices) = \ 2b3u' • ^ ne nex ^ term in the sequence, which we 

denote Pq, is obtained by solving the following equations for F and Pq in V: 

BF = do\dF=^-\^—, 
u" 5uhb 3 u' 

o at? ( p2 1 a 2 + a 3 (u') 2 A 
P = AF = [ aid o —d H d - a 3 u )F . 

It is easy to check that a solution is given by F = 2 ^ u , , and 
ai u'" 3ai (u") 2 a 2 1 a 3 



5-18) P = - 



If a\ ^ 0, we have dord(Po) = 3, which is greater than max{dord(^4) — 
\H\ + \K\, dord(P) + \K\, dord(C), dord(D) + \K\}, which is at most 2. 
Therefore, by Corollary 7.16, each Hamiltonian PDE ^~ = P n , n e Z + , is 
integrable, associated to an integrable Lenard-Magri scheme of S-type. 

After rescaling x and t appropriately in the equation 4| = Pq, we con- 
clude that the following bi-Hamiltonian equation is integrable, associated to 
an integrable Lenard-Magri scheme of S-type: 

, x du v!" {u") 2 1 

(8 - 19) ^ = RF~ 3 W + W + a ' aeC - 

As explained in [MSS90], by a point transformation one can reduce this 
equation to an equation of the form (4.1.4) in their list. 



92 



Remark 8.6. Note that equation ^ = Pq with Pq given by (8.18) is trans- 
formed, by the hodograph transformation x >—> u, u >— * —x, to the equation 
with Pq given by (8.16), after exchanging a 2 and 62 with 03 and 63 respec- 
tively. Equivalently, equation (8.19) can be transformed to equation (8.17) 
up to rescaling of x and t. 



8.4 Integrable Lenard-Magri schemes of C-type with a± = 

61 = 

Cases (i), (ii), (Hi), (iv) 

As pointed out above, in all the sequences (i)-(iv) we arrive, after one or 
two steps (and after shifting indices), at J/t-i = + b 3 (u') 2 . In the case 

a\ = we can actually find an explicit solution for the sequences {\h n } ne z + 
and {P n }neZ + satisfying the recursive formulas 

(8.20) P n \h n P n+1 , n e Z + . 

It is given by (n ^ 0): 



_ A fn\ (2n - 1 - 2fc)H A n+1 - fe 4 

(8.21) 



-u' 



h 



y kJ (2n-2k)U b% (bi + b^u'yfi+n-^ 
n / n\ (2n - 1 - 2fc)!! A™* 1 "^ -u' 



, (2n-2fc + 2)!! ^ (&2 + 6 3 (^)2)|+— ^ ' 

where A = 0263 — 03621 which is non-zero unless the operators H and if are 
proportional. Here and further we let (—1)!! = 1. 

First, note that for n = 0, the above expression for Pq is the same as the 
one in (8.14) with a\ = 0. Hence \h—i <-^-> Pq. Next, we check that indeed 
the sequences {\h n } ne z+ and {P n }neZ + solve the recursive relations (8.20). 
For this, we fix the fractional decompositions H = AB' 1 and K = CD^ 1 
given by 

02 + 03(^0% , n b 2 + b 3 {u')\ , 1 

A = d-a 3 u , C= d-b 3 u , B = D = d o — d . 

u" u" u" 

Since B = D, the relations (8.20) hold if there exists an element F n such 
that 

Sh 

(8.22) CF n = P n , BF n = — — , P n +i = AF n . 

ou 
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A solution F n to equations (8.22) is given by F n = —h n . Since h n depends 
only on v! , we have 

B(- M = - a „- aft „ = (- a) ^ = ^, 

hence F n = —h n satisfies the second equation in (8.22). The first and third 
equations in (8.22) can be easily checked using the following straightforward 
identities (m e Z+): 

-b 3 mu' 



ib 2 + b^u'Y)T (b 2 + b 3 (n')^ ' 
l . —v! 

A— — ; ; = mA- — ; / , snsl n +1 + {m + l)a 3 - 



(b 2 + b 3 (u')i)- {b 2 + b 3 {u'Y)- +1 {b 2 + b 3 {u'Y)- 

If A # 0, all the elements P n are linearly independent, therefore, by 
Theorem 7.15 and Lemma 7.20, each Hamiltonian PDE % = P n , n e Z+, is 
integrable, associated to an integrable Lenard-Magri scheme of C-type. 

Cases (v), (vii) 

In the sequences (v) and (vii) we arrive, after one or two steps (and after 
shifting indices), at J/i-i = $ ^ . In the case a\ = we can actually find 
an explicit solution for the sequences {^h n } ne z + and {P n }neZ + satisfying the 
recursive formulas (8.20). It is given by (n e Z+): 



A /n\ (2fc-l)!!^- fc a§ 2fc+1 

> oq^l fc=0 V 7 V ; 2 

' , v M(2*-i)l!«rM, ., 



First, note that the above expression for Pq is the same as the one in 
(8.16) with ai = 0. Hence J/i-i *— * i~b- Next, we check that indeed the 
sequences {\h n } ns i + and {P n }neZ + solve the recursive relations (8.20). For 
this, we fix the fractional decompositions H = AB^ 1 and K = CD^ 1 given 
by 

A = ° 2 + a l {u ' )2 d - a 3 u> , B = do±-d, C = b 2 , D = d. 
u u 

The relations (8.20) hold if there exist elements F n ,G n e V such that 

Sh Sh 
(8.24) CF n = P n , DF n = , BG n = — ^ , P n+ i = AG n . 
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Solutions F n , G n to equations (8.24) are given by F n = j^P n and G n = —h n . 
The first and third equations in (8.24) are immediate. The third equation 
follows from the immediate identity y 1 = — • Finally, the fourth identity 
in (8.24) is easily checked using the Tartaglia-Pascal triangle. 

Clearly, if 03 ^ 0, all the elements P n are linearly independent, therefore, 
by Theorem 7.15 and Lemma 7.20, each Hamiltonian PDE ^~ = P n , n e Z + , 
is integrable, associated to an integrable Lenard-Magri scheme of C-type. 



Cases (ix), (x) 

In the sequences (ix) and (x) we arrive, after one or two steps (and shift of 
indices), at = • ^ ai = ^ we can ^ nc ^ an ex P nc it solution for the 

sequences {\h n } ne i + and {P n }neZ + satisfying the recursive formulas (8.20). 
It is given by (n e Z + ): 

= y (n\ (2fe-l)!!ag- fc q§ 1 

18 bJ A (2fc - 1)!! ag- fc q* 1_ 

~ U / (2A; + 2)!! (w') 2fc+1 ' 

First, note that the above expression for Pq is the same as the one in 
(8.18) with ai = 0. Hence Jo_i <-^-> P - Next, we check that indeed the 
sequences {J/in}neZ + and {P n }neZ+ solve the recursive relations (8.20). For 
this, we fix the fractional decompositions H = AB~ l and K = CD^ 1 given 
by 

a 2 + a 3 (-u') 2 , 1 , 1 

^4 = o — a-su , B = oo—d, 6 = 030 , D = —0 . 

u" u" w 

Since equations (8.24) hold with F n = j^jP n and G n = —h n (a fact that 
can be easily checked directly), it follows that the recursive relations (8.20) 
hold. 

Clearly, if 02 =£ 0, all the elements P n are linearly independent, therefore, 
by Theorem 7.15 and Lemma 7.20, each Hamiltonian PDE ^ = P n , n e Z+, 
is integrable, associated to an integrable Lenard-Magri scheme of C-type. 



8.5 Integrable Lenard-Magri scheme of C-type with b\ ^ 

As we did in the previous sections, we study here the integrability of the 
Lenard-Magri scheme when b\ =£ 0. We will consider separately the various 
cases, depending on the parameters 62,63,02,03 being zero or non-zero. 
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Case 1: b 2 b 3 =£ 



Let us consider first the case when b 2 and 63 are both non-zero. If $/io e V/dV 

and Pq e V satisfy the relations $0 * — * Pq < — » $/io, then, by Lemmas 
8.3 and 8.4, we necessarily have Pq e C © Cii' and ^ = 0. Hence, any 
infinite sequence extending the given finite one will have \h n e Ker (4r;) and 
P n e C@Cu', for every n £ In other words, the Lenard-Magri scheme 
repeats itself and integrability does not occur. 



Case 2: b 2 ¥= 0, 63 = 0, a 3 = 

In the case when b\b 2 0, 63 = and a 3 = 0, we can find explicitly all 
possible solutions for the sequences {S j h n } ne z + and {P n }neZ + satisfying the 
Lenard-Magri recursive relations (7.23). 

In order to describe such solutions, we need to introduce some polynomi- 
als. We let p n (x; A, e), q n (x; A,e), n e Z+, be the sequences of polynomials, 

depending on the 2x2 matrix A = \ °^ f 2 ) , and on the sequence of 



b\ b 2 

constant parameters e = (eo,ei, ■ ■ ■ ), defined by the following recursive re- 
lations: po(x; A, e) = 0, and 



Pn+\{x; A, e) = ^-p n (x; A, e) + Ql&2 g n ( x; A, e) 



(8.26) d2 

(-To + 2bi 2 — )q n (x;A,e) = p n (x;A,e) . 
V dx A ax / 

Here and further, as before, we use the notation (8.9) with X{ and Xj replaced 
by bi and bj. It is easy to see that, if p(x) is a polynomial of degree n, then 
a solution q(x) of the differential equation q"(x) + 2bi 2 q'(x) = p(x) is a 
polynomial of degree n + 1, defined uniquely up to an additive constant 
eo- Hence, at each step in the recursion (8.26), the resulting polynomial 
Pn+i(x) depends on the previous step p n (x) and on the choice of a constant 
parameter e n+ \. 

With the above notation, all sequences {\h n } ne i + , {P n }neZ + , satisfying 
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the Lenard-Magri recursive relations (7.23), are as follows: 

P n = Pn (x; A, e+)e bl2X + Pn(~x; A, e") e - b ^ x + a 2 5 n , 

(8.27) h n = 1 (^(x; A, e+) + 6 12<?n (x; A, e + )) e bl2a: u 

- ^-(^(-^AO + b 12 q n {-x; A,e~))e~ bl2X u, 

where e- = (eg", ef, . . . ) and (5 = (<5o, <Si, . . . ) are arbitrary sequences of 
constant parameters. 

It is not hard to check that the sequences {$/i n }neZ + and {P n }neZ + indeed 
solve the recursive relations (7.23), and any solution of the recursive relations 
(7.23) is obtained in this way. To conclude, we observe that, since A = 
a 2 b\ — «i&2 is non-zero (unless the operators H and K are proportional), 
all the elements P n are linearly independent, therefore, by Theorem 7.15 
and Lemma 7.20, each Hamiltonian PDE % = P n , n e Z+, is integrable, 
associated to an integrable Lenard-Magri scheme of C-type. 

Case 3: b 2 ¥= 0, b 3 = 0, a 3 # 

In the case when b\b 2 ¥= 0, 63 = and 03 ^ 0, we have by Lemma 8.3(b) 
that %q{H) = a 2 C®Cu'. On the other hand, by Lemma 8.4(vi) there is 

no element J/ e J~(K) such that J/ <_ ^' it'. Similarly, by Lemma 8.3(a) 
we have T (K) = C\e bl2X u + C\e~ bl2X u + Ker (^). On the other hand, by 

Lemma 8.4(i) there is no element P e F{H) such that \e ±bl2X u <^> P. 

In conclusion, the Lenard-Magri recursion scheme, in this case, cannot 
be applied, since the following finite sequences cannot be extended to infinite 
sequences satisfying (7.23): 

JO u' tlf , JO JU \e ±b ™ x u $P . 
Case 4: b 2 = 0, 63 ^ 0, a 2 = 

In the case when 6163 ^ 0, 62 = and a2 = 0, we can find explicitly all 
possible solutions for the sequences {$/i n }neZ+ and {P n }neZ + satisfying the 
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Lenard-Magri recursive relations (7.23). They are as follows: 



P n = p n {u- A, e + )e fel3 " + Pn(-u; A, e -)e~ b ^ u + a 3 6 n u' , 
^ = A,e+) + b 13 q n (u; A,e + ))e 6l3U 



- i (<£(-«; A, e") + h 3 q n (-u; A, e")) e" 6 "" , 

where p n (-u;A,e + ) and g n (it;A, e + ) are the polynomials defined in (8.26), 

depending on the matrix A = ( °. 1 ? 3 ) , and on the sequences of constant 

V oi 03 / 

parameters e- = (eg- , ej- , . . . ) and 5 = (5q,5i, . . .). 

It is not hard to check, as in case 1, that the sequences {\h n } ne j J+ and 
{P n }neZ+ solve the recursive relations (7.23), and any solution of the recur- 
sive relations (7.23) is obtained in this way. We also observe that, since 
A = 0361 — 0163 is non-zero (unless the operators H and K are propor- 
tional), all the elements P n are linearly independent, therefore, by Theorem 
7.15 and Lemma 7.20, each Hamiltonian PDE % = P n , n e Z + , is inte- 
grable, associated to an integrable Lenard-Magri scheme of C-type. 

Case 5: b 2 = 0, 63 ^ 0, a 2 ¥= 

In the case when 6163 ^ 0, 62 = and 03 = 0, we have by Lemma 8.3(b) 
that T-Lq(H) = C © a 3 Cu' . On the other hand, by Lemma 8.4(v) there is 
no element §f e J~(K) such that J/ 1- Similarly, by Lemma 8.3(a) we 
have Tq{K) = CS i e bl3U + C\e- bwu + Ker (^) . On the other hand, by Lemma 

8.4(ii) there is no element P e J~(H) such that ^e- bl3X u +—> P. 

In conclusion, the Lenard-Magri recursion scheme, in this case, cannot 
be applied, since the following finite sequences cannot be extended to infinite 
sequences satisfying the relations (7.23): 

Case 6: b 2 = 63 = 

In the case when b\ =£ 0, which we set equal to 1, and b 2 = 63 = 0, we have 
different possibilities according to the constants a 2 and 03 being zero or not. 

If 0203 7^ 0, the Lenard-Magri recursion scheme cannot be applied. In- 
deed, by Lemma 8.3 we have Hq(H) = C®Cu' and F Q (K) = Cjit ©Ker (^), 
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§XU < 
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SO- 


H 


< 


K 


f 1 



and, whichever way we start the finite sequences {S^n}^=o; {-Pnj^Lo as i n 
(7.23), there is no way to extend them to non-trivial infinite sequences: 



—* ci2X + a^uu' < — > $^h\ . 

Next, we consider the cases when exactly one of the elements <Z2 and 
ci3 = is zero. Recall the sequence of polynomials p n (x; A, e) defined by the 
recursive equations (8.26). In the case 61 = 1,62 = 0, such equations reduce 
to po(x; ai, a,2, e) = and 

fOLl 0,2 ( d \ ~ 2 \ 

(8.29) p n+1 (x;ai,a 2 ,e) = [— + —[ — ) )p n {x;ai,a 2 ,e) 



b\ b\ \dx 

Here (^^J means integrating twice with respect to x, which is defined 
uniquely up to adding a linear term €2 n + ^2n+\x. In particular, at each step 
the degree increases by two. 

In the case a 2 0, (13 = 0, it is not hard to prove that all the sequences 
(S^n}neZ + , {Pn}neZ + , satisfying the Lenard-Magri recursive relations (7.23), 
are as follows: 

(8.30) P n = p n (x;ai,a2,e) + 8 n , J/i n = S(^J Pn( x ; ai, a 2 , e)u , 

where e = (eo,ei, ...) and (6q,6i, . . .) are arbitrary sequence of constant 
parameters. Since, obviously, all the elements P n are linearly independent, 
we conclude that each Hamiltonian PDE ^ = P n , n e Z+, is integrable, 
associated to an integrable Lenard-Magri scheme of C-type. 

Similarly, in the case 02 = 0,03 # 0, all the sequences {S l h n } ne i + , 
{P n }neZ + , satisfying the Lenard-Magri recursive relations (7.23), are as fol- 
lows: 

/ d \-2 

(8.31) P n = p n (u;ai,a 2 ,e)u' + 5 n u' , \h n = lyj-j Pn(u; at, a 2 , e) , 

where e = (eo,ei,...) and {5q,5i, . . .) are arbitrary sequences of constant 
parameters. Again, we conclude that each Hamiltonian PDE ^ = P n , n e 
Z+, is integrable, associated to an integrable Lenard-Magri scheme of C- 
type. 
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8.6 Summary 

Let us summarize the results from the previous sections by listing all the 
possibilities for the pairs H and K as in (8.1), and specifying, using the ter- 
minology of Section 7.5, whether the corresponding Lenard-Magri sequence 
(7.23) is integrable of S-type, i.e. the orders of the elements P n 's and T^f's 
tend to infinity, {b\ = 0, a\ # 0), whether it is integrable of Ci-type, i.e. the 
C-span of the elements P n 's and Jn n 's is infinite dimensional and the orders 
of H and K are both equal to —1 (b\ = a\ = 0), whether it is integrable 
of Ci-type, i.e. the C-span of the elements P n 's and §h n 's is infinite dimen- 
sional and H has order less than or equal to K and K has order 1 (pi =£ 0), 
whether it is of finite type, i.e. the C-span of the elements P n 's or $o n 's is 
necessarily finite dimensional, or whether it is blocked, i.e. there are choices 
of \h n or P n for which the scheme cannot be continued. This is the list of 
all possibilities: 

• integrable of S-type: 

(a) bi = 0, (62,63) # (0,0), a ia2 a 3 ^ 0; 

(b) b\ = 0, a± ^ 0, and either 62 ¥= 0, 02 = and (63,03) # (0,0), or 
h ± 0, a 3 = and (6 2 ,a 2 ) 7^ (0,0). 

• integrable of C\ -type: 

b\ = a>\ = 0, and either 6203 and (63,02) arbitrary, or 6302 ¥= 
and (62,03) arbitrary. 

• integrable of Ci-type: 

(a) b\a\ ^ 0, and either 62 = a 2 = and (63,03) (0,0), or 63 = 
a 3 = and (62,02) (0,0); 

(b) 61 7^ 0, oi = 0, and either 62 = 02 = and 03 (with 63 
arbitrary), or 63 = 03 = and 02 7^ (with 62 arbitrary). 

• finite type: 

(a) 6^263 ^ 0, oi = 0, (02,03) 7^ (0,0); 

(b) 61 = 0, ai 7^ 0, and either 62 = 02 = and 63 7^ (with 03 
arbitrary), or 63 = 03 = and 62 7^ (with 02 arbitrary). 

(cl) b\ = a\ = 0, and either 62 = 02 = and 6303 ^ 0, or 6202 7^ 
and 63 = 03 = 0; 

(c2) 616263 7^ 0, 010203 7^ 0; 
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(d) 6i6 2 ^3 ^ 0, a\ 0, a 2 a 3 = 0. 



• blocked: 

(a) b\ =£ 0, a\ = 0, and either b 2 = 0, a 2 =£ and (63,03) ^ (0,0), or 
h = 0, a 3 =£ and (6 2 ,a 2 ) 7^ (0,0); 

(b) 61 ^ 0, 6 2 6 3 = 0, aia 2 a 3 0; 

(c) b\a\ ^ 0, and either b 2 a 3 ^ and b 3 = a 2 = 0, or b 2 = 03 = 
and 6302 ¥= 0. 

8.7 Going to the left 

Suppose we have an integrable Lenard-Magri sequence (7.23) (of S, Ci or 
Qrtype). A natural question is whether this sequence can be continued to 
the left: 
(8.32) 

In this way we get some additional equations compatible with the given 
hierarchy Jr- = P n , n e Z+, and additional integrals of motion §h n , n = 
— 1, —2, . . . , in involution with the given $/i n 's, with n ^ 0. 

Clearly, trying to extend the Lenard-Magri scheme (8.32) to the left 
amounts to switching the roles of the non-local Poisson structures H and 
K, and to constructing the "dual" Lenard-Magri sequence 

( 8 .33) JO P_! \h-x P-2 ~ S h -2 ^ P-3 " ■ ■ ■ ■ 



So, we need to study, for each possible choice of the parameters aj,6j,i = 
1,2,3, what type of Lenard-Magri scheme we get when we switch all the 
coefficients o^'s with the b^s. 

By looking at the list of all possibilities in the previous section, after 
switching the roles of H and K we have the following: 

• integrable of S-type (a) ^— * finite type (a); 

• integrable of S-type (b) * blocked (a); 

• integrable of Ci-type ^— ? integrable of Ci-type; 

• integrable of C2-type (a) * integrable of C2-type (a); 

• integrable of C2-type (b) ^— * finite- type (b); 
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• finite-type (c) 



finite- type (c); 



• finite-type (d) 



blocked (b); 



• blocked (c) 



blocked (c). 



We are only interested in the integrable (S or C-type) Lenard-Magri 
schemes. We see from the above list that, after exchanging the roles of H 
and K, three things can happen. The "dual" Lenard-Magri scheme (8.33) 
can be of finite-type (this happens in the cases S(a) and 02(b)). In this 
situation continuing the Lenard-Magri scheme to the left we never get any 
new interesting integrals of motion or equations. 

The second possibility is that the "dual" Lenard-Magri scheme (8.33) is 
of integrable-type (this happens in the cases Ci and 02(a)). In this situation 
we can continue the Lenard-Magri scheme to the left indefinitely. In other 
words, in each of these cases we can merge two integrable systems, "dual" 
to each other, to get one integrable system with twice as many integrals of 
motion and equations. 

The most interesting situation is when the "dual" Lenard-Magri scheme 
(8.33) is blocked (which happens in the case S(b)). In this case, if the 
sequence (8.33) is blocked at P-k, k > 0, we obtain an integrable PDE 
which is not of evolutionary type. 

8.8 Non-evolutionary integrable equations 

According to the previous discussion, we need to consider the case of in- 
tegrable Lenard-Magri scheme of S-type (b), which means the following 5 



1. 


h 


= 0,62 


^0,63 


^ 0,oi ¥= 0,a 2 + 0,a 3 = 0; 


2. 


h 


= 0,62 


^0,63 


^ 0,ai * 0,a 2 = 0,a 3 # 0; 


3. 


h 


= 0,62 


^0,63 


# 0, ai # 0, a 2 = 0, a 3 = 0; 


4. 


h 


= 0,62 


^0,63 


= 0, ai # 0,a 2 = 0,a 3 # 0; 


5. 


h 


= 0,62 


= 0,63 


^ 0,ai # 0,a 2 ^ 0,a 3 = 0; 



cases: 



102 



Case 1: b\ = 0, b 2 ¥= 0, b 3 =£ 0, a x 0, a 2 =t 0, a 3 = 

This case gives, to the right, the Lenard-Magri scheme listed as case (ii) in 
Section 8.3, while, after exchanging the roles of H and K we get, to the left, 
the "blocked" Lenard-Magri scheme listed as case 3 in Section 8.5. Hence, 
overall, we get the following scheme: 

$P ^ le ±a ^ x u JUoJ^loJUlJ^ W b 2 + h(u') 2 ^ 

( g > aib 2 b 3 u"' ^ aib 2 blu'(u") 2 a 2 b 3 jc^ 

^ (b 2 + b 3 (u')i)l [b 2 + b 3 {u'Y)i y/h + h(u') 2 " " 

Here and further + means that we take arbitrary linear combination of the 
above expressions with + and with — . Trying to solve naively for P in the 
above scheme, we get the following expression 

P = ±h- e ±ai2X + b 3 u'd- l (e ±ai2X u') . 

The meaning of the above expression for P is that the following partial 
differential equation is a member of the integrable hierarchy associated to 
the Lenard-Magri scheme of S-type (ii): 

(8.34) (—) = +—( — e ±ai2X ) +b 3 e ±ai2X u x . 

V u x ) x a\ 2 \u x ) x 

Case 2: &i = 0, b 2 ^ 0, b 3 0, a\ ^ 0, a 2 = 0, a 3 # 

This case gives, to the right, the Lenard-Magri scheme listed as case (hi) in 
Section 8.3, while, after exchanging the roles of H and K we get, to the left, 
the "blocked" Lenard-Magri scheme listed as case 5 in Section 8.5. Hence, 
overall, we get the following scheme: 



$P ^ \e ±a ^ u JU^JU^JUu'^ Wb 2 + b 3 {u>y J± 

( g > aib 2 b 3 u"' 3 aib 2 blu'(u") 2 a 3 b 2 ^ k 

^ (b 2 +b 3 (u') 2 )l {b 2 + b 3 {u'Y)i V^ + WF" 

Trying to solve naively for P we get 

P = +a 13 b 2 d~ 1 e ±ai3U + b 3 u'e ±ai3U . 
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This means that the following hyperbolic partial differential equation is a 
member of the integrable hierarchy associated to the Lenard-Magri scheme 
of S-type (iii): 

(8.35) u tx = ±a l3 b 2 e ±a ^ u ± (e ±ai3U ) xx . 

013 xx 

Case 3: 61 = 0, b 2 ^ 0, 63 # 0, a\ ^ 0, a 2 = 0, a 3 = 

This case gives, to the right, the Lenard-Magri scheme listed as case (iv) in 
Section 8.3, while, after exchanging the roles of H and K we get, to the left, 
the "blocked" Lenard-Magri scheme listed as case 6 in Section 8.5. Hence, 
overall, we get, depending on how we choose to continue the scheme to the 
left, the following two possibilities (or any their linear combination): 

^ P !j M JU 1 50 JL JL Wb 2 + ~ 

ai 

H aib 2 b 3 u'" aib 2 blu' {u") 2 K 



or 



{b 2 + b 3 {u'Y)l (b 2 + b 3 (u') 2 Y- 



H aib 2 b 3 u"' aib 2 blu' {u") 2 k 



(b 2 + b 3 (u>) 2 )2 (b 2 + b 3 (u>) 2 )2 
Trying to solve naively for P we get, in the first case 
d b 2 2 b 3 , b 3 , x 

P = X -\ XUU U V u , 

2a\ al a\ 

which corresponds to the following integrable non-evolutionary partial dif- 
ferential equation: 

(8.36) m =M-) 



-u x ) x 2a± \u x J x ai 
In the second case we get 

D b 2 x b 3 2 , 
P = — u -\ u u 
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which corresponds to the following integrable hyperbolic partial differential 
equation: 



J. 37) u tx = ^u + -^(u 3 



a\ 6a 



I XX 



In conclusion, both equations (8.36) and (8.37) are members of the integrable 
hierarchy associated to the Lenard-Magri scheme of S-type (iv). 

Case 4: h\ = 0, b 2 0, b 3 = 0, a x # 0, a 2 = 0, a 3 ^ 

This case gives, to the right, the Lenard-Magri scheme listed as case (vii) in 
Section 8.3, while, after exchanging the roles of H and K we get, to the left, 
the "blocked" Lenard-Magri scheme listed as case 5 in Section 8.5. Hence, 
overall, we get the following scheme: 

K r +ai 3 « H r, K (Y, H i K c-{u') 2 H 



b 2 2b 2 



2bo 



Trying to solve naively for P we get 

p = ±a 13 b 2 d- 1 e ±ai3U . 

This means that the following hyperbolic partial differential equation is a 
member of the integrable hierarchy associated to the Lenard-magri scheme 
of S-type (vii): 

(8.38) u tx = +a 13 b 2 e ±ai " u . 

As expected, equation (8.38) is obtained by (8.35) letting 63 = 0. 



Case 5: 61 = 0, b 2 = 0, 63 j= 0, a x 0, a 2 ¥= 0, a 3 = 

This case gives, to the right, the Lenard-Magri scheme listed as case (x) in 
Section 8.3, while, after exchanging the roles of H and K we get, to the left, 
the "blocked" Lenard-Magri scheme listed as case 3 in Section 8.5. Hence, 
overall, we get the following scheme: 









► 0^ 




H 


ai 


v!" x 3ai 


K) 2 


o 2 1 




b 3 


(W) 3 ' ^3 




f 2b 3 {u'Y 



K r 1 H 
^ * 26 3 It' 
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Trying to solve naively for P in the above scheme, we get the following 
expression 

p = b 3 u'd- 1 (e ±ai2X u') , 
and the associated non-evolutionary partial differential equation is 

(8.39) (H) = b 3 e ±a ^ x u x . 

\U X / x 

In conclusion, equation (8.39) is a member of the integrable hierarchy asso- 
ciated to the Lenard-Magri scheme of S-type (x). Note that this equation is 
obtained letting 62 = in equation (8.34). 



Conclusion 

After rescaling the variables u, x and t, or replacing x by x+ const., or u 
by u+ const., in equations (8.34)-(8.39), we conclude that the following are 
all the integrable non-evolutionary partial differential equations which are 
members of some integrable hierarchy of bi-Hamiltonian equations, with H 
and K as in (8.1): 

(8.40) utx = e u - ae~ u + e(e u - ae~ u ) xx , 

' e x — ae~" 



.41) (*) = (e x - ae-*)u x + e 

\ U x J x 

•42) u tx = u + (u 3 ) 



x / x v U x / x 

3^ 

XX 5 

2 , 



(8.43) (^) =(-) +xu x , 

V U x / x V U x / x 

where a and e are or 1. 

Recall that the case when e = equation (8.40) is the Liouville equation 
when a = 0, and the sinh-Gordon equation when a = 1, cf. [Dor93]. 
Equation (8.41) (respectively (8.43)) can be obtained from equation (8.40) 
(resp. (8.42)) by the hodograph transformation u *-* x, x >— > —u. Equation 
(8.42) is called the "short pulse equation", [SW02], and its integrability was 
proved in [SS04]. Equations (8.40) with e = 1 was studied in [Fok95]. 



9 KN type integrable systems 

In this section V is a field of differential functions in u, and, as usual, we 
assume that V contains all the functions that we encounter in our compu- 
tations. 
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Recall from Example 4.14 that the following is a pair of compatible non- 
local Poisson structures: 

L\ = u'd^ 1 o u' (Sokolov) , L2 = d^ 1 a u'd^ 1 o u'd^ 1 (Dorfman) . 

We consider two non-local Poisson structures H and K which are linear 
combinations of L\ and L2: H = a\L\ + (I2L2 and K = b\L\ + 62-^2- As 
we have seen in the example of Liouville type integrable systems, discussed 
in Section 8, integrable hierarchies associated to Lenard-Magri schemes of 
C type are usually not very interesting (cf. Sections 8.4 and 8.5). Hence, in 
this section, we will only consider integrable Lenard-Magri schemes of S-type 
(in the terminology of Section 7.5), which is possible only when the order 
of the pseudodifferential operator H is greater than the order of K, namely 
when ai 7^ and b\ = 0. Therefore, we consider the following compatible 
pair of non-local structures: 

(9.1) h = u'r 1 on' + or 1 oii'r 1 ott'r 1 , K = r 1 oii'r 1 oM , r 1 , 

with a e C. We want to discuss the integrability of the corresponding Lenard- 
Magri scheme. 

9.1 Preliminary computations 

Note that K is the inverse of a differential operator, hence its minimal 
fractional decomposition is K = ID -1 , where 



We next find a minimal fractional decomposition for H. It is given by the 
following 

Lemma 9.1. For every a e C, we have H = AB^ 1 , where 



(9.2) 



D = do —do—d. 

11' ni' 



u u 



(9.3) 




Here and further, we have, recalling (9.2) 



(9.4) 
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The above fractional decomposition is minimal only for a =£ 0. For a = 0, 
the minimal fractional decomposition for H is H = where 

(9.5) <? = 2_doi. 

w w 

Proof. We need to prove that AB -1 = S^ 1 + aD^ 1 . By looking at the 
coefficient of a in AB^ 1 , we get 

1 dl do —Bo ^-do — l—d I = d- 1 ou'd- 1 ou'd- 1 = D- 1 . 



D{u') V v! u> D(u') 
Letting a = in AB^ 1 , we have 



u' \u'J J D[u') J V it' u' U(it') 



-u'^ 1 o D(u')dr 1 u'dr 1 o u'^ 1 = u'^ 1 ou'+f u'^ 1 o — 

\ it 

_ u "5-i + VVcT 1 - m'cT 1 o D^-Vd" 1 ) o u'r 1 . 

In the last identity we used the Leibniz rule for d: d o / = fd + /'. To 
conclude the proof, we need to check that the expression in parenthesis in 
the RHS is zero: 



(9.6) u'd~ l o — - u"d- x + { —\d- l u'd- 1 - u'd- 1 o DOu'^u'd- 1 = o . 

This identity is obtained applying repeatedly the commutation relation (/ e 
V), 

(9.7) d- 1 o / = /a- 1 - cT 1 o , 

which is a consequence of the Leibniz rule for d, and using the expression 
(9.4) for D{u'). □ 

In order to apply successfully the Lenard-Magri scheme of integrability 
we need to compute the kernel of the operator B. 
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Lemma 9.2. The kernel of the operator B in (9.3) is a ^-dimensional vector 
space over C, spanned by 



/l = l, h = ~, 
u 2 /u"y „ u" 





, /3 = 








"l 

W 






Ku'J 



Lb I Lai \ _ LL I 

U = — — -2u— + 2u'. 
w \u' J w 

Proof. It is immediate to check that all the elements fi are indeed in the 
kernel of B. On the other hand, since B has order 4, its kernel has dimension 
at most 4. □ 

9.2 Applying the Lemard-Magri scheme for a^O 

According to the Lenard-Magri scheme of integrability, starting with J/i-i = 
JO, we need to find sequences {\h n }^ =0 and {P n }^L solving the recursion 
relations (7.23). 

Since C = 1, in order to find solutions of the scheme (7.23) up to N = 3, 
we need to find elements F n , h n , P n e V, n = 0, . . . , 3, such that 

5h n —i 5h n 
PPn ? j Pn = , — = DP n , 
ou ou 

for all n = 0,1,2, 3 (we let, as usual, J/i-i = JO). Recalling the expressions 
(9.2) and (9.3) of A, B, D, and using Lemma 9.2, it is a straightforward but 
lengthy calculation to find solutions: 

F ° = : T = — (-) . P o = l, fa = ^, 

SO, 



k 






)'. 


a 

k 


au' 
= —{ 


KU '„ J 
'U N 


J u" 


a 

k 


an' 
= ~,\ 


v vl > 

'u" ^ 


' u> ' 


a 


au' 


\ u' > 


au' 



2 , 
—u 



F s = -fi = -1 



Po 


= 1 , 


fa 


Pi 


= u , 


fa 


P-2 


= u\ 


fa 


P, 


= w , 


fa 



So, 

i r /u"\2 



2. Vn 



Hence, we get the following Lenard-Magri scheme 

H 1 K r „ K r n H 2 K r n H 



JO 1 JO it JO A u 2 *±+ JO 
( 9 - 8 ) if L/A 2 x 



u 

2' 



We next prove that the scheme (9.8) can be extended indefinitely, pos- 
sibly going to a normal extension V of V. According to Theorem 7.15 and 



109 



Remark 7.10, this is the case, provided that the orthogonality conditions 
(7.16) hold. Since C = 1, the first condition in (7.16) is trivial. As for 

the second orthogonality condition, let <p e ( Span c {Po> -Pi) P2, P3}) ■ Since 
p _L Pq, we have that ip = p>' 1} for some p\ e V. Since 99 _L Pi, we have 

that pi = jjr, for some p>2 £ V. Since p L P2, we have that p>2 = -^f, for 

some 993 e V. And, finally, since 99 _L P3, we have that p% = , for some 
V?4 £ V. In conclusion, ip = Bp A , proving the second orthogonality condition 
(7.16). 

We compute explicitly the next element P 4 in the Lenard-Magri scheme, 
which gives the first non-trivial equation of the corresponding bi-Hamiltonian 
hierarchy. For this, we need to solve, for Fa,Pa e V, the following equations 

BF A = -P- = D(u) , P 4 = AF 4 . 
du 

The general solution is: 



u"\' l/^"\ 2 , w a 2 „ a 3 a 4 

+ a "«l + — /2 + — / 3 + — /4, 

2\u'y a a a 



where Qj, i = 1, ... ,4, are arbitrary constants. Hence, the first non-trivial 
integrable equation in the hierarchy has the form: 

du „ ,„ 3 (u") 2 , 9 

(9.9) — = Pa = u — j- + a\u + a 2 + a 3 u + a 4 ?r . 

dt 2 u 

In order to prove that equation (9.9) is indeed integrable, we are left to 
prove that the sequences {\h n } ne i + and {P n }neZ + are linearly independent. 
For this, we use Lemma 7.19. 

Since dord(-D(«')) = 4, we have dord(A) = 6, dord(£?) = 7, dord(C) = 
— 00 and dord(D) = 3. Moreover, \H\ = —1 and \K\ = —3. Hence, the RHS 
of inequality (7.20) is 4. 

Next, we compute the differential order of the next element P5 in the 
Lenard-Magri scheme. It is obtained by solving, for £ 4 = ^^,^5,^5 e V, 
the following equations: 

(9.10) £ 4 = DP A , BF 5 = £ 4 , P 5 = AF 5 . 
From the first equation in (9.10) we get 

£4 = d^d-^u'" + p) , 
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where p e V has dovd(p) = 2. Hence, the second equation in (9.10) gives 



— -dF 5 = u"' + Pl , 
D(u) 

where dord(pi — p) = 0. In particular, P5 has differential order less than or 
equal to 3. It follows by the third equation in (9.10) that 



Hence, = 1, and = for every n > 5. In particular, dord(Ps) = 5. 

According to Lemma 7.19, since we have dord(Ps) = 5 > 4, we obtain: 
dord (^f ) = 2n — 2, and dord(P n ) = 2n — 5, for every n ^ 3. In particular, 
all the elements {S l h n } nG i + and {P n }neZ + are linearly independent. As a 
consequence, every equation of the hierarchy Jr- = P n , n e Z+, including 
equation (9.9), is integrable of S-type. 

Note that, since the kernels of B* and D* have non-zero intersections, 
we cannot conclude that [P m ,P n ] is zero for every m,n e Z + . In fact, we 
have [P ,Px] = P , [P ,P 2 ] = 2Pi, [Pi,P 2 ] = P% and Ker(P*) nKer(Z)*) = 
Span{Po, Pi, P2} (which is isomorphic to SI2), in complete agreement with 
our Theorem 7.15 (and in disagreement, for example, with [01v93, Thm.5.36] 
and [Bla98, Thm.3.12]). 

When all constants ctj are equal to zero, equation (9.9) is usually called 
the Schwarz KdV equation, see e.g. [MS12] (in [Dor93] it is called the 
Krichever-Novikov (KN) equation, since it is a degeneration of the KN equa- 
tion). As explained in [MS12], equation (9.9) can be reduced to the Schwarz 
KdV equation by some point transformation. 

Remark 9.3. By Remark 7.23, all £ n ,'s and P n 's constructed in this section 
have coordinates in V = ¥[u, u'- 1 ,u", u'", . . . ]. By Example 4.5, this algebra 
is contained in a normal extension V = V[logu'], and all conserved densities 
h n 's can be chosen in V. 

9.3 The case a = 



In the case when a = all the computations are much easier. Since A = C 
1, the recursive ( 
to the equations 



1, the recursive conditions J/i n -i <r ~-^ Pn ~ * \h n , n e Z + , are equivalent 



BP = ^-J- - = DP 

ou ou 
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It is easy to find the first few steps of the Lenard-Magri scheme: 
(9.11) 

for arbitrary «i e C. 

As before, the scheme (9.11) can be extended indefinitely. Indeed, since 
C = 1, the first orthogonality condition in (7.16) is trivial, while the second 
one holds since Po 1 =ImB. 

Moreover, in this case dord(A) = dord(C) = — oo, dord(P) = 2, and 
dord(D) = 3, so the RHS of inequality (7.20) is 0. Since dord(P ) = 1 > 0, 
we can apply Lemma 7.19 to deduce that all the elements {$/i n }neZ + and 
{P n }neZ + are linearly independent. 

In conclusion, every equation of the hierarchy j^- = P n , n e Z + , is 

integrable of S-type. Note that the first non-trivial equation is ^ = Pi, 
which is the same as equation (9.9) with a<i = 03 = 04 = 0. Note also 
that, since Keri?* n KerD* = in this case, we have [P m ,P n ] = for all 
m,n e Z + . 

9.4 One step back 

As we did in the example of Liouville type integrable systems, we can ask 
whether the Lenard-Magri schemes (9.8) and (9.11) can be continued to the 
left. This amounts to finding P n e V and ^h n e V/dV, with n sj —1, such 
that 

(9.12) ...if/i^i^ JL^O 

We consider separately the cases a and a = 0. When a 0, the 
conditions (9.12) give the following equations for P_i, F and £_2 = ^j^-'- 

(9.13) DF-i = , AF = P_i , £_ 2 = PF , 

where A, P, and P are as in (9.2) and (9.3). All solutions P—\ of the first 
equation in (9.13) are 

P-i = Co + C\U + C2U 2 , 

with co,ci,C2 £ C. Next, we want to find all solutions F of the second 
equation in (9.13). Applying , with n > 4, to both sides of the equation 
AF = P_i we immediately get that dord(P) 3 and dF = fD(u'), with 
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dord(/) < 1. Hence, the second equation in (9.13) can be rewritten as the 
following system of equations, 

/ 9 14 ) d 2 f - 2^-5/ + f + af - u'F = c + c x u + c 2 u 2 , 

dF = fb(u') , 

for F, f e V with dord(i ? ) ^ 3 and dord(/) ^ 1. Applying to both sides 
of the first equation in (9.14) and to both sides of the second equation 
in (9.14), we get 

dF f ^L_ a 
du 1 " ~ J^f ' du 7 ~ ' 

Hence, dord(/) ^ 0. Next, applying to the first equation in (9.14) and 
-^3j to the second equation in (9.14), we get 

dF u" 1 _ df , 

Hence, / is a function of u only. Using the above result, we can rewrite the 
second equation in (9.14), after integrating by parts twice, as 

dF = d(l(-)' - d -l- + ^Lv!) - —(u'f . 
V u' V u' / du v! du 2 J dvfi 

In particular, it must be ^-^(n') 2 e dV, which is possible only if 44 i u ') 2 = 0, 
see [BDSK09]. In conclusion, / must be a quadratic polynomial in u with 

constant coefficients, and F = 77 (-77-) — §^^7 + §^ u '+ const. Plugging 
these results back into equation (9.14) we finally get that 

dF co ci c 2 2 
- / = 1 u H u . 



D{u') a a a 

Hence, the third equation in (9.13) gives £_2 = 0. In conclusion, in this 
case the "dual" Lenard-Magri sequence, obtained by exchanging the roles 
of H and K, is of finite type, namely it repeats itself with ^h n e Ker (4-) 
and P n e Ker(L>) for every n ^ — 1, and we don't get any new interesting 
integrals of motion or equations. 

Next, we consider the case a = 0. In this case, relations (9.12) give the 
following equations for P-i, P-2, and \h-2'- 

Sh 

(9.15) DF-i = , DF_ 2 = — - = , 

du 
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where S, and D are as in (9.2) and (9.5). As before, P-\ = Co + c\u + C2U 2 , 
with cq,ci,C2 e C. Hence, the second equation in (9.15) reads 

, ft1fi . (l(dP- 2 \'\' 1 fC + C 1 U + C 2 U 2 \' 

(9-16) 

For every n e Z+, we have the identity 

1 /u n \' 1 ii n nu 71 ^ 1 
v/\vf) ~ 2 '(u') 2 + 2 u' 

It follows that the RHS of (9.16) cannot be a total derivative unless ci = 
C2 = (cf. [BDSK09]). Moreover, if c\ = c 2 = equation (9.16) reduces to 

dP- 2 \' C 

— ) = TT7 + const -« > 

which, for the same reason as before, has no solutions unless c$ = 0. In 
conclusion, for every non-zero P_i 6 Ker D, the "dual" Lenard-Magri scheme 
is blocked at P_2- In this case, as we saw in Section 8.8, we obtain integrable 
PDE's which are not of evolutionary type. 

In particular, for (01,02) ^ (0,0) we get the following non-evolutionary 
integrable PDE: 



J_(^tx\ \ 1 / cp + Q\u + c 2 u 



(9.17) 

v ^ ^X X / x ^ x % 

while for Co = 1 and c\ = C2 = 0, we obtain the following integrable equation 
(9-18) (^) =^+jn x , 

V U X / X- 2^ 

where 7 is a constant. Note that, if we apply the differential substitution 
v = logu' to equation (8.40) with e = 0, we get equation (9.18). 

10 NLS type integrable systems 

Recall from Example 4.17 that the following is a triple of compatible non- 
local Poisson structures in two differential variables u, v: 

/ — 1 \ r / vd~ l ov —vd^ou 
L x = dl , L 2 = , L 3 = x j 

y 1 y y — w o « uo o u 

We want to consider two non-local Poisson structures H and A' which are 
linear combinations of them: H = a\L\ +a 2 L 2 +a^L^ and K = b\L\ + b 2 L 2 + 
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b 3 L 3 , where ctj's and frj's are constants. As in the previous section, we are 
only interested in integrable Lenard-Magri schemes of S-type. In particular, 
we assume that the order of the pseudodifferential operator H is greater 
that the order of K, and so we consider only the case when b\ = 0. Note 
that when b 2 = 0, we get K = 63L3, which is a degenerate pseudodifferential 
operator. In this case, we cannot apply Theorem 7.15 and we do not know 
how to prove integr ability. Hence, we assume that b 2 = 1- And, since we 
want that the order of H is greater that the order of K, we assume also that 
a\ = 1. 

In conclusion, we consider the following compatible pair of non-local 
Hamiltonian structures: 



(10.1) 



TT / — 1 \ / vd 1 o v —vd 1 ou 

H = dl + o 2 . n + °3 a-l a-l 

y 1 J \ —ud o v uo o u 

-1 \ , / #'o!) —vd~ 1 ou 
K = 1 1 n + M p-l fl-i 

10/ V —7X0 Ofl UO O U 



Note that if a 3 = 63 = 0, the above pair is such that H is "local" 
differential operator, and K is invertible. In this case, the Lenard-Magri 
recursion relations give J/i-i = and H Sh ?~ 1 = for every n ^ 0, hence 

= for every n. Therefore, in this case, the corresponding Lenard-Magri 
scheme is of finite type, and we don't get any integrable system. Hence, we 
assume that (03,63) ^ (0,0). 

Next, we need to find minimal fractional decompositions for H and K. 
This is given by the following 

Lemma 10.1. We have the following minimal fractional decomposition for 
the operator L3: 

vd^ 1 o v —vd^ 1 o u \ ( —uv \ / 1 
-ud^ 1 o v ud^ 1 o u J \0 u 2 J \ 7; u 

The rational matrix pseudodifferential operator H admits the fractional de- 
composition H = AB" 1 given by 

(10.2) A=f/-< -c*go«-asut,\ / 1 

\do^ + a 2 do±dou + a 3 u 2 J \% \dou 

which is minimal for a% =£ 0, while, for 03 = 0, H is a matrix differen- 
tial operator. The rational matrix pseudodifferential operator K admits the 
fractional decomposition K = CB^ 1 with B as in (10.2) and 



(10.3) t - 



-l —dou-bzuv 
1 b 3 u< 
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This decomposition is minimal for 63 ^ ; while, for 63 = 0, K = L 2 is an 
invertible matrix. 

Proof. Straightforward. □ 

As usual, in order to apply the Lenard-Magri scheme it is convenient to 
find the kernels of the operators B and C: 

KerB = C ( ° ) , KerC = c( ~Y 



We next compute the first few steps in the Lenard-Magri scheme. We 

have the following H and .fT-associations: $0 *—> Pq * — * $/io <h ~~~ > Pi < ~~~ > 

J/ii P2, where (a e C) 
(10.4) 

^0 = aa 3 ( "J ) , = ^(r 2 1 ^ 



[u + v 



Pi 



1 



P 



U'-a 2 V \ f f/ , 02, 2 . 2n , & 3, 2 , 2n2A 

+ 2a 2 u' -a|^+| (u(u 2 + v 2 ))' + 2a=°&v(u 2 + v 2 ) 



2 \ -u" + 2a 2 v' + a 2 2 u + ^ (v{u 2 + v 2 ))' - ^^ u {u 2 + v 2 ) J ' 
Indeed, we have P = AF , BF = 0, for F = a ( x ) . We have P = CFi, 



tjj^ = BF\, for F\ = ( ^ j , where a, (3 e C are chosen so that 003-/363 = 1 
(we can always do so, since, by assumption, (03,63) # (0,0)). We have 



Pi = AF 2 , = BF 2 , for P 2 = ( JJ ) ■ We have P = CP 3 , ^ = BF 3 , 



for P 3 = ^ + ^Xi^t +V2) j. And, finally, we have P 2 = AF„ 
= PP 4 , for P 4 = P 3 . 

Next, we check that the orthogonality conditions (7.16) hold for N = 0. 
We have F = f ^ J e if and only if J(it/ + vg) = 0, namely if 

f = — -q + — , for some h 6 V. But in this case 

c( 9 + b f 



- V u9+^\ =r ( 9 + huh 
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proving the first orthogonality condition (7.16). As for the first orthogonality 
condition, if 03 = there is nothing to prove since H is a matrix differential 
operator (i.e. the denominator is 1 in its minimal fractional decomposition). 

If 03 7^ 0, we can choose a = and we have F = ^ ^ j e if and only 
if £(— vf + ug) = 0, namely if g = —f + for some h e V. But in this case 

f =( E f f + hL)= B (i ) eimfl . 

V u " u / \ u / 

proving the second orthogonality condition (7.16). Therefore, by Theorem 
7.15 and Remark 7.10, we deduce that the elements (10.4) can be extended, 
possibly going to a normal extension V of V, to infinite sequences {^h n } ne z + , 

{Pn}neZ+, Such that J/in-1 «— » Pn < * P*n- 

Finally, we have = 1, |K| = 0, dord(A) = 2, dord(S) = dord(C) = 
dord(D) = 1, and dord(P 2 ) = 2. Hence, the inequality (7.20) holds. There- 
fore, by Lemma 7.19 we have dord(P n ) = dord(^J-) = n for every n e Z+. 
In particular, all the elements $/i n 's and P n 's are linearly independent. 

In conclusion, each equation of the hierarchy = P n is integrable, and 
the local functionals $ft, n 's are their integrals of motion. The first "non- 
trivial" equation of this hierarchy is for n = 2. Letting a 2 = 0, 03 = 2a, and 
63 = 2/3, it takes the form 

v" + av{u 2 + v 2 ) + (3{u{u 2 + v 2 ))' 
-u" - au{u 2 + v 2 ) + (3{v{u 2 + v 2 ))' 

If we view u and v as real valued functions, and we consider the complex 
valued function ip = u + iv, the system (10.5) can be written as the following 
PDE: 

which, for (3 = 0, is the well-known Non-Linear Schroedinger equation (see 
e.g. [TF86, Dor93, BDSK09]). The case /3 ^ has been studied by many 
authors as well, see [KN78, CLL79, WKI79, CC87]. 

It is not difficult to show that, when going back, the Lenard-Magri 
scheme is "blocked" at $/i-2 when 03 = 0, and it is of finite type when 
03 # 0. Hence, we don't get any non-evolutionary PDE in this case. 

Remark 10.2. By Remark 7.23, all £ n 's and P n 's constructed above have 
coordinates in V u = ¥[u^- 1 ,v,vf,v',u",v",...'\. Moreover, using a different 



(10.5) 



da 



v 
~dt 
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fractional decomposition, this time over V v = F[ii, v- 1 , v! , v' , u" , v", . . . ], we 
can show that all coordinates of the £ n 's and P n 's lie in V v , hence they actu- 
ally lie in the algebra of differential polynomials V = F[u, v, u 1 , v',u", v", ... ]. 
This is a normal algebra of differential functions, therefore all conserved den- 
sities /i n 's can be chosen in V. 
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